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1.0  SUMMARY 


1.1  Contract  Description 

This  contract  studies  unique  optical  computing  concepts  which  use  nonlinear 
optical  phenomena  to  perform  matrix  multiplication  and  to  provide  reconf igurable  opti¬ 
cal  interconnection.  The  study  focuses  on  the  use  of  real-time  holography  in  nonlinear 
media  such  as  photorefractive  crystals  for  optical  computing. 

1.2  Scientific  Problem 

By  incorporating  the  parallel  nature  of  optics  in  nonlinear  media,  it  is  possible 
to  perform  parallel  matrix  multiplication  using  four-wave  mixing.  In  addition,  the 
dynamic  holography  in  nonlinear  optical  media  provides  a  natural  candidate  for  the 
reconf  Igurable  interconnection.  The  general  problem  in  this  program  is  to  generate  and 
investigate  new  concepts  which  use  these  nonlinear  optical  phenomena  for  optical 
computing. 

Specifically,  this  program  investigates  experimentally  and  theoretically  the 
multiplication  of  rviatrices  using  optical  four-wave  mixing  in  nonlinear  media,  and  the 
possibility  of  using  such  matrix  multiplication  and  wave  mixing  for  reconfigurable 
interconnection. 

1.3  Progress  Summary 

There  are  several  areas  of  significant  progress  achieved  under  this  contract 
that  are  directly  related  to  the  development  of  optical  matrix  multiplication  and  recon¬ 
figurable  optical  interconnection.  These  include: 

1.  First  experimental  demonstration  of  parallel  matrix-vector  multiplier 
using  optical  four-wave  mixing  in  a  barium  litanate  (BaTiO^)  crystal. 

2.  First  experimental  demonstration  of  summation  process  inside  the  non¬ 
linear  media  in  matrix-vector  multiplications. 

3.  Experimental  demonstration  of  2  x  2  matrix-matrix  multiplication  using 
optical  four-wave  mixing  in  a  BaTiO^  crystal. 
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4.  Experimental  demonstration  of  matrix-matrix  multiplication  using  color 
multiplexing. 

5.  Experimental  demonstration  of  matrix-matrix  multiplication  using  convo¬ 
lution. 

6.  Experimental  demonstration  of  matrix-vector  multiplication  using  a 

spatial  light  modulator  and  a  phase  conjugator. 

7.  Development  and  experimental  demonstration  of  a  new  concept  of 

reconf igurable  optical  interconnection  using  photorefractive  holograms. 

8.  Experimental  demonstration  of  high-eff  iciency  interconnection  using 

BaTiO^  crystals. 

9.  Development  of  matched  amplification  at  the  Fourier  plane  to  achieve 
maximum  efficiency. 

10.  Experimental  demonstration  of  high  efficiency  in  photorefractive  optical 
interconnection  using  matched  amplification. 

1  1.  Experimental  demonstration  of  reconfigurability  using  a  liquid  crystal  TV 
in  conjunction  with  a  photorefractive  barium  titanate  crystal. 

1.4  Publications  and  Presentations 

1.  "Reconf igurable  Interconnection  Using  Photorefractive  Holograms,"  Proc. 
SPIE,  Vol.  1151,  Paper  No.  03,  to  be  presented  at  SPIE's  Annual  International 
Sumposium  on  Optical  and  Optoelectronic  Applied  Science  and  Engineering 
(August,  1989,  San  Diego,  CA). 

2.  "Energy  Efficiency  of  Optical  Interconnections  Using  Photorefractive 
Holograms,"  submitted  to  Appl.  Opt.  (1989). 

3.  "On  Image  Amplification  by  Tw'o-^^ave  Mixing  in  Photorefractive  Crystals," 
submitted  to  Appl.  Opt.  (1989). 

4.  "Photorefractive  Nonlinear  Optics  and  Optical  Computing,"  Opt.  Eng.  28(4), 
328  (1989). 

5.  "Energy  Efficiency  of  Optical  Interconnection  Using  Photorefractive  Dynamic 
Holograms,"  Optical  Computing.  1989  Technical  Digest  Series,  Vol.  9  (Optical 
Society  of  America,  \tashington,  D.C.),  pp.  128. 
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b.  "Fnerg\  Ftficient  Optical  Iniercontieciion  Using  Dynamic  Holograms  in 

Photorelractive  \tedia,"  OSA  Annual  Meeting.  1988  Technical  Digest  Series, 
Vol.  1  !  (Optical  Society  ol  America,  Washington,  D.C.,  1988),  pp.  178. 

7.  "Optical  Interconnection  Using  Photorelractive  Dynamic  Holograms,"  Appl. 
Opt.  2093  (1988). 

8.  "Optical  Matrix-\'ector  Multiplication  Using  a  Spatial  Light  Modulator  and  a 
Phase  Conjugator,"  Spatial  Light  Modulators  and  Applications,  1988  Technical 
Digest  t^erios,  Vol.  8  (Optical  Society  of  America,  Washington,  D.C.,  1988),  pp. 
208. 

9.  "Optica!  Matrix-Matrix  Multiplication  Using  Multicolor  t'our-Wave  Mixing," 
Proc.  SPIF,  Vol.  881,  paper  38,  presented  at  OU/LASU'88  (3an.  10-17,  1988,  Los 
Angeles,  CA). 

1C.  "Optical  5(atrix-\ ector  Multiplication  via  i'our  Waxe  Mixing  m  Photorefrav- 

tive  ^(edla,"  0})t.  Lett.  ^2,  138  (1987);  Opt.  Lett.  [2,  373  (1987). 

11.  "Optical  Matnx-\ector  Multiplication  Using  Four-Wave  Mixing,"  Paper  MM5, 

OSA  Annual  Meeting  (Oct.  1986,  Seattle,  W  A),  3.  Opt.  hoe.  Am.  ^  (13),  ib 
(198b). 

1 . 5  FMlcnt  Discloc'ircs 

"Multicolor  Matrix-Matrix  Multiplier  v'.ith  Paraiieii.on,"  M.  Khoshnevisai.,  A.L.T. 

Chiou  and  P.  Veh,  Rockv,ell  Patera  Disclosure  87SC37. 

"Matrix-Vector  Multiplier  using  Photorefractive  Phase  Conjugators,"  P.  3  eti,  .A.L.T. 

Chiou  and  \U  Khoshnevisan,  Rockwell  Patent  Disclosure  87SC96. 

"Optical  Matrix-Matrix  Multiplier  with  \'^  Parallelism  by  Spatial  Convolution  via  Four- 

Wave  Mixirig,"  A.L.T.  Chiou,  Rockwell  Patent  Disclosure  87SCb3. 

"t^econfigurabie  Optica!  Interconnect  using  Dynamic  Holograms,"  P.  Veh,  Rockweli 

Patent  Disclosure  87SC5L. 
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2.0  TFCHNICAL  DISCUSSION 


The  technical  problem  addressed  in  this  contract  is  the  study  of  optical  matrix 

multipiicaMon  and  recent  igurabie  optical  interconnection.  By  incorporating  the  parallel 

1  2 

naij’-e  o!  optici>  m  nonlinear  media  such  as  photorefractive  crystals,  ’  it  is  possible  to 

perform  parallel  matr  ix  multiplication  using  four-wave  mixing.  In  addition,  the  dynamic 

holographs  in  noniinear  optical  media  provides  a  natural  candidate  for  the  reconf igurat le 

3  U 

optical  interconnec  tion. 

hpecilical!) ,  ti.is  program  m  aimed  at  the  uiilixalioii  of  nonlinear  optica! 
ptmiiomeria  bijc  h  as  phase  conjugat  lori,^  '  lw,o-v.ave  mixing, and  tour-wav  e  mixing^  ^ 
in  nonlinea’  media  ’  to  perform  tne  rni  tnx-vector  and  matrix-matrix  multiplications 
afic:  rec o'i :  onrarable  optical  interconnection.  Photoref''a<' tiv e  crvsials  such  as  BaTiO^ 
tmc:  stror,ti'.!!i.  ija.rium  niohate  (SBN)  we-e  used  as  the  norilmear  media. 

Ciei.erallv  speiiKing,  tour-v%ave  mixing  is  a  nonlinear  optical  pri.vcess  iii  winch 
three  ir.put  waves  mix  1^  vieid  a  fourtf,  wave.  In  phase-matched  four-wave  mixing,  the 
three  input  waves  consist  ot  two  i  ouriterpropagaiing  pump  waves,  Ej  and  r.2"  an  arbi- 
trar;>  n'-otte  wavcc  1. Al:  thoee  couple  through  the  third-order  suscepubi  1 1 1  , ,  >  to 

vi.hd  a  to.,’  t'  vv.-.ve.  1.^^.  w  ‘  a  n  is  proporiior.tii  to  the  produr  i  of  Ej,L2  ‘iiid  complex  conju- 
.1-11 


hijcii  a  l^tur  Wave  mixma  e  ai  be  understood  in  terms  of  the  t  ocording  arid  readout  pr  oces¬ 
ses  wiiuA;  occur  iri  hoiogr ani'v .  li:  trie  nonlinear  media,  one  of  the  pump  waves  and  ttie 
probe  Wave  form  an  mterlereiK  e  pattern  which,  in  turn,  induces  an  index  grating.  The 
oiner  pumij  wdv>.‘  is  diffracted  from  this  grating  and  generates  tfte  fourtli  wave.  Tfie 
toruiatior:  of  the  gra'ing  (or  hoUvgr.tm)  and  readout  processes  take  place  at  the  same 
l  m  Inns,  f  our  x  av  e  m  l«  mg  is  sometimes  referred  to  as  a  rea  1-t  ime  ho  log  r  a  pin  . 

liv  usuig  four  wave  mixing  in  nonlinear  media,  multiplication  o!  signals  can 
take  plcK.e  in  a  subpic  osecoisi  lime-  scale.  In  addition,  it  we  use  the  paiallel  nature  of 
i.'ptira>  wi.'',.  eat.  h  wave  t  ,i'i  i  arr,  suatial  iri  for  ma  t  lof  i  fc^r  tiic*  purjx'se  cn*  image*  [xrc’icess- 
ine  uang  the  iwe  t  r  a '  c.v  er  i  ■  d  imenaor  ;s.  Sut.ii  an  ap[;lualK'n  ot  io‘ir  wave  mixing  to 
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real-tune  iinaf  processing  was  recently  reported  using  BSO  crystals. ^ ^  In  what  follows, 
v.e  will  describe  some  unique  concepts  which  use  the  two  transverse  dimenoions  to  carry 
the  matrix  information  for  the  purpose  of  matrix  multiplication  and  reconfigurable 
optical  interconnection. 

2. 1  Matrix-Matrix  Multiplication 

The  matrix  multiplication^^  between  two  N  \  N  matrices  can  be  stated  as 

fol  lows: 


C  =  AE  (1) 

where 

Note  that  a  matrix  multiplication  consists  of  two  mam  operations,  a  parallel  multiplica¬ 
tion  and  an  addition. 

Referring  to  Fig.  la,  let  us  consider  a  four-wave  mixing  configuration  which  is 

1  3 

sijitablt'  for  matrix  multiplication.  Hearn  1  and  Beam  2  contain  the  information  about 
the  two  matrices  A(x,?,)  and  respectively.  Beam  1  is  propagating  along  the  y-axis, 

and  Beam  2  is  propagating  along  the  x-axis.  These  directions  are  chosen  for  the  sake  of 
clarity  in  introducing  the  concept.  They  are  not  the  only  directions  for  u'atrix 
multiplication.  Also,  these  matrices  can  be  either  continuous  or  discrete.  In  the  discrete 
case,  each  beam  consists  of  a  matrix  of  beamlets,  as  shown  in  Fig.  lb. 

In  the  nonlinear  medium,  these  two  matrix-carrying  beams  form  an  interfer¬ 
ence  pattern.  As  a  result  of  the  nonlinear  response  of  the  medium,  a  volume  grating  is 
formed.  This  grating  contains  information  about  the  product  of  the  matrix  elements  of 
ihiese  two  matrices,  and  can  be  written  as: 


cn  = 


A(x,z)B*(z,y)e 


+  c.c. 


(3) 
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[•ig.  1  \tatrix-niatrix  fnulliplication  using  four-wave  mixing. 

-  I^is  ihe  diMerence  of  lf>e  uavevectors  of  the  ma Irix-carrying  beams  and  cc  repre- 
complex  conjuga'iion.  The  parameter  02  ii>  the  Kerr  coefficient  and  is  proportional 


6 

CI0178TD/bje 


Rockwell  International 

Science  Center 


SC5502.fr 

to  the  third-order  susceptibility  of  the  medium.  Note  that  the  nonlinear  response  of 
the  medium  serves  the  purpose  of  parallel  multiplication. 

The  volume  grating  is  then  read  out  by  a  third  beam  which  can  simply  be  a 
plane  wave.  The  diffracted  beam  consists  of  the  integrated  contribution  from  each  part 
of  the  grating  along  the  beam  path,  and  thus  can  be  written 

C(x,y)  -  I  A(x,it)B*(z,v)dz  (4) 

where  the  integration  is  carried  out  along  the  beam  path.  Note  that  the  integration  com¬ 
pletes  the  operation  of  matrix  multiplication.  The  information  about  the  product  of 
these  two  matrices  is  nov.  impressed  on  the  transverse  spatial  diatributiori  of  the  dif¬ 
fracted  beam. 

Due  to  the  phase  matching  requirement,  the  readout  beam  must  be  incident 
along  the  directions  which  satisfy  the  Bragg  diffraction  condition  to  achieve  high  effi¬ 
ciency.  In  anisotropic  nonlinear  media,  the  polarization  states,  as  well  as  the  direction 
of  propagation,  can  be  chosen  such  that  the  largest  of  the  nonlinear  susceptibilities  is 
fully  utilized. 

The  four-wave  mixing  can  be  either  degenerate  or  nondegenerate.  In  the 
degenerate  case,  al!  the  beaiiii  have  the  same  frequency.  In  the  nondegenerate  case,  the 
frequencies  of  the  beams  can  be  slightly  different.  This  may  be  useful  for  the  purpose  of 
separating  the  diffracted  beam  from  the  undiffracted  portion. 

To  illustrate  the  information  capacity  of  such  a  nonlinear  optical  matrix 
multiplication,  let  us  consider  a  four-wave  mixing  process  using  an  Ar  ion  laser  at  48S0A 
in  a  medium  of  a  1  cm  cube.  The  grating  space  is  of  the  order  of  0.5  sm.  Thus,  10  um  x 
10  Mill  IS  enough  for  each  pixel  of  information.  In  other  words,  a  1  cm  cube  is  cable  of 
handling  1000  x  1000  matrices.  \t  ith  a  material  response  time  of  1  ns,  such  a  matrix 
multiplication  has  a  potential  data  throughput  rate  of  quadrillion  bits  per  second 
do'^  bits/s)! 
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2.2  Matrix-Vector  Multiplication 

Another  scheme  as  shown  in  Fig.  l  is  suitable  for  matrix-vector  multiplica- 

11-1^ 

tion.  Here,  as  an  example,  let  us  consider  a  discrete  case  in  which  we  need  to  carry 

out  the  multiplication  of  an  N-elernent  vector  and  an  N  x  N  matrix.  The  vector  is  fanned 
out  into  N-rows  of  identical  vectors.  These  N  x  N  small  beams  are  directed  to  a  non¬ 
linear  medium.  The  matrix  which  also  contains  N  x  N  small  beams  is  also  directed  to  the 
medium  in  such  a  way  that  each  beam  of  the  matrix  is  counterpropagaling  in  a  direction 
relative  to  the  corresponding  beam  of  the  vector.  Thus,  in  the  medium,  there  are  N  x  N 
spatially  separated  regions  which  are  pumped  each  by  a  pair  of  counterpropagating 
beams.  Nou,  N  x  N  probing  beams  are  directed  into  the  medium  in  such  a  way  that  each 
probe  beam  will  propagate  through  an  intersection  region.  The  probe  beams  will  be 
"plane  wave"  beamlets  propagating  in  parallel.  As  a  result  of  the  four-wave  mixing,  each 
probe  beam  will  generate  a  phase-conjugated  beam  which,  within  a  proportional  factor, 
can  be  written  M(i,j)a(j).  By  using  a  cylindrical  lens,  a  summation  over  j  can  be  ob¬ 
tained.  Thus,  we  have 

b(i)  =  V  M(i , j)a( j)  (5) 

J 

where  a(j)  is  the  j-th  element  of  the  vector  a  and  M(i,j)  is  the  matrix  element.  Such  a 
scheme  for  matrix-vector  multiplication  can  also  be  used  for  matrix-matrix  multiplica¬ 
tion  by  decomposing  a  matrix  into  column  vectors  and  then  multiplying  the  matrix  with 

each  of  the  column  vectors.  Using  color-multiplexing,  one  can  perform  matrix-matrix 

3  1 5 

multiplication  with  N  -parallelism. 

2.3  Reconfigurable  Optical  Interconnection 

Optical  interconnection  will  play  a  key  role  in  both  the  optical  computing  and 
VLSI  systems. ^  There  are  many  advantages  of  optical  interconnect.  These  include 
high  space-bandwidth  and  tirne-bandwidth  products,  thus,  many  independent  channels 
could  be  exploited  for  demanding  computations.  Optical  processors  are  inherently  two- 
dimensional  and  parallel.  Optical  signals  can  propagate  through  each  other  in  separate 
channels  with  essentially  no  interaction.  Optical  signals  can  also  propagate  m  parallel 
channels  without  any  interference  and  crosstalk.  In  the  VLSI  systems,  optical  inter 
connect  can  be  used  to  solve  the  problem  of  communication,  as  well  as  the  clock 
distribution. 
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Fig.  2  Matrix-vector  multiplication  via  four-wave  mixing. 

Computer-generated  holograms  (CGH)  can  be  used  for  fixed  optical  intercon¬ 
nect.  The  most  general  CGH  consists  of  a  two-dimensional  array  of  subholograms.  Each 
subhologram  is  capable  of  diffracting  a  gate  output  to  any  gate  or  combination  of  gate 
inputs.  Such  a  fixed  interconnect  pattern  is  adequate  for  many  applications,  including 
the  point  operation,  matrix  operation,  the  clock  distribution  and  a  globally  synchronous 
systolic  processor.  However,  there  are  many  situations  in  which  a  dynamic  optical  inter¬ 
connect  is  required. 

In  Fourier  transform  and  sorting,  every  element  of  an  output  array  is  affected 
by  all  the  elements  of  the  input  array,  and  conversely,  each  element  of  the  input  array 
affects  all  elements  of  the  output  array.  In  addition,  the  interconnections  change  at  dif¬ 
ferent  stages  of  the  computation.  Thus,  these  operations  require  global  and  dynamic 
interconnection  between  different  elements  of  the  input  array. 

In  image  restoration,  pattern  recognition,  and  neural  network  systems,  the 
interconnect  patterns  could  be  data-dependent,  making  it  impossible  to  foresee  the 
interconnect  requirements  at  different  stages  of  processing  without  having  foreknow¬ 
ledge  of  the  input.  The  computational  throughput  of  parallel  processor  implementing 
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these  types  of  operation  v,iU  be  critically  affected  by  the  availability  of  a  dynamic  and 
global  interconnect  network. 

In  Fourier  transform  and  sorting,  the  requirement  for  dynamic  interconnect 
can  be  avoided  by  resorting  to  a  fixed  but  global  interconnect  pattern  known  as  the  "per¬ 
fect  shuffle".  In  the  area  of  image  restoration  and  pattern  recognition,  there  is  a  mini¬ 
mum  amount  of  regularity  and  structure.  In  addition,  there  is  a  possible  data  and  time 
dependency  in  the  interconnect  requirements.  A  global  and  reconfigurable  interconnect 
network  will  be  vital  to  achieving  high  throughput  and  high  efficiency  with  parallel  proc¬ 
essors  implementing  these  operations. 

The  most  general  interconnect  system  is  one  in  which  any  gate  output  can  be 
connected  to  the  input  of  any  gate  or  combination  of  gates  (see  Fig.  3).  The  effect  of 
such  an  interconnect  can  be  repres«.nted  by  the  matrix  equation 

O  =  M  1  (6) 

w'here  I  is  a  vector  representing  the  two-dimensional  input  array,  M  is  the  matrix  repre¬ 
senting  the  interconnect,  and  O  is  a  vector  representing  the  output  array.  In  digital  opti- 
cal  computing,  ’  the  input  array  is  actually  the  gate  output  array.  Each  matrix 
element  (i,j)  is  nonzero  if,  and  only  if,  there  is  a  connection  between  pixel  j  of  the  input 


non;.in&ar  medium 


GATE 

OUTPUT 


GATE 

INPUT 


I'lg.  3  Reconfigurable  optical  interconnection. 
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array  and  pixel  i  of  the  output  array.  The  matrix-vector  multiplication  technique 
described  in  Section  2.2  can  be  used  as  an  optical  interconnect,  which  provides  both  local 
and  global  communication  between  the  gate  input  and  gate  output.  In  the  VLSI  systems, 
the  vectors  I  may  consist  of  N  laser  beams,  each  containing  a  stream  of  data  from  a 
processor,  and  the  vector  O  represents  the  output  laser  beams,  each  feeding  into  a  proc¬ 
essor.  The  main  advantage  of  using  such  nonlinear  matrix-vector  multiplication  for 
interconnection  is  that  the  matrix  M  can  be  easily  changed  for  the  purpose  of  reconfigur¬ 
ing  the  interconnect. 

Equation  (6)  can  be  generalized  such  that  both  the  input  array  1  and  the  output 
array  O  are  two-dimensional  matrices.  In  that  case,  the  most  general  interconnect 
matrix  M  must  be  four-dimensional  (i.e.,  tensor  of  rank  four).  Each  matrix  element  (ijkl) 
represents  the  connection  between  pixel  (ij)  of  the  output  array  and  pixel  (ki)  of  the  input 
array.  Because  of  the  two-dimensional  nature  of  optical  waves,  it  is  desirable  to  repre¬ 
sent  the  tensor  M  by  a  two-dimensional  array  of  two-dimensional  submatrices.  Thus,  the 
matrix-matrix  multiplication  described  in  Section  2.1  can  be  used  for  reconfigurable 
optical  interconnect  which  provides  the  most  general  interconnection  pattern.  While  this 
interconnection  scheme  allows  complete  generality,  a  price  is  paid  in  terms  of  the  space- 
bandwidth  product  requirements  on  the  nonlinear  media. 
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3.0  PROGRESS 

3.1  Matrix-Vector  Multiplication 

13  19 

In  our  earlier  demonstration  of  optical  matrix-vector  multiplication  ’  (see 
Appendix  5.1),  we  have  used  optical  phase  conjugation  via  nonlinear  four-wave  mixing  in 
a  photorefractivc-  BaTiO^  crystal  to  perform  parallel  multiplication.  The  summation 
required  to  obtain  matrix-vector  products  is  performed  subsequently  by  using  a  cylin¬ 
drical  lens  external  to  the  nonlinear  medium.  With  a  slight  modification  of  the  experi¬ 
mental  geometry,  we  have  also  successfully  demonstrated  that  both  the  parallel 
multiplication  and  the  summation  can  be  done  inside  the  nonlinear  medium  to  achieve  the 
desired  matrix-vector  multiplication  without  an  external  cylindrical  lens. 

The  top  view  of  the  new  experimental  geometry  is  illustrated  schematically  in 
Fig.  4.  The  readout  beam  in  this  case  consists  of  a  vertical  column  of  beamlets  with 
equal  intensity.  This  beam  enters  the  crystal  at  an  oblique  angle  such  that  each  beamlet 
traverses,  inside  the  crystal,  all  the  counterpropagating  pumping  beamlets  at  the  same 
horizontal  plane.  As  the  proper  elements  of  the  matrix  and  the  vector  are  encoded  in  the 
two  counterpropagating  beams  via  the  appropriate  masks,  the  phase-conjugate  output  of 


PHOTOREFRACTIVE 


Fig.  4  Optical  matrix-vector  multiplication  using  photorefractive  four-wave  mixing. 
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each  beamlet  consists  of  the  sum  of  the  product  resulting  from  each  individual  encoun¬ 
ter.  From  a  pure  geometrical  point  of  view,  the  probe  beam  (i.e.,  the  readout  beam)  can 
conveniently  be  injected  (from  c-face  of  the  crystal)  perpendicular  to  the  counterpropa- 
gating  pumping  beams. 

The  experimental  results  for  the  multiplication  of  a  4  x  4  binary  matrix  and  a 
4  X  1  binary  vector  are  shown  in  Fig.  5.  The  magnified  images  of  the  masks  used  to 
encode  the  matrix,  the  vector  and  the  probe  are  shown  in  Figs.  5a,  b  and  c,  respec¬ 
tively.  The  experimental  results  representing  the  final  matrix-vector  product  is  shown  in 
Fig.  5d.  The  relative  intensity  of  each  output  spot  representing  the  element  of  the 

product  vector  is  shown  in  Fig.  6.  Due  to  the  slow  response  time  (or  the  order  of  second 

2 

for  optical  intensity  of  the  order  of  a  few  tens  of  miliiwatt/cm  )  and  the  consequential 
sensitivity  to  environmental  changes  of  the  photorefractive  process,  the  output  intensity 
fluctuates  significantly. 


(a) 

INPUT  MATRIX 


(b) 

INPUT  VECTOR 
(EXPANDED  INTO 
MATRIX  FORM) 


(c) 

READ  OUT 
(PROBE) 
VECTOR 


(d) 

OUTPUT 
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Fig.  5  .  xperimental  results  of  optical  matrix-vector  multiplication  using  phoiore- 

fractive  four-wave  mixing. 
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Fig.  6  Output  of  the  detector  monitoring  the  intensity  of  the  optical  beam  represent¬ 
ing  the  product  vector. 
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Fig.  7  Decomposi tion  of  matrix-matrix  multiplication  into  matrix-vector  multipli¬ 
cation. 

3.2  Matrix-Matrix  Multiplication  via  Color  Multiplexing 

Vf  e  have  also  extended  our  experimental  demonstration  on  matrix-vector 
multiplication  to  matrix-matrix  multiplication  by  color  multiplexing  ^  (see 

Appendix  5.2).  The  basic  idea  is  to  decompose  the  problem  into  matrix-vector 
multiplications,  as  shovxn  in  Fig.  7,  and  carry  out  all  the  parallel  matrix-vector 
multiplications  simultaneously  by  using  color  multiplexing.  Note  that  each  of  the 
operations  shown  on  the  right-hand  side  of  the  equation  in  Fig.  7  is  in  fact  a  matrix- 
vector  multiplication.  The  basic  principle  of  color  multiplexing  used  to  encode  the 
component  vectors  with  different  colors  is  illustrated  in  Fig.  8  for  the  case  of  4  .  4 
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Fig.  8  Optical  niainx-matrix  multiplication  via  color  multiplexing. 

matrices.  Each  row  vector  of  the  matrix  is  illuminated  by  one  color,  and  all  the  color 
components  are  then  combined  by  the  prism  (angular  multiplexing)  into  a  single  row  prior 
to  further  expansion  by  anamorphic  optics  (not  shown  in  the  figure)  to  match  the  mask 
representing  the  second  matrix.  After  proper  element-by-element  multiplications  and 
summation  (summing  optics  omitted)  into  the  column,  the  resulting  multicolor  output  is 
then  demultiplexed  into  different  color  components  that  together  represent  the  final 
product. 

Using  an  argon  ion  laser  that  oscillates  simultaneously  at  five  colors  in  the 
blue-green  and  a  BaTiO^  crystal,  we  have  demonstrated  the  principle  described  above  for 
the  case  of  2  x  2  matrices.  Details  are  given  in  Appendix  5.2. 

3.3  Matrix -Matrix  Multiplication  using  Spatial  Convolution 

In  addition  to  the  matrix-matrix  multiplication  using  color-multiplexed  four- 
wave  mixing  described  in  the  previous  section,  we  have  also  successfully  demonstrated 
another  scheme  for  matrix-matrix  multiplication  using  spatial  convolution  via  four-wave 
mixing.  From  the  experimental  point  of  view,  the  key  difference  between  this  approach 
and  the  others  described  above  is  that  the  nonlinear  crystal  is  now  located  at  the  com¬ 
mon  Fourier  plane,  rather  than  the  common  image  plane,  of  the  input  matrix  masks.  The 
encoding  scheme,  as  explained  in  Appendix  5.3,  is  also  different  from  those  used  in  the 
earlier  approaches.  From  the  conceptual  point  of  view,  matrix-matrix  multiplication  in 
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full  parallelism  is  achieved  by  space-multiplexing  via  spatial  convolution  using 
degenerate  four-wave  mixing. 

3.4  Matrix-Vector  Multiplication  Using  a  Spatial  Light  Modulator  and  a  Phase 

Conjugator 

Ue  have  developed  and  demonstrated  another  new  scheme  for  optical  matrix- 

vector  multiplication  that  uses  a  phase  conjugator  (with  a  finite  storage  time)  in 

conjunction  with  a  spatial  light  modulator  (SLM)  (see  Appendix  5.4).  The  optical  system 

involved  is  relatively  simple  compared  with  the  other  approaches.  Without  any  modified- 

2 

tion,  such  a  scheme  can  also  perform  matrix-matrix  multiplication  with  N  /2  parallelism. 

Referring  to  Fig.  9,  we  use  a  SLM  to  impress  the  matrix  and  vector  informa¬ 
tion  in  sequence  to  an  input  laser  beam.  This  beam  is  directed  toward  a  phase  conjugator 
which  has  a  finite  storage  time  (a  photorefractive  barium  titanate,  for  example).  A 
cylindrical  lens  is  inserted  in  the  phase  conjugate  beam  path  to  perform  the  summation. 


COLLIMATED 
LASER  BEAM 

JJ 


VECTOR  X 

T  =  A^ 

Fig.  9  A  schematic  diagram  illustrating  the  basic  concept  of  optical 
matrix-vector  multiplication  using  a  photorefractive  phase 
conjugator  in  conjunction  with  a  spatial  light  modulator. 

The  principle  of  operation  is  as  follows:  the  SLM  first  impresses  the  matrix 
information  onto  the  input  laser  beam.  This  beam  is  then  incident  into  a  phase 
conjugator  which  stores  the  matrix  information  after  a  finite  grating  formation  time. 
When  the  matrix  information  is  removed  from  the  SLM,  say  by  turning  all  the  pixels  into 
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maximum  transmission  condition,  the  phase  conjugate  beam  which  contains  the  recon¬ 
struction  of  the  matrix  information  exists  for  a  finite  duration.  This  finite  storage  time 
depends  on  the  strength  of  the  input  (read)  beam.  During  this  time,  if  the  next  frame  of 
the  SLM  carries  the  vector  information,  parallel  multiplication  is  performed  as  the 
phase  conjugate  beam  propagates  back  through  the  SLM.  Here,  the  vector  is  represetned 
as  a  two-dimensional  array  of  N  identical  column  vectors,  where  N  is  the  dimension  of 
the  vector.  A  cylindrical  lens  in  the  output  port  is  used  to  performed  the  summation. 
The  dark  storage  time  during  which  the  matrix  information  can  be  retrieved  is 

determined  by  the  photorefracti ve  material  and  the  pumping  configuration.  It  ranges 

2 1 

from  seconds  to  microseconds. 

The  system  can  also  perform  matrix-matrix  multiplication  by  time  multi¬ 
plexing.  In  this  case,  each  column  vector  V|  (i  =  1  to  N)  v,hich  constitutes  the  second 
matrix  M2  is  sequentialU  impressed  onto  the  beam  to  multiply  with  the  first  matrix  Mj 
according  to  the  matrix  vector  multiplication  scheme  described  below.  To  avoid  the 
degradation  of  the  information  of  M[  stored  in  the  photorefractive  hologram  during  th 
readout,  it  is  necessary  to  refresh  the  holographic  memory  with  Mj  to  restore  its 
diffraction  efficiency.  This  can  be  done  by  re-impressing  M]  onto  the  beam  after  each 
readout  cycle.  Consequently,  a  total  of  2\'  clock  cycles  consisting  of  N  cycles  of  write 
and  N  cycles  of  read,  will  be  required  to  carry  out  the  inultiplication  of  two  N  ■  N 
matrices. 

Using  a  photorefractive  barium  tilanate  crystal  as  a  phase  conjugator  in 
conjunction  with  a  .  L8  magneto-optic  spatial  light  modulator  (SIGHT-MOD  SMDLSI 
from  Semetex  Corp.),  we  have  demonstrated  the  basic  principle  described  above.  The 
virtue  of  optical  phase-conjugator  ensures  that  precise  pixel  by  pixel  alignment  is 
achieved  automatically. 

3.5  Optical  Interconnection  using  Dynamic  Photorefractive  Holograms 

3.5.1  Principle  of  Operation 

During  the  first  year  of  the  program,  we  have  developed  a  new  concept  of 

22 

reconf  igurable  optical  interconnection  using  photorefractive  holograms  (see 
Appendix  5.5).  This  new  scheme  provides  an  optical  interconnection  between  an  array  of 
laser  sources  and  on  array  of  detectors  with  a  very  higti  energy  efficiency. 
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Reconfigurable  interconnection  linking  laser  arrays  and  detector  arrays  plays  a 
key  role  in  optical  computing.  Conceptually,  such  interconnection  can  be  achieved  by 
using  an  optical  matrix-vector  multiplication. 

v'  ^  M  V  (7) 

where  v  is  the  input  vector  representing  the  signals  carried  by  the  array  of  lasers,  and  v' 
IS  the  output  vector  representing  the  signals  carried  by  the  array  of  detectors.  The 
matrix  M  represents  the  interconnection  pattern.  Uhen  a  transparency  or  spatial  light 
modulator  (SLM)  is  used  as  the  interconnection  pattern,  a  large  fraction  of  energy  is 
absorbed  by  the  transparency  or  SLM.  This  energy  loss  increases  as  the  dimension  of  the 
array  increases.  Often,  the  fractional  energy  loss  can  be  as  large  as  (N-l)/N,  where  N  is 
the  dimension  of  the  array.  For  a  1000  >  1000  crossbar  suitch,  the  loss  due  to  fanout  can 
be  as  big  as  99.9%.  This  is  not  acceptable  in  high-speed  computing  because  signals  are 
passing  through  the  SLM  billion  times  per  second  ano  the  energy  loss  can  be  enormous. 

l^efernng  to  lugs.  10  and  11,  we  describe  a  new  method  of  reconfigurable  opti¬ 
cal  interconnection  which  uses  the  nonreciprocal  energy  transfer  in  photorefractive  tv.o- 
wave  mixing  to  itnprove  the  energy  efficiency.  Figure  10  describes  a  one-dimensional 
case  for  the  sake  of  clarity  iri  explaining  the  concept.  A  small  fraction  of  a  laser  beam  is 
coupled  out  of  the  beam  by  using  a  beam  splitter.  This  small  fraction  (called  probe 
beam)  is  then  expanded  by  using  a  cylindrical  lens  and  passes  through  the  SLM.  In  the 
example  shown,  the  laser  beam  is  to  connect  to  Detectors  b  and  d  as  prescribed  by  the 
SLM.  The  transmitted  beam  is  then  recombined  ixith  the  main  beam  inside  a  photore¬ 
fractive  crystal.  As  a  result  of  nonreciprocal  energy  coupling,  almost  all  the  energy  in 
the  mam  beam  is  transferred  to  the  probe  beam  which  carries  the  interconnection  pat¬ 
tern.  The  results  is  an  optical  interconnec tion  with  high  energy  efficiency.  Figure  11 
describes  the  reconfigurable  interconnection  for  laser  arrays  and  detector  arrays.  In  the 
example  (a  4  ■  4  interconnection)  shown.  Laser  1  is  to  connect  to  Detectors  b  and  c. 
Laser  2  is  to  conrn.ct  to  Detectors  a  and  d,  Laser  3  is  to  connect  to  Detectors  c  and  d. 
Laser  4  is  to  connect  to  Detector  a  and  c.  A  cylindrical  lens  is  used  to  focus  the  two- 
dimensional  array  of  beams  into  a  vector  (one-dimensional  array).  As  a  result,  Detector 
a  receives  signals  from  Lasers  2  and  4,  Detector  b  receives  signals  from  Laser  1,  Detec¬ 
tor  c  recei'  es  signals  from  Lasers  1,  3  and  4,  and  Detector  d  receives  signals  from  Lasers 


18 

CI0178TD/b)e 


Rockwell  International 

Science  Center 


SC5502.fr 


•C41«l I 

MIRROR 


Fig.  10 


A  schematic  drawing  of  a  1  ^  N  optical  inierconnei-lion  using  dynamic 
pholorefractive  holograms.  N  =  4. 


Fig.  1  1  A  schematic  dravi.  ing  of  a  N  «  N  optical  interconnection  using  dynamic 
pholorefractive  holograms.  N  =  4. 

2  and  3.  Such  a  concept  can  be  extended  to  interconnect  Nj_  lasers  with  Nq  detectors 
where  N[_  and  are  two  large  numbeis. 

The  interconnection  can  be  reconfigured  by  using  a  different  SLM  pattern, 
Tfie  energy  los^  due  to  SLM  is  no  more  than  the  reflectivity  of  the  beam  spTtter  plus  the 
maieridl  bulk  absorption  due  to  the  pholorefractive  crystals.  The  beam  splitter  can  be 
chosen  such  that  the  refloctivitv  is  sriiall  (e.g.,  5%  or  less)  so  that  such  energy  loss  is 
m iri 111 !  1 /ed.  The  hologram  formation  in  the  piiotorefraclive  medium  will  limit  the  recon- 


19 

C1017.STD/b|e 


Rockwell  International 

Science  Center 


SC5302.fr 

figuration  time.  Once  the  interconnection  pattern  is  formed  inside  the  photorefractive 

crystal  as  a  hologram,  such  a  scheme  is  capable  of  providing  the  interconnection  for  high 

data  rate  transmission.  In  such  an  interconnection,  the  output  of  each  laser  can  be  input 

to  any  one  or  all  of  the  detectors.  Optical  phase  conjugation  can  also  be  used  in 

conjunction  vnth  the  two-v.ave  mixing  to  correct  for  any  phase  aberration  that  may  be 

23 

caused  bv  the  crystal  imperfection. 

3.5.2  Experimental  l3emonstration 

An  experimental  configuration  for  a  1  x  32  interconnection  using  this 
holographic  approach  is  shown  in  Fig.  12.  A  Ronchi-Ruling  (5C  line-pairs  per  inch)  is  used 
as  a  32-eiemeni  binar\  vector  mask.  5%  of  the  collimated  input  beam  transmitted 
through  the  beam  splitter  BS  is  expanded  vertically  by  a  pair  of  cylindrical  lenses  (Lcl 
and  Lc2)  to  illuminate  the  vector  mask.  95%  of  the  input  energy  reflected  by  the  beam 
splitter  BS  serves  as  the  pump  beam  which  interacts  with  the  spatially  modulated  signal 
beam  insicio  the  barium  titanate  crystal.  The  crystal  and  the  mask  are  placed  at  the 
bdCk-fcxial  plane  and  the  front-focal  plane,  respectively,  of  a  Fourier  transform 
cvlmdrica!  lens  (Lc3).  The  crystal  is  oriented  so  that  energy  is  transferred  from  the 
pump  beam  to  the  signal  beam.  The  amplified  signal  output  carrying  the  interconnection 
pattern  is  rmmaged  by  the  cylindrical  lens(Lc4)  on  the  screen  SC. 

CYllfglrl'ipAI  H(  AM  I  It  '  " 


I  ml  '(  tf'  r  n. 


1  ig.  1  2  A  scheii  aiic  diagram  of  the  experimental  configuration  lor  a 
1  -  32  interconnection  using  a  photorefractive  barium  titanate 
crystal. 
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The  intensity  distribution  at  the  output  image  plane  is  shown  in  Fig.  13.  The 
experimental  results  is  shown  in  Fig.  14  for  the  case  when  only  three  out  of  the  32 
channels  are  "ON".  The  table  in  the  right  hand  side  of  the  figure  compares  the  energy 
efficiencies  achieved  with  and  without  the  photorefractive  holograms.  For  the  latter 
case,  the  beam  splitter  BS  is  removed  (see  Fig.  12)  so  that  all  the  input  energy  is  directed 
into  the  signal  channel.  Note  that  the  energy  efficiency  is  improved  by  more  than  a 
factor  of  7. 


.C.IJ. 

32  ELEMENT  VECTOR  MASK 


(bl  OUTPUT  THROUGH  BaTiOa  CRYSTAL  WITH  PUMP  OFF 
(cl  OUTPUT  THROUGH  BaTi03  CRYSTAL  WITH  PUMP  ON 
(I  a  .  AMPLIFIED  VIA  TWO  WAVE  MIXING) 


Fig.  13  Intensity  distribution  at  the  output  image  plane. 


3.5.3  Demortstration  of  High  Data  Rate  Transmission 

Ue  have  experimentally  demonstrated  that  the  holographic  interconnection 
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using  photorefractive  holograms  can  accept  very  high  data  rate  signals  [see 

Appendix  5.5J  even  though  the  reconfiguration  time  is  limited  by  the  photorefractive 

21 

grating  formation,  which  is  typically  of  the  order  of  milliseconds  at  modest  intensity. 

To  demonstrate  this  fact,  temporal  modulation  was  impressed  on  a  laser  beam  (argon, 

X  -  514.5  nm)  using  an  acoustooptic  device  to  simulate  a  signal,  which  is  to  be 
interconnected  Vvith  some  output.  The  signal  used  was  a  pulse  tram  of  frequency 
f^  -  0.S33  MHz  with  each  pulse  being  ~  0.2  us  wide.  This  rate  is  clearly  much  higher  than 
the  reciprocal  of  the  photorefractive  response  time.  The  modulated  laser  beam  was  then 
split  into  two  beams  and  mixed  in  the  crystal  as  described  before,  and  the  amplified 
probe  beam  was  monitored  with  a  photodetecior.  The  upper  oscilloscope  trace  in  Fig.  15 
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SC««  1*0 

32  ELEMENT 
VECTOR  MASK 


Fig.  14  Intensity  distribution  at  the  output  image  plane  v/hen  all  but 
three  of  the  signal  channels  are  "ON."  The  table  on  the  right- 
hand  side  compares  the  energy  efficiencies  achieved  using  the 
conventional  approach  and  the  holographic  approach. 
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Fig.  15  (a)  Signal  carried  by  the  probe  beam  and  (b)  amplified  signal  carried  by  the 

probe  beam  after  undergoing  two-wave  mixing. 


shows  the  input  probe  signal,  and  the  lower  one  shows  the  amplified  probe  signal.  The 
results  show  a  steady-state  response  in  which  the  temporally  modulated  pump  and  probe 
beams  interact  simply  by  diffracting  off  the  stationary  index  grating  that  is  created  in 
the  crystal  after  the  photorefractive  response  time.  As  this  experiment  was  performed 
merely  to  demonstrate  the  high  signal  bandwidth  of  the  system,  optimization  of  the 
parameters  was  not  done,  and  the  results  in  Fig.  15  can  clearly  be  improved. 


22 

C10178TD/bje 


Rockwell  International 

Science  Center 


SC5502.fr 

3.5.4  Energy  Efficiency  of  Photorefractive  Interccxinection;  TheoretiCciJ  Calculation 
and  Experimental  Measurement  f^e  Apperxjtx  5.5,  5.6,  and  5.7] 

Referring  to  Fig.  16,  we  can  estimate  the  energy  efficiency  of  photorefractive 
crossbar  switch  as  follows.  Let  1^  be  the  input  beam  intensity  and  R  be  the  reflectance 
of  the  beam  splitter.  It  is  legitimate  to  assume  that  the  beam  splitter  is  practically 
lossless.  Ue  may  also  assume  that  the  surface  of  the  photorefractive  crystal  is 
antireflection-coated  so  that  the  Fresnel  reflection  loss  can  be  neglected.  Under  il.ese 
conditions,  the  two  beams  that  arrive  at  the  photorefractive  crystal  have  energies 
1^(1  -  R)  and  I^Rt/N,  respectively,  where  t  is  the  intensity  transmittance  of  the  "on”  cells 
of  the  mask.  Inside  the  crystal,  these  two  beams  undergo  photorefractive  coupling.  As  a 
result,  most  of  the  energy  of  the  pur>ip  beam  1^(1  -  R)  is  transferred  to  the  probe  beam 
Ij^Rt/N,  which  contains  the  intec connection  pattern.  The  energy  efficiency  can  be  easily 
derived  and  is  given  by" 


out 


_ 1  +  m 

N  1  +  m  exp(-’rL) 


exp  (-aL) 


(ii) 


where  m  is  the  beam  intensity  ratio, 


m 


(1  '  R)N 

tR 


(9) 


and  L  is  the  interaction  lengfh,  ,  is  the  coupling  constant,  and  a  is  the  bulk  absorption 
coefficient,  f'or  photorefractive  crystals  such  as  BaTiO^  and  SBN,  the  coupling  constant 
IS  very  large  (i.e.,  yL  >>  1),  The  efficiency  can  be  written  approximately 


r: 


(1  -  R)l 


exp(-aL) 


(10) 
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Fig.  16  A  schematic  diagram  showing  the  intensities  of  various  beams  in  a  N  .  N 
permutational  interconnection. 

We  notice  that  for  large  N,  the  energy  efficiency  is  limited  by  the  crystal  bulk  absorption 
exp(-aL)  and  is  maximized  by  using  a  beam  spliter  with  a  very  small  reflectance  R  (i.e., 
R  ~  0). 

In  Fig.  17,  the  maximum  energy  efficiency  (with  optimum  beam  splitting  ratio) 
predicted  by  Eq.  (8)  (asstiming  that  the  insertion  loss  of  the  mask  is  negligible,  i.e.,  t  =  1 
and  that  the  absorption  coefficient  a  =  0)  is  plotted  against  the  dimension  of  the  array  N 
for  different  coupling  strength  yL.  The  diagonal  straight  line  represents  the  intrinsic 
fan-out  loss.  For  cases  where  a  f  0,  the  curve  remains  valid  provided  that  the  vertical 
scale  is  multiplied  by  exp(-aL). 

Using  an  experimental  configuration  illustrated  schematically  in  Fig.  18,  we 

have  measured  the  energy  efficiency  of  the  photorefractive  interconnection  holograms  in 

1 3 

a  barium  titanate  crystal  as  a  function  of  the  fan-out  loss.  Instead  of  expanding  the 
signal  beam  through  a  vector  mask  as  described  earlier  (see  Fig.  12),  a  neutral  density 
filter  (NDF)  in  the  signal  arm  is  used  to  simulate  the  fan-out  loss.  The  experimental 
results  are  shown  in  Fig.  19.  The  solid  line  is  a  theoretical  fit  using  Eq.  (8)  with  aL  and 
-fL  as  the  fitting  parameter.  The  results  (oL  =  9.2,  yL  =  1.1)  agree  with  those  obtained 
from  other  independent  measurements. 

3.5.5  Maximizing  Energy  Efficiency  by  Matched  Amplification  via  Fourier 
Transformation  (sec  Appendix  5.7) 

In  the  experiments  discussed  above,  both  signal  and  pump  beams  are  Gaussian 
beams  with  no  spatial  intensity  pattern.  In  the  interconnection  applications,  these  beams 
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IC4}>Tt 

ENERGY  EFFICIENCY  vs  ARRAY  DIMENSION 


Fig.  17  The  theoretical  energy  efficiency  as  a  function  of  the  dimension  of  the  array 
as  predicted  by  Eq.  (8)  for  t  =  I,  a  =  0,  and  an  optimum  R  --  0. 


SC44329 


fhg.  18  A  schematic  diagram  of  the  experimental  configuration  for  measurement  of 
the  energy  efficien>_y  of  photorefractive  interconnection  holograms  in  a 
barium  litanate  crystal. 
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N 

Fij,.  19  Energy  efficiency  (n)  of  photorefraclive  two-beam  coupling  in  a  barium 
titanaie  sample  as  a  function  of  the  transmittance  of  a  neutral  density  filter 
placed  in  the  signal  input  arm.  The  transmittance  is  labeled  1/N  to  relate  it  to 
the  fanout  loss  of  a  N  x  N  permutation  crossbar  network. 

are  spatially  modulated.  The  energy  efficiency  of  the  photorefractive  interconnection 
depends  on  the  spatial  overlap  of  the  beams.  We  have  developed  the  utilization  of 
f'ourier  transform  to  achieve  maximum  beam  overlap  and  thus  to  achieve  maximum 
energy  efficiency. 

It  is  important  to  note  that  although  the  beamlets  in  the  signal  beam  and  the 
pump  beam  are  Intrinsically  different  because  of  the  interconnection  pattern,  the 
individual  pixels  can  have  identical  shape  (e.g.,  square).  Complete  overlap  is  possible  in 
the  Fourier  domain,  provided  all  the  pixels  are  identical.  Let  s(x,y)  be  the  aperture  of  an 
individual  pixel,  and  o(u,v)  be  the  Fourier  transform.  Fourier  transformation  of  a  beam- 
let  of  the  pump  beam  and  the  corresponding  column  of  beamlets  of  the  signal  beam  can 
be  performed  by  a  lens.  The  resulting  amplitude  distributions  at  the  rear  focal  plane  are 
given  by 


Pump:  P(u,v)  =  o(u,v) 

(11) 

Signal;  S(u,v)  :  o(u,v)  exp  [  i(t  (u,v)  1 
n 

(12) 
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respeciivel> ,  where  u,v  are  coordinates  in  the  Fourier  plane.  The  summation  in  Eq.  (12) 
is  over  all  the  windows  open  in  the  column,  and  it)p|(u,v)  is  a  phase  which  depends  on  the 
position  of  the  window. 


Ue  note  that  the  Fourier  transform  of  a  colunm  of  the  signal  beam  consists  of 
a  linear  superposition  of  identical  patterns  each  with  a  different  phase  factor.  Such 
phase  factors  (expfiyj^)  are  due  to  the  relative  positions  of  the  windows  in  each  column. 
Although  these  phase  factors  may  lead  to  interference  structures,  all  the  energies  are 
confined  under  the  same  envelop  c(u,v).  A.s  a  result  of  the  shift  invariance  in  Feurie^ 
transformation,  each  of  the  signal  beamlets  overlaps  completely  with  the  pump  beam  at 
the  Fourier  plane.  Thus,  maximum  energy  coupling  is  achieved. 


The  experimental  configuration  is  illustrated  in  Fig.  20.  The  output  from  an 
argon  laser  (514.5  nm)  is  spatial-filtered  and  expanded  to  form  a  collimated  beam  with 
2  c:n  diameter.  A  variable  beam  splitter  consisting  of  a  half-wave  plate  and  a  polarizing 
beam  splitter  cube  is  used  to  split  the  incoming  beam  into  two,  a  pump  beam  and  a  signal 
beam.  Another  half-wave  plate  is  used  to  rotate  the  polarization  of  the  reflected  (signal) 
beam  by  90°  from  the  direction  perpendicular  to  the  plane  of  the  paper  to  the  in-plane 
direction.  A  polarizer  placed  downstream  further  filters  out  the  residual  perpendicular 
component.  A  10  x  10  binary  matrix  mask  is  used  to  encode  a  spatial  pattern  onto  the 
signal  beam. 

To  maximize  the  spatial  overlap  of  the  pump  and  the  signal  beams  inside  the 
crystal  via  the  properties  of  Fourier  transform,  a  pump  mask  with  a  single  aperture 
identical  to  the  unit  cell  of  the  signal  mask  is  used.  The  two  masks  are  placed  at  the 
front  focal  planes  of  two  Fourier  transform  lenses  FL^  and  FL2  (of  identical  focal  length 
f  =  25  cm),  respectively.  The  crystal  is  located  at  the  conunon  Fourier  plane  of  the  two 
masks  (Fig.  20).  The  shift-invariant  property  of  Fourier  transformation  ensures  that, 
apart  from  a  phase  factor,  the  diffraction  pattern  from  each  signal  channel  overlaps  with 
that  from  the  pump  inside  the  crystal.  The  intensity  patterns  of  the  pump  and  the  signal 
at  the  image  plane  and  the  Fourier  plane  are  shown  in  Fig.  21.  Note  that  the  two 
intensity  distributions  at  the  Fourier  plane  differ  significantly  in  fine  structures  due  to 
multiple  beam  interference  among  the  various  signal  channels. 


27 

CI0178TD/b)e 


Rockwell  International 

Science  Center 


SC5502.fr 


Fig.  20  An  experi menial  configuration  for  a  I-to-N  x  N  (for  N  =  10)  broadcasting 
network  using  photorefractive  holograms  at  the  Fourier  domain. 


Intensitity  Patterns 


Image  Plane 


Fourier  Plane 


Probe 


Pump 


F-ig.  21  The  intensity  patterns  of  the  masks  for  the  probe  and  the  pump  beams  at  the 
image  plane  and  the  Fourier  plane. 
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To  investigate  how  the  energy  efficiency  depends  on  the  number  of  signal 
channels,  we  used  a  variable  rectangular  aperture  in  front  of  the  signal  mask  to  vary  the 
number  of  channels  from  100  to  1.  In  each  case,  we  also  adjusted  the  laser  power  and  the 
variable  beam  splitter  so  that  both  the  pump  power  and  the  total  signal  input  power  are 
fixed  at  600  uW  and  6  u\^,  respectively.  Figure  22  is  a  plot  of  energy  efficiency  vs 
number  of  signal  channels.  The  energy  efficiency  turns  out  to  be  fairly  insensitive  to  the 
number  of  signal  channels.  The  crystal  orientation  used  in  this  particular  experiment  is 
far  from  optimum,  and  one  can  expect  significant  improvement  in  energy  efficiency  once 
optimized. 
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Fig.  22  Energy  efficiency  (n)  as  a  function  of  number  of  signal  channels  (N)  m  a  1-to- 
N  .  N  broadcasting  configuration  using  photorefractive  holograms  in  a  barium 
titanate  sample. 

3.5.6  Contrast  Enhancement  by  Double  Passage  (see  Apperxlix  5.9) 

All  of  our  experimental  results  reported  earlier  are  obtained  with  photographic 
masks  replacing  the  SLM.  To  demonstrate  the  reconfigurability  and  to  study  the 
reconfiguration  lime,  one  can  no  longer  avoid  the  use  of  SLM  which,  in  general,  have 
relatively  poor  contrast  ratio  (of  the  order  10  to  100).  In  the  parallel  matrix-vector 
multiplicdliun  approach,  the  maximum  number  of  fan-out  channels  of  the  optical  cross¬ 
bar  can  be  limited  by  the  finite  contrast  of  the  SL\(. 
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As  a  simple  example,  consider  a  M  x  M  permutation  matrix  mask  with  contrast 
ratio  C  =  T  |/Tq,  where  T  ^  and  Tq  are  the  intensity  transmittance  of  the  "ON"  state  and 
the  "OFF"  state,  respectively.  U'e  assume  that  the  signal  from  each  source  is  intensity 
modulated  between  zero  and  peak  value  P,  the  signal  level  received  by  each  detector 
through  the  "ON"  cell  is  PTp  whereas  the  noise  level  received  through  the  (M-l)  "OFF" 
cell  is  P(M-1)Tq.  The  signal-to-noise  ratio  (S/N)  is  therefore  given  by 
S/N  =  Tj^/(M-1)Tq  =  C/(M-1)  =  C/M  for  M  >  1  .  The  condition  S/N  >  1  is  satisfied 

provided  that  the  number  of  fan-out  channels  (M)  is  less  than  the  contrast  ratio  (C). 

A  simple  method  to  improve  the  contrast  is  by  double  passage  through  the 
same  SLM  via  retrof lection  or  phase  conjugation^.  In  principle,  if  the  contrast  ratio  for 
signal  passage  is  C,  the  contrast  ratio  for  double  pasage  will  become  C  .  The  experi¬ 
mental  configuration  used  to  verify  this  fact  is  shown  schemaiicaly  in  Fig.  23.  An 
ex  panded  and  collimated  laser  beam  from  an  argon  laser  (514.5  nm)  is  transmitted 
through  a  circular  aperture  (diameter  =  14  mm)  located  at  the  input  plane.  Half  of  the 
aperture  is  masked  by  a  NDF  to  form  a  binary  amplitude  mask  at  the  input  plane.  The 
transmitted  beam  is  collected  by  a  lens  (L)  and  directed  into  a  barium  titanate  crystal. 
The  crystal  is  oriented  so  that  the  incoming  beam  is  phase-conjugated  back  via  self- 
pumped  mode.  After  reversing  through  the  aperture,  the  phase-conjugated  beam  is 
sampled  by  a  beam  splitter  through  imaging  optics  which  reimages  the  binar>’  mask  on 
the  output  plane.  The  input  intensity  distribution  (i.e.,  single-passage  through  the  mask) 
and  the  output  intensity  distribution  (i.e.,  double-passage  through  the  mask)  are  measured 
by  scanning  a  detector  with  a  small  aperture  (diameter  =  1  mm)  across  the  beam 
diameter  at  the  input  and  the  output  planes,  respectively.  The  experimental  results  (for 
a  mask  with  optical  density  OD  =  0.5)  are  shown  in  Fig.  24(a)  and  (b).  Apart  from  some 
imperfection  at  the  edge  of  the  output  image,  the  output  contrast  ratio  is  approximately 
the  square  of  that  of  the  input  as  expected.  Instead  of  a  phase-conjugate  mirror,  one  can 
also  place  a  plane  mirror  at  the  input  plane  immediately  after  the  mask  to  retroreflect 
the  beam  back  through  the  binary  mask.  The  resulting  intensity  distribution  is  shown  in 
Fig.  24(c).  The  phase  conjugation  has  the  key  advantage  of  being  intrinsically  self- 
aligned.  The  grating  formation  time  for  self-pumped  phase  conjugation  using  a  typical 
barium  titanate  crystal  at  a  few  tens  of  milliwatt  of  optical  pow'er,  is,  however, 
relativelv  ilow  (a  fraction  of  a  second  to  several  seconds).  Note  that  the  output  obtained 
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3.5.7  Experimental  Demonstration  of  Reconfigurability 

Using  an  experimental  configuration  similar  to  the  one  illustrated  in  Fig.  20, 
with  the  photographic  mask  (transparency)  replaced  by  an  SUM,  we  have  demonstrated 
the  reconfigurability  of  the  interconnection  by  writing  new  interconnection  matrix 
pattern  on  the  LCTV.  The  pattern  is  generated  by  an  IBM  PC  and  sent  to  the  LCTV, 
'A  ith  optical  power  of  the  order  of  tens  of  milliwatts,  the  frame  rate  is  currently  limited 
by  tr.e  photorefractive  response  time  to  a  few  frames  per  second.  As  explained  in 
Section  3.5.3,  the  data  rate  is  not  limited  by  the  photorefractive  response  time  and  can 
be  higher  than  several  megahertz.  To  enhance  the  contrast  of  the  binary  pattern,  the 
signal  beam  is  passed  through  the  SUM  twice  via  a  retroref lecting  mirror  (as  described  in 
Section  3.5.6)  prior  to  interacting  with  the  pump  beam  m  a  photorefractive  barium 
titanate  cr^  stal. 
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APPENDIX  5.1 


Optical  Matrix-Vector  Multiplication  Through  Four-Wave  Mixing 
in  Photorefractive  Media 
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Optical  matrix-vector  multiplication  through  four-wave  mixing 
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We  propose  and  describe  a  new  method  of  optical  matrix-vector  multiplication  by  usin^t  four-wave  mixinp  in 
photorelractive  media  I  siiip  a  BaTiO  crystal,  we  have  demonstrated  such  a  parallel  mullipiicatior.  Ttie  re>ult> 
are  presented  and  discussed. 


A  large  number  of  signal-  and  image-processing  algo¬ 
rithms  can  be  expressed  in  terms  of  matrix  operations. 
The  multiplication  of  two  matrices  is  one  of  the  most 
basic  operations  in  matrix  algebra.  Optics,  with  its 
inherent  parallelism,  can  offer  great  improvement  in 
the  speed  of  these  operations.  Optical  processors  for 
multiplying  two  matrices  have  been  described  in  the 
literature.'  Considerable  work  has  also  been  report¬ 
ed  on  performing  optical  matrix-vector  multiplica¬ 
tion.  ’  In  this  Letter,  we  describe  a  new  method  of 
performing  matrix-vector  multiplication  by  using  op¬ 
tical  four-wave  mixing  in  nonlinear  media.  In  addi¬ 
tion,  we  also  report  the  first  experimental  demonstra¬ 
tion  to  our  knowledge  of  such  an  optical  matrix  proces¬ 
sor  by  using  a  photorefractive  BaTiO,  crystal. 

Optical  four-wave  mixing  has  been  a  subject  of  con¬ 
siderable  interest  during  the  past  several  years.  Much 
attention  has  been  fticused  on  wave-front-correction*'  ' 
and  phase-conjugate  interferometry.^ In  the  area  of 
signal  processing,  optical  four-wave  mixing  has  been 
used  to  perform  spatial  information  processing.'-  im¬ 
age  subtraction.'-  '  and  logic  operations."  However, 
little  attention  has  been  paid  to  the  use  of  the  inher¬ 
ently  parallel  multiplication  nature  of  four-wave  mix¬ 
ing  for  optical  matrix  multiplication.  In  this  Letter 
we  describe  a  new  method  of  performing  matrix-vec¬ 
tor  multiplication  using  four-wave  mixing  in  nonlinear 
media. 

In  what  follows,  we  briefly  describe  the  basic  princi¬ 
ples  of  real-time  matrix  multiplication  using  optical 
four-wave  mixing  in  nonlinear  media.  For  the  sake  of 
convenience,  we  will  limit  ourselves  to  the  case  of 
square  matrices,  although  extension  to  nonsquare  ma¬ 
trices  will  be  straightforward. 

Referring  to  Fig.  1.  we  consider  a  scheme  that  is 
suitable  for  a  matrix-vector  multiplication.  Here,  as 
an  example,  let  us  consider  a  discrete  case  in  which  we 
need  to  carry  out  the  multiplication  of  an  .V-elemeni 
vector  and  an  .V  X  .V  matrix.  In  this  scheme,  the 
vector  i.t,  fanned  i/Ut  into  .V  rows  of  identical  vectors. 
These  .V  x  ,V  beamlets  are  directed  to  a  nonlinear 
medium.  The  matrix,  which  al.-'O  contains  ,V  X  ,V 
beamlets.  is  also  directed  to  the  medium  in  such  a  way 
that  each  heamlet  of  thi  matrix  is  counterpropagating 
in  a  direction  relative  to  the  corresponditif,  heamlet  of 


the  vector.  Thus,  in  the  medium,  there  are  .V  X  A' 
spatially  separated  regions,  each  of  which  is  pumped 
by  a  pair  of  counterpropagating  beamlets.  Now  a 
probing  beam,  which  consists  of  .V  X  A'  beamlets,  is 
directed  into  the  medium  in  such  a  way  that  each 
probe  heamlet  will  propagate  through  the  correspond¬ 
ing  intersection  region,  as  shown  in  Fig.  1.  The  probe 
beam  will  be  plane-wave  beamlets  propagating  in  par¬ 
allel.  All  the  beamlets  in  this  probe  beam  are  of  equal 
intensity.  As  a  result  uf  the  four-wave  mixing,  each 
probe  heamlet  will  generate  a  phase-conjugated  beam- 
let,  which,  within  a  proportional  factor,  can  be  written 
as  A/(i,  ilafjl-  By  using  a  cylindrical  lens,  a  summa¬ 
tion  over  j  can  be  obtained.  Thus  we  have 

bii)  =  X  aij),  (1) 

where  aij)  is  the  ;th  element  of  the  vector  a  and 
is  the  matrix  element.  Such  a  scheme  for  matrix- 
vector  multiplication  can  also  be  used  for  matrix-ma¬ 
trix  multiplication  by  decomposing  a  matrix  into  col¬ 
umn  vectors  and  then  multiplying  the  matrix  with 
each  of  the  column  vectors. 

The  probe  beam  can  also  simply  be  a  uniform  plane 
wave,  without  beamlets.  The  phase-conjugated  h^am 
will  consist  of  A'  X  A'  beamlets,  because  only  a  matrix 
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Eig.  2.  Experimenlal  setup  used  lu  demonstrate  llie  optical 
matrix -vector  multiplication. 


of  N  X  N  regions  is  pumped  by  counterpropagating 
beamlets.  The  matrix  and  the  vector  information  can 
actually  be  carried  on  any  two  of  the  three  incident 
beams  in  the  four-wave  mixing  process. 

The  summation  process  can  also  be  obtained  with¬ 
out  the  external  cylindrical  lens  by  using  a  different 
scheme,  described  as  follows.  Consider  a  probe  beam 
that  consists  of  a  column  of  N  equal  beamlets.  The 
probe  is  incident  into  the  nonlinear  medium  in  such  a 
way  that  each  of  the  probe  beamlets  is  made  to  propa¬ 
gate  through  a  row  of  A’  intersection  regions.  The 
phase-conjugation  process  automatically  performs  the 
summation  as  well  as  the  multiplication.  The  phase- 
conjugated  beam  is  thus  the  product  of  the  matrix- 
vector  multiplication. 

As  a  result  of  the  nature  of  the  four-wave  mixing 
process  in  the  medium,  this  new  matrix-vector  multi¬ 
plier  operates  on  the  field  amplitudes  and  thus  can  be 
used  to  handle  matrices  and  vectors  with  complex 
elements.  It  is  a  coherent  device  rather  than  an  inco¬ 
herent  one.  This  aspect  is  distinctly  different  from 
most  of  the  earlier  approaches,  which  are  all  incoher¬ 
ent.  When  the  device  is  operated  in  the  coherent 
mode,  the  phase  of  each  beamlet  must  be  maintained 
uniformly  over  the  transverse  dimension  of  the  beam- 
let.  In  addition,  such  phases  must  also  be  maintained 
fixed  in  the  summation  process.  If  these  phases  are 
not  uniform  over  the  beamlets,  the  final  step  becomes 
an  incoherent  summation  as  a  result  of  the  spatial 
averaging.  Under  such  circumstances,  this  matrix- 
vector  multiplier  operates  on  the  intensities  and  thus 
handles  only  positive  numbers. 

The  matrix-vector  multiplication  is  demonstrated 
experimentally  by  using  four-wave  mixing  in  a  photo- 
refractive  BaTiO.'  crystal.’  The  crystal  is  cut  in  such 
a  way  that  the  faces  are  all  perpendicular  to  the  princi¬ 
pal  axes.  All  the  beams  are  polarized  extraordinary  in 
the  xz  plane  and  are  incident  onto  the  a  faces  of  the 
crystal. 

In  our  experiments,  an  Ar-ion  laser  beam  at  wave¬ 
length  514. .5  nm  with  an  output  power  of  a  few  hun¬ 
dred  milliwatts  is  used  as  the  light  source.  The  laser 
beam  is  expanded  and  collimated  into  a  beam  size  of 
approximately  1  cm.  Figure  2  shows  our  experimental 
setup.  The  expanded  beam  is  then  split  into  three 
beams  by  using  beam  splitters  BSi  and  BS_.  To  dem¬ 


onstrate  the  principle  of  operation,  we  chose  the  fol¬ 
lowing  matrix  and  vector: 
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The  matrix  and  the  vector  information  is  imprinted 
onto  two  of  the  beams  by  using  transparencies  (MAi 
and  MA;  in  Fig.  2)  that  consist  of  circular  dots.  White 
circular  dots  represent  the  I’s,  and  the  dark  regions 
the  O’s.  Instead  of  using  a  cylindrical  lens  to  fan  out 
the  vector  into  five  identical  rows  as  described  earlier, 
we  simply  use  a  transparency  of  a  matrix  that  consists 
of  five  identical  rows  of  (1  0  1  1  1).  The  incident 
matrices,  corresponding  to  the  matrix  M  and  the  vec¬ 
tor  a,  are  shown  in  Figs.  3(a)  and  3(b).  Four  lenses,  Lj, 
L_.,  L3,  and  (with  focal  lengths  /j  =  /s  =  0.5  m  and  f  - 
=  (4  =  0.2  m)  are  used  to  image  the  transparencies  MAi 
and  MA.'  onto  the  center  of  the  nonlinear  crystal, 
where  the  diameter  of  each  beamlet  is  O.lo  mm  with  a 
center-to-center  distance  of  0.51  mm.  The  size  of  the 
whole  image  inside  the  crystal  is  2.2  mm  X  2.2  mm. 
The  uniform  beam  is  also  reduced  (and  collimated)  to 
a  diameter  of  3.6  mm  inside  the  crystal  by  lenses  L5 
and  Lr,  (with  focal  lengths  /,-  =  0.6  m  and  =  0.25  m, 
respectively).  The  angle  $  between  the  uniform  beam 
and  one  of  the  pump  beams  is  11°.  The  two  beams 
that  carry  the  matrix  and  the  vector  information,  re¬ 
spectively,  are  incident  onto  the  opposite  a  faces  of  a 
photorefractive  BaTiO:\  crystal.  The  incident  angle  d) 
is  approximately  20°  with  respect  to  the  surface  nor¬ 
mal.  As  a  result  of  the  four-wave  mixing,  a  phase- 
conjugated  beam,  consisting  of  a  5  X  5  beamlet  pat¬ 
tern,  is  generated  and  is  shown  in  Fig.  3(c).  A  cylin¬ 
drical  lens  is  then  used  to  perform  the  summation. 
The  resulting  product  vector  is  shown  in  Fig.  3(d).  By 
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Fip.  3.  Photographs  showing  the  matrix  and  the  vectors: 
(ai  input  matrix;  (b)  input  vector,  expanded  into  a  matrix 
with  five  identical  rows:  (c)  resultant  phase-conjugate  reflec¬ 
tion  iii<  trix;  (d )  resultant  vector  of  a  matrix -vector  multipli- 
catiiin  after  summation  of  tc)  by  a  cylindrical  lens. 
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measuring  the  intensity  of  the  individual  dots,  we 
found  that  the  product  vector  is  (3  2  0  1  4),  in  good 
agreement  with  our  prediction;  in  other  words. 
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The  intensities  of  the  product  vector  elements  are 
measured  by  integrating  over  the  cross  section  of  each 
beamlet.  The  experimental  error  is  less  than  5%. 
This  small  error  is  due  to  environmental  perturbations 
and  the  slow  response  of  the  photorefractive  crystal. 
Although  our  experiment  demonstrates  the  operation 
for  matrices  and  vectors  with  binary  numbers  (1,  0), 
this  matrix-vector  multiplier  can  also  be  operated  in 
the  analog  mode  to  handle  elements  with  gray  levels. 
Such  an  analog  operation  in  photorefractive  crystals 
requires  that  the  probe-beam  intensity  be  weak  so  that 
the  error  due  to  the  intensity  denominator  in  the  index 
modulation  can  be  neglected.  This  material-related 
error  disappears  when  the  matrix  and  the  vector  con¬ 
tain  only  binary  numbers  (1,  0). 

In  conclusion,  we  have  proposed  and  demonstrated 
a  new  method  of  matrix  operation  by  using  optical 
four-wave  mixing  in  nonlinear  media.  In  particular, 
we  have  demonstrated  experimentally  matrix-vector 
multiplication  in  a  photorefractive  BaTiO.)  crystal. 

The  authors  acknowledge  helpful  discussions  with 
M.  D,  Ewbank  and  M.  Khoshnevisan  (Rockwell  Inter¬ 
national)  and  J.  Neff  (Defense  Advanced  Research 
Projects  Agency). 
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ABSTRACT 

The  concept  of  using  color-multiplexed  four-wave  mixing  to  achieve  optical  matrix-matrix 
multiplication  with  parallelism  is  described.  Experimental  results,  for  the  case  of 
2*2  matrices,  using  two-color  four-wave  mixing  in  a  photore f ract i ve  crystal  are  presented. 

1  -  INTRODUCTION 

In  this  paper  we  propose  and  demonstrate  the  concept  of  color-multiplexed  four-wave 
itixing  to  achieve  optical  ma tr ix-ma t r i x  multiplication  (MMM)  with  parallelism.  This  is  an 
extension  of  our  earlier  work  on  optical  matrix-vector  multiplication  (MVH)^  using  rhotore- 
fractive  four-wave  nixing.  The  additional  degree  of  parallelism  required  by  M.MM  (N^  for  MMM 
versus  for  MVH)  is  provided  by  color-multplexing  the  nonlinear  wave-mixing  process.  In 
Section  2  the  background  for  optical  matrix  processors  is  briefly  reviewed,  and  the  general 
idea  of  using  color-multiplexing  to  convert  a  MVM  scheme  into  MMM  is  described.  In  Sec¬ 
tion  3  we  present  a  specific  approach  to  apply  the  color  multiplexing  technique  to  photore- 
fractive  four-wave  mixing  for  MMM  with  full  parallelism.  The  experimental  details  and  re¬ 
sults  are  discussed  in  Section  4.  In  Section  5  we  summarize  our  results  and  discuss  some  of 
the  practical  difficulties  and  limitations,  as  well  as  some  advantages  of  this  approach. 

2.  BACKGROUND 

Matrix  algebra  has  been  used  extensively  as  a  convenient  mathematical  language  to  formu¬ 
late  the  problems  frequently  encountered  in  signal  and  image  processing  and  numeric 
computing.  Sooe  of  the  more  recent  applications  in  optical  computing  include  optical 
i nterconnections , 2  optical  neural  network  processors,^  and  optical  artificial  intelli¬ 
gence.'’'^  In  particular,  two-dimensional  format  of  a  matrix  can  be  conveniently  represented 
as  a  spatial  intensity  distribution  of  an  optical  beam  over  a  plane.  In  this  context,  a 
matrix  can  be  viewed  as  a  two-dimensional  image  and  can  be  transferred  to  any  designated 
plane  and  scaled  to  appropriate  size  by  conventional  optical  technique.  The  two 
computational  primitives  involved  in  matrix  algebra,  namely  multiplication  and  addition,  can 
also  be  carried  out  optically  using  simple  reflective,  refractive,  or  diffractive  optics  or 
their  combinations.  Various  aspects  of  optical  matrix  processors  have  been  discussed  fairly 
extensively  in  some  recent  review  articles. In  the  following,  we  describe  briefly  only 
those  directly  relevant  to  this  work. 

A  basic  architecture  for  an  incoherent  optical  MVM  scheme  for  performing  discrete 
Fourier  transforms  was  reported  by  Goooman  et  al.  A  schematic  diagram  of  this  architec¬ 
ture  fcr  multiplying  a  4  *  4  matrix  and  a  4  »  1  vector  is  illustrated  in  Fig.  1.  In 
general,  a  linear  array  of  "N"  sources  (light  emitting  diodes  or  semiconductor  lasers,  for 
instance)  is  used  to  represent  the  input  vector  (dimension  =  N);  the  value  of  each  element 
of  the  input  vector  is  encoded  as  the  optical  power  of  each  of  the  source.  The  source  array 
can  be  addressed  electronically  in  parallel.  A  transparency  or  a  spatial  light  iDOCulator 
(SLM)  is  usee  to  encode  the  input  matrix  with  the  value  of  each  matrix  element  represented 
by  the  (local)  optical  intensity  transmittance  at  the  corresponding  location  on  the  matrix 
mask.  The  apwiopriate  pixel-by-pixel  multiplication  as  prescribed  by  the  definition  of  MVM 
is  achieved  by  projecting  the  source  array  to  match  the  matrix  mask,  using  anamorphic  optics 
(omitted  in  Fig.  1  for  the  sake  of  clarity).  The  required  summations  are  then  carried  out 
by  another  set  of  anamorphic  optics  (also  omitted  in  Fig.l)  that  projects  the  transaitted 
light  onto  a  linear  detector  array.  The  output  of  the  detector  represents  the  product. 

When  the  sources  are  addressed  in  parallel,  the  two  dimensional  matrix  mask  allows  one  to 
carry  out  the  MVM  operation  with  full  parallelism  (of  order  N^). 

The  computational  complexity  for  the  multiplication  of  two  N  x  H  matrices  is  of  the 
order  N^.  In  other  words,  it  requires  multiplications  and  summations  to  multiply  two 
IJ  X  t;  matrices.  The  problem  can  also  be  decomposed  into  N  sub-problems  of  MVM.  This  is 
illustrated  in  Fig.  2  for  the  special  case  of  N  =  3 .  Since  the  columns  with  all  tf,e  ele¬ 
ments  equal  to  zero  can  be  ignored,  each  sub-problem  essentially  involves  only  one  MVM. 

Using  the  configuration  shown  in  Fig.  1,  one  can,  therefore,  complete  the  multiplication  of 
two  IJ  *  N  matrices  Hi  and  M2  in  N  clock  cycles  by  clocking  into  the  source  array  a  sequence 
^f  H  vectors  (i  =  1  to  N)  which  together  constitute  the  matrix  Ml;  the  other  matrix  M2  is 
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MATRIX  MASK 


Ficure  1.  The  incoherent  matr ix-vector  multiplication  scheme; 

the  cylindrical  optics  are  omitted  for  clarity. 
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Figure  2.  Decomposition  of  a  problem  of  matr i x-ca t r ix  multiplication 
into  the  sub-problems  of  matrix-vector  multiplication. 

represented  by  the  two-dimensional  mask  as  described  earlier.  In  this  approach,  ^  storage 
device  is  required  at  the  output  end  since  the  final  result,  which  takes  the  form  of  N 
column  vectors,  is  obtained  sequentially  in  H  steps.  In  addition  to  the  t ime-oul t iplexi ng 
approach  with  parallelism  described  above,  one  can  also  use  a  color  (wavelength)  multi¬ 
plexing  scheme  to  achieve  MMM  with  full  (N^)  parallelism. 

A  schematic  diagram  illustrating  the  basic  principle  of  the  color  multiplexing  MMM 
scheme  is  shown  in  Fig.  3.  Each  constituent  column  vector  Vj^  of  the  matrix  M|  is  color 
encoded  by  illuminating  it  with  a  unique  color  at  a  unique  input  angle.  Chromatic  (as  well 
as  angular)  multiplexing  is  achieved  by  combining  the  output  color  beam  strips  into  a  single 
column  with  the  aid  of  color  dispersion  in  a  prism.  Multiplication  in  parallel  of  each 
color  encoded  vector  with  matrix  Mt  is  accomplished  as  described  earlier  (see  Fig.  1). 

The  output,  after  chromatic  (and  angular)  demultiplexing  via  another  prism,  is  resolved  into 
the  product  matrix  which  can  be  image  onto  a  two-dimensional  array  of  detectors.  In 
practice,  one  nay  have  to  calibrate  the  output  to  correct  for  the* spectral  non-uniformities 
of  the  souc'e,  the  matrix  mask  and  the  detectors. 

The  color  cultiplexing  technique  is  by  no  mean  limited  to  the  incoherent  matrix  process¬ 
ing  scheme  described  above.  The  incoherent  approach  is  chosen  as  an  example  in  this  section 
to  explain  the  basic  idea  of  color  multiplexed  matrix  processors  mainly  because  of  its 
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Figure  3.  A  color-mul t i olex i ng  scheme  for  incoherent 
optical  matrix-matrix  multiplication. 

simplicity.  It  also  serves  as  a  reference  for  comparison  with  the  coherent  processor  using 
pho tore f rac t i ve  four-wave  mixing  to  be  described  in  the  next  section. 

3.  OPTICAL  MATRIX  PROCESSORS  DASED  0»  NONLINEAR  FOUR-WAVE  MIXING 

The  investigations  of  nonlinear  optical  four-wave  mixing  for  optical  phase  conjugation 
and  spatial  convolution/correlation  were  repor ted  * ^ *  I  2  several  years  ago.  Recently,  the 
applications  of  nonlinear  four-wave  mixing  for  other  optical  processings  such  as  MVM, ^ 
parallel  logic  operations ,  ^  ^  optical  pattern  recognization  for  symbolic  substitution,  l**  and 
residue-arithmetic  computation^^  have  also  been  demonstrated.  Most  of  these  applications 
utilize  the  same  basic  multiplicative  properties  of  four  wave  mixing.  They  differ  only  in 
the  input  spatial  patterns  and  in  their  particular  uses  of  lenses,  some  for  imaging  and 
others  for  spatial  Fourier  transform.  A  configuration  for  MVM  is  illustrated  schematically 
in  Fig.  4.  A  catrix-mask  M  is  illuminated  with  a  uniform  coherent  beam  and  imaged  onto  the 
nonlinear  medium  from  one  side.  The  mask  representing  the  vector  (ai,  a2,  a  3)  is 
illuminated  from  the  other  side.  The  transmitted  beams  are  expanded  by  cylindrical  optics, 
into  three  identical  row  vectors  which  align  pixel-by-pixel  with  the  matrix  M.  In  practice 


•3 


Figure  4.  A  schematic  diagram  illustrating  the  basic  concept 
of  optical  matrix-vector  multiplication  using 
nonlinear  four-wave  mixing. 
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one  can  avoid  tj'.e  cylindrical  optics  by  using  a  two-dimension  mask  to  encode  directly  tliree 
(or  N,  in  general)  identical  rows  cf  the  input  vector.  This  allows  one  to  use  identical 
imaging  optics  in  both  input  ports  and  thus  achieve  better  beam-overlap.  Inside  the 
nonlinear  crystal,  there  are  nine  pairs  of  counter-propagating  beamlets,  representing  the 
elements  of  the  matrix  M  and  the  vector  a.  When  a  uniform  probe  beam  is  injected  into  the 
nonlinear  crystal,  optical  phase  conjugation  takes  place  via  the  four-wave  mixing  process. 
Under  appropriate  conditions,  the  phase-conjugate  output  intensity  is  proportional  to  the 
products  of  the  three  ir.out  intensities.^*  A  beam  splitter  is  used  to  sample  the  phase- 
conjugate  output,  which  consists  of  nine  beamlets  representing  the  pixel-by-pixel  product  of 
the  matrix  and  the  vector.  A  cylindrical  lens  is  then  used  to  perform  the  row-wise  addi¬ 
tions  of  the  beamlets  to  obtain  the  product  vector  (bi,  b2,  bj)  at  the  output  plane.  The 
experimental  detail  and  results  demonstrating  the  multiplication  of  a  binary  matrix  (5  *  5) 
and  a  binary  vector  using  a  photoref ract ive  barium  titanate  crystal  as  the  nonlinear  medium 
were  reported  recently.^ 

Experimental  observation  of  simultaneous  optical  phase  conjugation  of  several  colors  of 
input  from  a  si.ngle  barium  titanate  crystal  has  been  reported  by  Chang  et  al.l’  Such  multi¬ 
color  phase  conjugation  via  four-wave  mixing  in  photore f r ac t i ve  crystal  can  be  employed  to 
perform  MMH.  .“igure  5  shows  a  schematic  diagram  of  such  a  matrix  processor.  This  configu¬ 
ration  is  basically  the  same  as  the  one  shown  in  Fig.  4  except  for  the  addition  of  two  dis¬ 
persive  elements,  one  in  the  input  for  color-multiplexing,  and  another  in  the  output  for 
color-demultiplexing.  Notice  that  the  input  matrix  A  is  decomposed  into  its  constituent 
vectors  encoded  with  different  colors  as  described  in  the  previous  section.  The  beam  illu¬ 
minating  the  other  matrix  M  and  the  probe  beam  consists  of  all  the  colors  uniformly  distri¬ 
buted  over  the  beam  profiles.  As  a  result  of  fogr-wave  mixing  of  each  of  the  spectral 
component  independently,  pixel-by— pixel  multiplication  with  parallelism  is  achieved.  The 
signal  to  noise  ratio  of  the  output  depends  critically  on  the  the  pixel  inteaction,  and 
hence  the  beam  alignment  of  each  of  the  spectral  components  inside  the  crystal.  In  the 
scheme  that  we  propose  (see  Fig.  5),  the  alignment  is  accomplished  in  two  steps.  The  input 
angles  for  each'  of'  the  spectral  components  and  the  associated  vector  masks  are  first  adjus¬ 
ted  so  that  the  output  appears  as  a  single  column  consisting  of  all  the  spectral  components. 
Cylindrical  optics  are  then  used  to  expand  this  single  column  to  coincide  with  the  matrix 
inside  the  nonlinear  medium.  This  procedure  ensures  that  all  the  spectral  components  are  in 
unison  and  car,  be  aligned  simultaneously  in  the  second  step.  Cross-talk  due  to  interactions 
among  different  colors  does  not  exist,  provided  the  frequency  separation  between  the  colors 
are  far  beyond  the  tempxbral  (frequency)  bandwidth  of  the  nonlinear  process.  This  is  usually 
the  case  for  most  of  the  photoref ractive  materials  since  their  response  time  typically 
ranges  from  micro-second  to  second. 


Figure  5.  A  schematic  diagram  illustrating  the  basic  idea  of 
optical  matrix-matrix  multiplication  using  color- 
multiplexed  four-wave  mixing. 


SPiE  Vol  8S1  Opiicat  CofT^puimg  and  NonhneBr  Materials  {1 988}  /  253 


4.  EXPERIMENTAL  CONFIGURATION  AND  RESULTS 


Using  a  photore f r act ive  barium  titanate  crystal  as  the  nonlinear  medium  and  an  argon  ion 
laser,  we  have  demonstrated  the  multiplication  of  2  *  2  binary  matrices.  The  experimental 
configuration  is  shown  schematically  in  Fig.  6.  Broad  band  reflectors  are  used  for  the 
laser  resonator  so  that  it  oscillates  at  more  than  one  wavelength.  Specifically,  we  operate 
the  laser  (Spectra  Physics  Model  171)  at  five  wavelengths  (514.5  nm,  496.5  nm,  488  nm,  476.5 
nin,  and  457.9  nm)  in  the  blue-green  region  of  the  visible  spectrum.  The  5-color  output  beam 
is  expanded  to  about  1  cm  in  diameter  and  collimated  by  the  beam  expander/col 1 ima tor  (BE). 

A  neutral  density  filter  (NDF)  with  OD  =  1  is  used  to  reduce  the  phase  conjugate  feedbac)t 
into  the  laser,  and  a  polarizing  beam  splitter  cube  (PBS)  is  used  to  filter  out  the  s- 
polarization  component.  The  beam  splitter  BS i  (R  =  50%),  lens  Li  (f  =  75  cm),  mirror  M2  and 
lens  L2  (f  =  33.3  cm)  constitute  the  optics  for  the  probe  beam.  A  2  x  2  matrix  mask  MM  1  is 
illuminated  with  one  of  the  pumping  beams  split  off  from  the  beam  splitter  BS2.  The  folding 
mirror  M2,  and  the  two  lens  imaging  system  (lenses  L3  (f  =  50  cm)  and  Lti  (f  =  30  cm)]  are 
used  to  project  a  reduced  image  of  the  mask  MMi  onto  the  crystal.  The  angle  6  subtended  by 
the  two  beams  is  about  7®.  The  beam  transmitted  through  BS 2  is  sent  through  a  prism  Pi  to 
angularly  separate  all  the  5-color  components  two  of  which  (  476.5  nm  and  514.5  rim)  are 
selected  to  illuminate  the  masks  MM2  and  MM3,  respectively.  The  two  color-encoded  patterns 
are  combined  at  bean  splitter  ES 3  and  imaged  onto  the  crystal  through  the  folding  mirror  M?, 
and  the  two  lenses  Ls  (f  =  50  cm)  and  L6  (f  =  30  cm).  For  proper  operation,  beam  alignm.ent 
and  pixel  interaction  inside  the  crystal  are  critical.  The  barium  titanate  crystal  (5  mim  x 
4  mm  X  4  mm)  is  cut  so  that  the  angle  «  between  the  c-axis  and  one  of  the  input  surface  is 
about  30®.  This  is  approx im, ate ly  equal  to  the  angle  between  the  c-axis  and  the  crating 
(normal)  vector  when  the  incidence  of  the  two  writing  beams  are  near  normal  (see  Fig.  6). 
Such  a  configuration  allows  us  to  maximize  the  effective  electro-optic  coefficient  without 
using  oblique  incidence,  which  reduces  the  useful  area  of  the  nonlinear  medium.  The  phase 
conjugate  beam  encoded  with  the  pixel-by-pixel  multiplication  of  each  of  the  colors  is 
sampled  by  the  beam  splitter  BS^,.  In  the  output  port,  the  folding  mirror  Ms  and  lens  L?  (f 
=  80  cmi)  im.age  the  center  of  the  crystal  onto  the  screen  S,  the  prism  P2  angularly  separates 
the  color  components  (color  cem.ul  t  iplex  ing )  of  the  output  beam,  and  the  cylindrical  lens  CL 
carries  out  the  row-wise  sumrr.ation  for  each  of  the  colors. 


Figure  6.  The  layout  for  demionstrat ing  the  multiplication  of  two  2x2 

m.atrices  using  color-multiplexed  photor  e  f  r  ac  1 1  vo  four-wave  mixing. 

Figure  7a  shows  the  two  m.atrices  M  and  A,  the  decomposition  of  m^atrix  A  into  column 
vectors,  and  the  color  encoding  scheme.  The  three  masks  M.M  1 ,  MM3  and  MM2  (see  rig.  6)  for 
encoding  tne  matrix  M  and  the  first  and  second  column  of  the  m.atrix  A,  respectively,  ate 
shown  in  Fig.  7b.  The  experimental  results  are  shown  in  Fig.  8.  Figure  8a  is  a  picture  of 
the  output  pattern  on  the  screen  S  when  the  summing  cylindrical  lens  (CL)  is  removed;  the 
tcl^afive  optical  power  of  each  pixel  monitored  by  a  detector  is  indicated  by  the  traces  from 
the  chart  recorder.  Figure  8b  is  the  final  result  representing  the  product  matrix  when  the 
cylindrical  lens  (CL)  is  in  place.  In  practice,  calibration  to  compensate  for  the  spectral 
characteristic  of  the  source,  masks,  and  detector  may  bo  necessary. 
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Figure  7.  (a)  The  two  matrices  (M  anc  A)  chosen  for  our  demonstration,  the  de¬ 

composition  of  A  into  column  vectors,  and  the  encoded  colors. 

(b)  The  three  matrix  masks  K.Mi,  KM3,  cr<3  M.M (see  Fig.  6)  for  encoding 
'  the  matrix  M  and  the  two  column  vectors  Ai  and  A2  derived  from  matrix  A. 
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Figure  8.  (a)  the  output  pattern  with  the  suttaing  cylindrical  lens  (CL)  (see  Fig.  6) 

removed;  the  relative  power  of  each  pixel  is  indicated  by  the  chart  re¬ 
corder  traces  and  the  color  of  the  pixel  is  indicated  below.  (b)  The  cor¬ 
responding  results  representing  the  product  of  M  and  A  when  the  suctming 
cylindrical  lens  (CL)  is  in  place. 
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SUMMARY  A.ND  CONCLUSION 


We  have  presented  the  basic  concept  and  experimental  results  for  optical  inatr  i  x-iaatrix 
multiplication  using  color-multiplexed  four-wave  mixing  in  a  nonlinear  medium.  Tc  Clarify 
the  concept  of  color-multiplexing,  an  incoherent  optical  approach  was  also  introduc<.d. 
Although  the  latter  is  relatively  simple  and  straight  forward  to  implement,  no  exf>er iroental 
results,  to  the  best  of  our  knowledge,  have  yet  been  reported.  For  multicolor  four-wave 
mixing,  a  photore f ract  i  ve  barium  titanate  crystal  is  chosen  for  our  experimental  demonstra¬ 
tion  because  of  its  strong  nonlinear  response  in  the  blue-green  region  of  the  visible  spec¬ 
trum  and  its  availability.  When  the  number  of  colors  is  limited  to  a  few,  no  significant 
reduction  in  phase-conjugate  reflectivity  of  each  color  component  is  observed.  However,  as 
the  number  of  colors  increases,  say  to  10  or  20,  it  is  anticipated  that  the  grating  con¬ 
trast,  and  hence  the  phase  conjugate  reflectivity  for  each  culor  will  be  reduced  because  of 
the  photo-excitation  of  charges  by  all  the  other  colors  sharing  the  same  volume.  Multiple 
quantum  well  structures  with  different  sets  of  layers  designed  to  enhance  the  optical  non¬ 
linearity  at  different  wavelengths  are  potential  candidates  for  color  multiplexed  four-wave 
mixing.  With  the  current  technology,  the  nonlinear  optical  approach  is  limited  by  the 
availability  of  an  optical  medium  with  sufficient  nonlinearity  (for  four -wave  mixing)  over 
a  wide  spectral  range  and  a  multi-color  coherent  source  that  matches  the  same  spectral 
range . 

In  our  exper  in.ental  demonstration,  we  chose  to  use  binary  matrices  and  monitored  the 
final  result  in  a  mixed  binary  format.  For  analog  operation,  several  important  issues  such 
as  t.he  dynamic  range  and  pi xe  1  - to-p i xe  1  uniformity  require  further  investigation. 

One  of  the  major  advantage  of  the  four  wave  mixing  configuration  is  the  availacility  of 
three  input  ports.  As  a  result,  the  system  is  sufficiently  flexible  that  it  car  be  adapted, 
with  little  or  no  mod i f icat ion ,  to  perform  other  types  of  processing  such  as  parallel  logic 
operations,  pattern  recognition  for  symbolic  subst  i  tut  ion ,  ^ and  residue  ar  i  th.T.c- 1  ic  compu¬ 
tation.  It  IS  quite  conceivable  that  the  incorporation  of  color  multiplexing,  '-.-.ich  adds 
another  cim.ension  to  t.-.is  approach  with  even  more  decrees  of  freedom,  will  equip  tee  system 
for  other  forms  of  image  or  data  processings. 
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Optical  Matrix -Matrix  Multiplication  with  N  Parallelism 
by  Spatial  Convolution  via  Four-Wave  Mixing 
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SUMMARY 

The  convolution  of  two  (two-dimensional)  spatial  patterns  using  nonlinear  opti¬ 
cal  four-wave  mixing  in  their  common  Fourier  plane  can  be  devised  to  achieve  matrix- 
matrix  multiplication  optically,  in  full  paralielism.  In  principle,  it  can  be  applied  to  the 
multiplication  of  a  matrix  A  (of  any  dimension  mxn)  and  a  compatible  matrix  B  (of 
diniension  nxp).  A  unique  feature  of  this  approach  is  that  the  output  beamlets  repre¬ 
senting  the  elements  of  the  product  matrix  lie  on  a  straight  line;  therefore,  a  linear  array 
of  detectors  (instead  of  a  two-dimensional  arra")  can  be  used  for  monitoring  the  result. 
The  major  advantage  of  this  scheme  is  the  full  parallelism  achieved  with  a  relatively 
simple  optica!  arrangement.  The  trade-off  is  the  resulting  spatial  complexity. 

Vlhite  and  Yariv^  have  demonstrated  that  (spatial)  convolution  and  correlation 
of  two  (two-dimensional)  patterns  can  be  achieved  in  real  time  by  four-wave  mixing  in 
the  common  Fourier  p'ane  of  the  input  patterns  arid  recording  the  phase-conjugated 
output  at  the  correspondin.  object  plane.  A  typical  experimental  configuration  is  shown 
schematically  in  Fig.  1.  Specifically,  if  U3  is  a  small  aperture  simulating  a  point  source 
(or  a  delta  function),  a  pattern  representing  the  convolution  of  U1  and  U2  (U1*U2)  is 
observed  in  the  output  plane.  An  illustrative  example  is  given  in  Fig.  2.  The  design  of 
U1  and  L’2  to  represent  two  matrices  so  that  their  product  is  represented  by  U1*U2  is 
explained  in  tfio  follov  mg  paragraphs.  Explicit  mathematical  analysis,  including  the 
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Fig.  1  A  schematic  diagram  illustrating  the  basic  idea  of  four-wave  mixing  in  a  non¬ 
linear  medium  located  at  the  common  Fourier  plane  of  the  input  spatial 
pattern. 
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U1  U2  U1  •  U2 

(a)  (b)  (c) 


Fig.  2  A  schematic  illustration  of  convolution  of  two  spatial  patterns:  (a)  and  (b)  are 
the  two  input  spatial  patterns;  (c)  represents  the  resulting  spatial  convolution 
of  the  two  inputs. 

effect  of  finite  aperture  dimensions  among  the  others,  of  a  similar  matrix  multiplication 

2 

concept  based  on  spatial  filtering  was  given  b\  R.A.  Heinz  et  al. 
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For  simplicity  and  clarity,  let  us  consider  the  simplest  cases  of  two  matrices  A 
and  B,  both  of  dimension  2x2,  and  their  product  C  =  AB  given  below. 


A  = 

^11 

ai2 

^21 

^22 

B  = 

t^ll 

bj2 

^21 

^22 

C  = 

AB  ^ 

^11 

*^12 

*^21 

C22 

(7) 


^11*^11  ^12^21  ^11^12  ^12^22 

^21^11  ^22*^21  ^21^21  ^22^22 


A  transparency  U1  consists  of  four  small  apertures,  each  with  its  intensity  transmittance 
proportional  to  each  of  the  matrix  element  ajj  is  shown  in  Fig.  3a.  A  similar  trans¬ 
parency  U2  corresponding  to  the  transpose  of  B  (i.e.,  rows  and  columns  interchanged)  is 
shown  in  Fig.  3b.  Note  that  the  vertical  distance  between  the  elements  in  U1  is  designed 
so  that  it  is  considerably  larger  than  that  in  U2,  while  the  horizontal  distances  between 
the  elements  are  identical  in  U1  and  U2.  Note  also  that  the  dashed  rectangular  box  con¬ 
necting  the  four  elements  in  U2  is  an  artifact  to  facilitate  the  following  explanation  and 
should  be  ignored  in  the  actual  design.  By  comparing  Figs.  2  and  3,  it  is  straightforward 
to  see  how  the  quadruplet  of  the  dashed  rectangular  box  is  formed  at  the  corresponding 
position  in  Fig.  3c  as  a  result  of  the  two-dimensional  spatial  convolution.  As  noted 
before,  all  the  dashed  lines  in  Fig.  3c  are  artifacts  and  the  actual  result  should  show'  only 
the  corners  of  each  box,  i.e.,  a  total  of  16  spots  with  the  intensity  of  each  proportional  to 
the  cross  product  of  all  the  elements  in  A  and  those  in  B.  The  four  doublets  lying  along 
the  Y-axis  in  Fig.  3c  are  drawn  slightly  offset  from  the  axis  to  expose  the  individual 
component.  In  practice,  the  two  components  of  each  doublet  are  spatially  overlapping  on 
the  Y-axis.  The  intensity  of  the  four  doublets,  from  top  to  bottom,  on  the  Y-axis  are 
proportional  to  Cj2)  C[j,  C22  and  C2[,  respectively.  The  eight  cross  terms  that  do  not 
contribute  to  the  matrix-matrix  multiplication  can  be  filtered  out  easily  as  they  are 
physically  separated  from  the  Y-axis. 
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Fig.  3  A  schematic  diagram  illustrating  optical  matrix-matrix  multiplication  by 
convolution  for  the  case  of  2  x  2  matrices:  (a)  matrix  A  and  the  mask  Uj  to 
encode  A;  (b)  matrix  B  and  the  mask  to  encode  B  ;  (c)  the  spatial  pattern 
resulting  from  convolution  of  U  i  and  U^. 

Another  example  illustrating  the  case  of  multiplying  a  2  x  3  matrix  with  a 
3  X  4  matrix  is  shown  in  Fig.  4.  Note  that  all  the  cross  terms  that  do  not  contribute  to 
the  matrix-matrix  product  are  omitted  in  Fig.  4c  for  the  sake  of  clarity. 

Using  an  argon  ion  laser  (5145A)  and  a  BaTiO^  crystal,  we  have  experimentally 
demonstrated  the  concept  described  above  for  the  case  of  2  x  2  matrices.  The  images  of 
the  masks  used  to  encode  the  two  matrices  A  and  B  are  shown  in  Figs.  5a  and  5b,  respec¬ 
tively.  The  experimental  result  representing  the  product  is  shown  in  Fig.  5c.  An  over¬ 
exposed  version  of  Fig.  5c  is  shown  in  Fig.  5d  to  reveal  the  noise  resulting  from  the  other 
cross  terms  that  do  not  contribute  to  the  matrix  product. 

In  principle,  the  basic  concept  described  above  can  be  appled  to  the  multiplica¬ 
tion  of  a  matrix  (of  any  dimension  mxn)  and  a  compatible  matrix  B  (of  dimension  nxp).  In 
practice,  the  signal-to-noise  (S/N)  ratio  is  expected  to  degrade  as  the  dimension  of  the 
matrix  increases.  The  problem  of  S/N  is  less  serious  in  the  muxed  binary  mode.  The 
major  advantage  of  this  scheii'e  is  the  full  parallelism  ac*^leved  with  relativey  simple 
optical  arrangenients.  From  the  technological  point  of  view,  the  fact  that  all  the  output 
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Fig.  4  An  illustration  similar  to  Fig.  3  for  the  case  when  the  dimensions  of  the 
matrices  A,  B  and  C  are  2  .  3,  3  *  4  and  2*4,  respectively. 
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Fig.  5  Experimental  results  for  optical  matrix-matrix  multiplication  by  convolution 
via  four-wave  mixing  in  the  spatial  frequency  domain,  (a)  and  (b)  are  the 
matrix  masks  used  for  encoding  the  matrices  A  and  (c)  shows  the  result  of 
their  spatial  convolution;  (d)  is  an  overexposed  version  of  (c)  to  reveal  the 
cross  terms  that  do  not  contribute  to  the  matrix  product. 


elements  lie  on  a  straight  line  can  be  significant,  since  this  allows  one  to  use  a  linear 
detector  array  instead  of  a  2-D  array  for  output  detection. 
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.-\  variety  of  optical  matrix  processing  schemes  [1,  2]  use  two-dimensional  spatial 
modulation  of  optical  intensity  to  represent  a  matnx  to  exploit  the  inherent  parallel  nature  vt 
optics.  Such  an  approach  tcpically  requires  the  projection  of  a  spatial  pattern  to  match  another 
pattern  to  perl’orm  the  element-by-element  multiplication.  The  basic  Lncoherent  matrix-vector 
multiplication  scheme  [3].  for  example,  requires  the  use  of  anamorphic  optics  to  project  a  line.ar 
array  of  sources  (or  a  one  dimensional  spatial  light  modulator)  to  precisely  match  a  two 
din'-ensiona!  m.atrix  masks.  In  the  matrix-vector  multiplication  scheme  using  four-wave  mixing 
in  nonlinear  media  [4],  simultaneous  alignment  of  all  the  pixels  of  the  matrix  and  vector  i>  ;; 
major  task,  paniculaily  for  a  large  num'oer  of  pixels.  Misalignment  of  the  pixels  may  lead  to 
severe  errors.  For  a  given  size  of  matrix  mask,  the  density  of  elements  increases  as  the 
dim.ension  (N)  of  the  matrix  increases.  As  a  result,  the  requirement  on  alignment  becomes 
more  and  more  stringent.  In  practice,  the  critical  alignment  required  is  likely  to  impose  a 
practical  limiit  on  the  optimum  dimension  of  the  matrix  (N)  to  be  of  the  order  of  one  hundred  or 
less  depending  on  the  specific  architecture. 

In  this  paper,  we  repon  a  new  scheme  for  optical  matrix-vector  multiplication  that  uses  a 
phase  conjugator  (with  a  finite  storage  time)  in  conjunction  with  a  spatial  light  rntdulator 
(SUM)  to  eliminate  the  pixel-by-pixel  alignment  requirement  at  the  cost  of  some  reduction  in 
parallelism  (N^/2  instead  of  N^).  Phase  aberration  due  to  imperfection  in  optics  is  also 
self-corrected  by  the  phase  conjugation  process.  The  optical  s}  tern  involved  is  relatively 
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Simple  compared  \>.ith  the  other  approaches.  Without  any  modification,  such  a  scheme  can  also 
perform  matrix-matrix  multiplication  with  N'^/2  parallelism. 


LASER  BEAM 


a 


Fig.l  A  schematic  diagram  illustrating  the  basic  concept  ct  optical  matrix-vectc^ 
.Tijitiplicaticn  using  a  photorefractive  phase  conjugator  in  conjunction  with  a 
scatia!  light  modulator. 


Re 


e  to  fig.l,  we  use  a  SIM  to  in. press  the  matrix  and  \'ec;or  information  in. 
input  laser  beam.  This  beam  is  directed  toward  a  phase  conjugator  which  h.;S  a 
me  (a  photorefractive  barium  titanate.  for  example).  A  cylind.ncai  lens  is  insened 


in  the  phase  conjugate  output  beam  path  to  perform  the  summation. 


Tr.e  principle  of  operation  is  as  follow,  the  SLM  first  impresses  the  matrix  informam. 
onto  the  input  laser  beam.  This  beam  is  then  incident  into  a  phase  conjugator  which  stores  me 
matrix  infomiation  after  a  finite  grating  formation  time.  When  the  matrix  information  m 
renm''\  ed  from  the  SLM,  say  by  turning  all  the  pixels  into  maximum  transmission  condition,  me 
pna-e  conjugate  beam  which  contains  the  reconstruction  of  the  matnx  infornaiion  exists  for 
finite  duration.  This  finite  storage  time  depends  on  the  strength  of  the  input  (read)  bc.-.m 
During  this  time,  if  the  next  frame  of  the  SLM  carries  the  vector  information,  par.-.i'.e. 
mui:;;  !;. .men  m  periormed  as  the  phase  conjugate  fvam  propagatCN  back  through  the  SLM 
L.er'.*  t.'.v  ',e‘,'t'’r  is  repreicnted  as  a  tvsi.t-uitiicimjtatal  areas'  ol  N  identical  uolunm  scciims.  ss  net; 
'■  ;  •  ne  tint  en- ton  ot  tiw  see  tern.  A  cylindrical  ietm  m  tiie  output  >n  i  ^  u  cc  t.e  jm:  ft  mm  t 
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sunimanon.  T^e  dark  storage  time  during  which  the  matrix  information  can  re  retrieved  :> 
detenrar.cd  h>  the  photorefractive  material  and  the  pumping  configuration.  I:  ranges  from 
'ceonds  to  microseconds. 

The  system  can  also  perform  matrix-matrix  mulriplication  by  time  multiplexing.  In  this 
case,  each  column  vector  (i=I  to  N;  which  constitutes  the  second  matnx  Mt  is  sequentially 
impressed  onto  the  beam  to  multiply  with  the  first  matrix  Mj  according  to  the  matrix-vector 
multiplication  scheme  described  above.  To  avoid  the  degradation  of  the  information  of  M, 
stored  in  the  photorefractive  hologram  during  the  readout,  it  is  necessary  to  refresh  the 
holographic  memory  with  Mj  to  restore  its  diffraction  efficiency.  This  can  be  done  oy 
re-impressing  .\1  j  onto  the  beam  after  each  readout  cycle.  Consequently,  a  total  of  2N  clods 
cycles,  consistmg  of  N  cycles  of  write  and  N  cycles  of  read,  will  be  required  to  carry  out  the 
multiplication  c:  two  NxN  matrices. 

Using  a  photorefraciive  barium  titanate  crystal  as  a  phase  conjugator  in  conjunction  uath 
a  -tSx-^S  magneto-optic  spatial  light  modulator  , SIGHT-MOD  SMD48I  from  Semetex  Core.'). 
V.  e  have  demor. ''mated  the  basic  principle  described  above.  Preliminary  experimental  results 
wil!  be  discussed.  Some  advantages  and  disadvantages  of  this  approach  will  be  compared  w'lth 
::u;'>e  >.  f  the  others. 

This  work  is  panially  supported  by  DARPA/AFOSR  contract  No.F49b20-S7-C-00i5. 
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Optical  interconnection  linking  laser  ar'ays  and  detector 
arrays  plays  a  key  role  in  opti-  al  computing.'  -  I Onceptual- 
iv,  such  an  interconnection  can  be  achieved  by  using  an 
optical  matrix  vecto.  multiplication. 

V  -  M\.  Ill 

where  V  is  the  in pu'  vector  representing  the sip.nt Is  carried  by 
an  arra\  of  .V  la'C-r'  and  v'  is  th‘-  output  vector  representing 
the  signals  received  bv  t  he  arra>  of  .\  detectors.  M  is  an  .\  x 
.V  ma.rix  representing  th'^-  interumnection  pattern 

When  a  triinspareniA  or  a  sp.it  lal  light  modulator '  SLM '  is 
used  a-  the  interi onnec t mn  ma-k.  parallel  matrix  vector 
multiplu  ate:  can  be  a(t.;ev>.ii  i.et  the  inlerci.n.ne.- tn-n 


mask  consist  of  an  .V  x  .V  array  of  transmission  windows. 
The  input  vector  (a  row  of  .V  elements)  is  fanned  out  by 
using,  for  example,  a  cylindrical  lens  pair,  so  that  the  light 
from  each  element  is  broadcast  over  a  corresponding  column 
of  the  window.-  .After  transmitting  through  the  intercon¬ 
nection  mask,  similar  optics  are  used  to  collect  the  light  from 
each  row  of  windows  and  to  sum  the ’output  into  a  column 
vector  of  .V  elements.  Such  an  architecture  provides  the  .V- 
parallelism.  However,  a  large  fraction  of  energy  is  absorbed 
by  the  transparency  or  SLM.  When  used  as  a  crossbar 
switch,  such  an  architecture  has  an  energy  efficiency  of  only 
1 ''.V,  where  .V  is  the  dimension  of  the  array.  '  This  occurs 
because  (.V  -  ll  '.V  of  the  light  energy  from  each  element  of 
the  input  vector  does  not  pass  through  the  crossbar  mask. 
The  energv  loss  increases  as  the  dimension  of  the  array 
increases.  This  is  sometimes  referred  to  as  the  fanout  ener¬ 
gy  loss.  For  a  1000  X  1000  crossbar  switch,  the  lo,  s  due  to 
fanout  can  be  as  high  as  99. 9T  This  is  not  acceptable  in  high 
speed  computing  because  signals  are  passing  through  the 
SLM  at  gigabit  rates  and  the  energx  loss  can  be  enormous. 
In  addition  to  the  inherent  tanoui  energy  loss,  all  SLMs  have 
finite  insertion  loss  dut'  to  imperfect  transmissions  such  as 
from  afisorptlon  and  scat  tering  It  we  include  such  insert  ion 
losses,  the  energx  effii  lem  for  a  c'ossbar  switch  would 
bfi  ome  t  X.  where  I  is  the  tran-mittance  if  <  D  of  each 
window 
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Fij;  !  >i.htmaiii  drawii'i:  n!  an  1  >.  ,V  (vpiical  iiiternT.neition 
usini;  dvr.amii  phol"rf!riU!i'>a  hulu^ram'  with  A  =  -1- 


K.a  -  Si  h'ln.i’.  ii  draw. Hi.'  an  \  •'''  .V  upti.a'  '.ntvrni.nnr-ati'.n 
ii-:n^  'i\  r.ainii  ph"i.  ir.n  h"‘.iaran:'  wiih  =  ; 


It  h.i'  been  km  'w  n.  !i  ir  s.  ntie  t  iine  t  hat  ah")'  >l:raI■lhi^  npt  w'a! 
eiernt-ni  '  Hb.'Ki  ean  be  enipiiwed  lur  tree  .'jtace  "iitieal  inter- 
ei 'nnec: n in.  Light  iruin  eat  h  la.^er  Miiirce  u!  the  in|>ut  array 
i>  Brat; a  >cat  tereci  and  redirected  tn  nne  nr  mure  deieiii'r>  >  a 
the  n'jtput  arra\.  There  are  i-eeeral  epeciln  reqi.iirenient- 
th.at  need  tn  be  met  inr  an  HtbF.  tn  he  u>ed  (nr  iiiterennnee- 
ti'in  iit'X'I,.''!  eirnnit'.  Ttu-'t-  include  near  prriect  tilignment 
and  high  diftrat  tii  -n  eft'ici eP.cv,  In  nririition.  a  new  hiilograir, 
1'  needed  tnreach  new  interciinnectinn  pattern.  In  view  c! 
t'ne.'e  ir'iie-,  i'  cie>irab!e  tn  use  dvnamic  hnlngraphic  media 
s'.K  h  ti'  photi  irelract  ive  crotal-  ’  Phntorelrat  tive  i  rv-ta;> 
sui  h  as  H.'iTiO.  and  SKX  art  b\  tar  the  nio.-t  etiuieni  re.ii- 
lime  hnb'graphic  meri.a.  Tnt-  rei  rdmg  nl  a  cnltinie  hnlc- 
gr;in,  in-trle  there  crv-t  lU  ref|iiire- tune  id  the  i-rder  ol  1  m> 
U'ing  light  inten.  ;ty  ul  tin-  ureter  id  1  N\'  i  m-,' 

In  thi-.  ( ''.'triiTiunicatKin.  we  pr'-pure  a  scheme  lUrming  a 
rivnnrnli  hubigrarn  insirie  a  ph' .ret ract ne  crystal  so  that 
the  hi 'log  ram  is  cap;iijl(-  ol  perlorming  rei  onfig'irable  optical 
intertop.neciii  r.  with  a  verv  high  energv  etliciencv  Such  a 
.'C  heme  utiii/er  the  it ’nrecinrocal  energy  trans'<-"  in  two- 
wa'.e  miM.'ig  to  ai  hie'.e  ati  exrrenielv  high  energe  eliicience 
Helerring  t"  Lig-  i  and  J.  we  den  ribe  ihm  new  method  o( 
ret  "Ill  igiir.ilile  opinal  intervonnet  lion  which  utilize-  tht 
n  .nret  iprot  a!  energe  tr,-.:i-ler  in  photi. retractive  two-wave 
mixing  to  improve  the  energv  ef'n  lent  v  Figure  I  describcr 
a  1  -  I.l  t  a=e  fi  ir  t  he  rake  oi  rlantv  in  explaining  the  tencept,  A 
-mall  !  rai  t  ion  of  thi  -  KT.  ifienl  la-f  r  beani  ir  t  oupled  out  of  the 
lieam  b\  uring  a  bean,  rplnter  i  HS  Tr.i-  sm;.'!  (rat  tion 
oabni  proiic  bt-anii  i-  lii!-n  ex|)anfl'-ti  b\  Uring  a  eyiindnea! 
len-  anu  p.C'-t--  thr  'iigti  I.’K  -pali.al  light  modulator  iSLMn 
I  r.  •in  (  xanipii-  -  how  t,  t  (.e  la-t-r  li-ap,  i-  f  -  be  ci  iiinected  to 
dt-tei  t'lr-  '  ami  a  a-  pre-'  ribed  Ic.  tin-  SLM.  The  iran-mit- 
ti-d  prob>.  tieam  ir  im-n  ri-i  "tnbim--!  with  th>-  mam  bean, 
in-itii-  a  ph'  .1 1  irt-t r.ii  1 1  \  t- 1  re  -  I  al  .-X-  a  re-uh  o!  iionrr  iprocal 
i-n.-rg'.  'o  ipbrig  ainc'-’  ali  ttm  em-rgv  in  ibt-  ni.'iin  boarn  i- 
I  r  1.-.  - ! <  r  r‘  fl  t . .  i  m-  pri  .1  a  O'  in.  w  Ii a  1;  o.i me-  I  la-  i r.t t-n  o: 
rm'iioj  p.itier;  I  i.a  ra-n''  ;-  .m  mien 'Ji.ta-it  i.  .n 

tipo,  I;  [  i  on''C vO  g  /  .  •! 


Figure  2  describes  the  reconfigurable  interconnection  for 
laser  arrays  and  detector  arrays  In  the  example  (a  -4  X  4 
interconnection  I  show  n,  laser  1  is  to  be  connected  to  detec¬ 
tors  h  and  c.  laser  2  is  to  be  connected  to  detectors  a  and  cl. 
laser  .3  is  to  be  connected  to  detectors  c  and  d.  and  laser  4  is  to 
be  connected  to  detectors  a  and  c.  In  terms  of  the  matrix- 
vector  multiplication,  such  an  interconection  can  be  written 
as 


where  i  ..  t  ..  i  and  14  are  signals  carried  by  the  laser- 
elements  of  the  arrav  r. 

A  cylindrical  lens  is  used  to  focus  the  2-1)  array  of  fieaii;- 
into  a  vector  (1-1)  array  1.  .As  a  result,  deiecinr  a  receives 
signals  from  lasers  2  and  4.  detector  h  receive-  signals  irom 
laser  1 .  detector  (  receiv e-  signals  Iroin  laser-  1 . :  .  and  4.  and 
detector  a  receives  signals  from  laser-  2  and  ..  ."sui  li  a 
concept  can  he  extended  to  in’ ert onnet t  .V  la.-t-r-  wn!.  .\.t 
detector-  where  .V  and  ,\/  are  f.vo  large  numbers. 

1  he  two-wavt  mixing  de-cribefi  in  Fig-  '  and  2  in.iv  lx 
viewed  as  a  real-time  holograi'lyv  in  which  the  rei'ordmg  and 
readout  occur  siniultaneousiv  inside  the  [ihotorei rai  t iv  c 
crv  .-tal.  The  beam  splitter  and  tht-  .''I.M  are  u-ed  to  rei  ord  .1 
volume  hologram  which  refireseni.-  the  iniereonnect loii  p.it- 
tern  as  iirescribed  by  the  SLM.  The  energv  coupling  in¬ 
volved  in  the  two-wave  mixing  ensure-  that  the  dittraction 
efluieiicy  during  the  readout  is  ainio.-t  loim.  This,  ot 
cour.-e.  requires  a  proper  (.'rientaiion  of  the  crvsial  -o  that  the 
energy  of  the  readout  beam  is  greatly  depleted.  The  high 
energy  efficiencv  result  -  from  the  I  act  that  mo-t  of  the  energy 
Is  carried  by  the  readout  beam,  w  hich  does  not  pass  through 
the  SLM  hut  is  difiracted  into  the  interconnectu  n  pattern 
by  the  hologram. 

The  energv  efficiencx'  of  such  an  interconnection  pattern 
can  he  estimated  tor  a  crossbar  switch  a-  follows  Let  h‘  lie 
the  reflectame  o)  the  beam  splitter.  It  is  legitimate  to  as¬ 
sume  that  the  tieam  splitter  is  iiraciicallv  los-ltr-s.  W'e  max' 
also  assume  that  tlie  surlCK  C-  of  the  photoreirat  lix  e  crx  sial  is 
antirenection-coated  so  that  the  Fresnel  reflection  loss  can 
he  neglected.  In  these  conditions,  the  two  beam-  thai  arrive 
at  the  (vhotorefractive  crx-ti-,!  haxe-  energie-'  '  1  -  ft'  and  ft; 
.V,  respiectivelx .  Inside  the  crx  sial.  the-e  txvo  beams  undx^rgo 
photoreiractive  coupling.  .As  a  result.  nTot  of  the  energx'  of 
the  pump  beam  1 1  -  fti  is  transtt  rred  to  the  prolie  beam  ft: 

.  xvhich  contains  the  interconnet  tior.  [uittern.  The  energx 
effix  lency  can  he  easily  derix  ed  and  is  given  h\ 

i.fl  1  -  .e. 

r  -  _  expi  -  i 1 .  1  ,- 

1  +  expi 

where  n,  i-  the  beam  intensity  ratio. 

I :  -  LiV 


and  /,  is  the  interaction  length.  '  is  tht  coupling  constant, 
and  o  is  the  bulk  absor|)t  ion  coef  fu  lent  f  or  photorcl  ract  ix  e 
crystal-  such  a-  HaTiO  anri  .'sBN.  the  cou))hng  lonstant  1- 
xerx  large  li.e,  y /,  >>  Ii  The  efluieiux  can  bt  xxritlen 
a[ii>ri.'Ximatt  lx 


W  e  not  III  that  lor  largi-  .tin-  energx  efiu  ii-iu  \  1-  linmed  b\ 
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h .  implif.ed  probe  imeasit>- !  verticaJ  5Cile  =  SOmv  div  honzonuJ  scale  =  l(i.sec  div . 

Kif  ■'  lai  Sit'nal  carried  ii'.  the  (ir.ihe  beair,  Wrticai  -ca)!  =  .■ 
rrA'  dn  .  hr  rizn:-,i;i!  >ra'i-  =  !  u<  ri:\  ■  ti'  Aniplitied  -iciia!  carried  t.v 
t  he  [jri  I  he  hea:;,  id :  er  anderL’'  ';r;c  t  s\  -  ■«  a\  e  mixihc  \  erticai  '( aie  = 

■'' '  ;t,\  di\.  huriziT.:  „a  =  1  ii>  (i;-. 


the  erC'tal  bulk  ati''  rjh ii>ii  e\f>'-.  /.  ■  anri  j>  niaximi/ed  lj\ 
i.i'inc  a  heatr, ']  h.:  it  r  witi.  a  vt-ry  email  refiec  tant  i-  Ii  u.e..  Ii 

(II. 

In  (ii.ir  prt-iiniip.ar'.'  experimenta!  wurk.  vet  inveeticated  i  he 
i"iie  I  ii  enerc.v  eiiicient.  \  and  the  cajialiilii  v  oi  hit'h  data  ra'e 
t  rar.'rni"ii  n.  In  examininc  the  enerev  e!iici(-n(  > .  ue  ti'eri  a 
la-er  beam  fnan  an  araun-inn  l.i-e-r  operalint'  at  'iH.h  nm. 
'khe  lii'er  iiean.  ie  (  ollim  'ted  int.ia  beam  nf  d-mm  diameir-r 
i);cii~im;a  leti'i'l  im  ai  lendt  h  ;' =  d  n..  A  beam  epblter  w  itn  a 
retledame  id  !i  =  n.i)"'  tieefi  tc  redirect  •’>'<  id  the  eiierc\ 
thruadh  the  vpfjtjal  heht  im 'ciuiatnr.  e-talti’ehins:  the  pri.he 
iiearn  The  remainder  ot  the  energy  tran-  nit'  thrmieh  thf- 
iieatii  'f.'htrer  a.nd  cinistituie.'  the  pump  beam,.  In  our  pre 
hnnnar'.  experiment,  the  SLM  wrt'  replaced  with  a  neutral 
den -it  >  finer  w  ;l  t;  a  \  anable  optic  a!  den~it>  t"  >1  mu  late  the 
lanii  n  <  nerd',  h.--  'I  he  twn  rieant'  inter'C-ct  in.'idea  liariuni 
lilatiaie  rx-tal  '.xitli  an  intf-ract  ioi,  lendth  of  .'.d  min  The 
[innpi  he  or,  v.  n  n  an  inten-if.  ot  |i  I  \\  r  m  ■  enter.'  the  crx  .'ta! 
a'  iie.i.i-  re  0  01,11  on  niem  e  'I  he  pr,  .be  b«-an:  CMI h  an  inteii'il  v 

Ii  '  W  I  n:  I'  ini  ider;  at  an  a.tidle  of  t"'  The  cr.V't.ii  m 

I  'ran.'  ei;  -  •  h.n  t  he  pr-  hr  beiin.  o  amp; it  led  'be  expeii'e  of 

if,e  p  ,n.[,  be. on  .Alter  ji,i--,nd  tliro  icp  jj,  ir-C'ta!.  thr- 
[c.n.p  t.e.oi.  I  -  ■■  ;rt  1,1  :i\  rlej,'.  ■  u  .hereri'  i  'a-  inten-m  of  thi 
pr,  t  ;e  ,  I  -  e  [I  ,n  I  d  i  ,  ,  ;  ' '  W  .  11.  'I  tK  -.all.-  e \ per ; 1 1 : e!  ' 

I  ■  ri  I n  1  !  n , .  'Id  la  ;-  r.i  d‘  t;- cc  n  t:  i  raii'init  tanc  e  oi  o  ] 

'!  r.i  -I  1  or.  t  ta  -on  inten-.-a  -  ot 


O.O."'  and  O.OO.h  VV/cm-,  respectively.  The  intensity  of  the 
amplified  probe  beam  remains  at  1.9  W/cm-'.  This  trans¬ 
lates  to  an  energy  efficiency  of  ~20T.  In  other  words,  if  the 
fanout  loss  is  99'’;  (e.g..  a  100  X  100  crossbar),  the  energy 
efficiency  of  our  new  scheme  can  be  at  least  20  times  better 
than  that  of  the  direct  approach.  Notice  that  our  crystal  is 
not  antireflection  coated,  and  "IST  of  energy  is  lost  at  the 
front  surface.  Bulk  absorption  in  this  particular  crystal 
accounts  for  ~60T  of  the  energy  loss. 

Even  though  the  reconfiguration  time  is  limited  by  the 
photorefractive  response  time,  which  is  typically  of  the  order 
of  milliseconds  at  modest  intensities,  the  photorefractive 
interconnection  system  can  accept  very  high  data  rate  sig¬ 
nals.  To  demonstrate  this  fact,  temporal  modulation  was 
impressed  on  a  laser  beam  (argon.  .\  =  .01-4..5  nml  using  an 
acfiustoopiic  device  to  simulate  a  signal,  which  is  to  be  inter¬ 
connected  with  some  output.  The  signal  used  was  a  pulse 
train  of  frequency /,  =  0.833  .MHz  with  each  pul>e  being  -'0.2 
uf-  wide.  This  rate  is  clearly  much  higher  than  the  reciprocal 
of  the  photorefractive  respoO'e  time.  The  modulated  laser 
beam  wa>  then  split  into  two  beams  and  mixed  it/  the  crx'.sta! 
as destrilied  before,  and  the  amplified  probe  beam  was  moni¬ 
tored  with  a  photodetector.  The  upper  oscilloscope  trace  in 
Fig.  3  show>  the  input  probe  signal,  and  the  lower  one  shoW' 
the  amplified  profie  signal.  The  re'ult^  show  a  steady-state 
re'poii'e  in  which  thi-  temporaiix  niodui.iled  [uimp  and 
probe  beam.'  interact  simply  iiy  rii:  iracimg  ,  31  the  stationary 
index  grating  that  is  created  in  the  cry>tal  after  the  photore¬ 
fractive  re^poit'e  time.  .As  this  experiment  wa.s  performed 
merely  to  demonstrate  the  high  signal  iiandwidth  of  the 
system,  optimization  of  the  parameter'  was  not  done,  and 
the  results  in  Fig.  3  can  clearly  lie  imitroved. 

Summarizing:  we  have  descriiieti  a  new  method  of  recon - 
figuralile  optical  interconnection  using  dynamic  holograms 
in  photorefractive  crvstals.  Such  a  new  method  provides  a 
very  high  energy  efficience .  The  interconnection  can  be 
reconfigured  in  using  a  different  SLM  pattern.  The  holo¬ 
gram  formation  time  will  limit  the  rec nnfiguratlon  time. 
Once  the  hologram  which  contains  the  interconnection  pat¬ 
tern  is  formed  inside  the  photorefractive  crystal,  such  a 
scheme  is  caftahle  of  providing  the  interconnection  for  high 
data  rate  transmission.  In  such  an  interconnection,  the 
out|)Ui  of  each  laser  can  lie  inpui  to  any  one  or  all  the 
detectors.  Optical  phase  conjugation  can  also  be  used  in 
conjunction  with  the  iwo-wa\e  mixing  to  correct  for  any 
phase  aberration  that  may  be  causc  d  tiy  the  crystal  imperfec¬ 
tions.' 
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Abstract.  This  paper  describes  various  nonlinear  optical  phenomena  in 
photorefractive  media  and  selected  applications  in  optical  computing.  These 
phenomena  include  optical  phase  conjugation,  two-  and  four-wave  mix¬ 
ing,  and  real-time  holography.  The  applications  include  image  amplifi¬ 
cation  and  subtraction,  logic  and  matrix  operations,  and  optical  intercon¬ 
nection. 

Subject  terms:  optical  computing:  optical  interconnections:  optical  image  pro¬ 
cessing,  nonlinear  optics. 
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1.  INTRODI  CTION 

The  held  of  optics,  speciticallx  nonlinear  optics,  is  cmcrcinp  as 
^  an  area  ot  increasine  imporiance  tor  applications  in  optical  coni- 
putinp  because  ot  its  inherent  parallelism  and  its  extreme  hand 
uidtti  ot  communication  '  In  addition,  the  nonlinc.ir  optical  iNI.O’ 
response  ot  materials  alloxxs  |or  the  possibihtx  ot  manipulatinL' 


tfi.il'.e  l'.i;s  :  I  II  1  i  i  re.',:  .1-.!  ( >,  :  V  t'l.s'  r-..-.  ih.iniis^  rip-  re.cOL-a  la?. 
13  1  US',  ..—  ‘  I'C  piO'll,  111,  ■■  J  il’  1  ‘o  . 

'  lus'j  S. t’lv  ,t, ,  <  Jp’:.  .11  inanjna-iii.iti,'!  t  neni,.--.  r- 


the  propagation  of  an  optical  beam  bx  another  beam  of  light 
inside  a  nonlinear  medium.  These  manipulations  include  changes 
of  phase,  polan/alion.  iniensitx .  direction  of  propagation,  and 
even  the  frequencx  of  light  itself."’  In  the  area  ot  optical  com¬ 
puting.  these  NLO  manipulations  ofter  unique  capabilities  in 
holographic  interconnection,  parallel  image  subtraction,  phase 
c(xniugatK>n.  optical  data  storage,  spatial  light  modulations,  and 
bistable  optical  dexices. 

Some  of  the  most  useful  NLO  phenomena  include  optical 
phase  coniugation.'*  four-xxaxc  mixing.'  and  ixxo-xxaxe  mixing/’ 
Traditionallx .  these  phenomena  are  obserxed  in  nonlinear  media 
that  exhibit  a  large  third-order  ruscepubtliix .  Txpicallx.  these 
phenomena  require  xerx  high  optical  intensities  for  efhcient  op- 
cratKxn.  L'sing  one  of  the  most  efhcient  NLO  materials  such  as 
carbon  disullidc.  an  optical  intensitx  lexel  on  the  order  of  1  M\\' 
enr  IS  still  needed  These  intensity  lexels  are  too  high  for  optical 
computing  applications  Photorefractixe  matenals  such  as  BaTiO.. 
SBN.  BSO.  and  BOO  arc  by  far  the  most  efhcient  media  for 
these  NLO  phenomena  using  relattxelx  loxx  intensitx  lexeK  te  g  . 
I  cm"i.  In  addition,  these  materials  are  acentric  and  exhibit 
nonrcciprocal  energy  transfer  in  a  degenerate  txxo-beam  coupling 
configuration.  Such  optical  nonrecipriK'itx  is  not  axailable  in  a 
third-order  nonlinearity 

in  this  paper,  xxe  rexiexx  some  of  the  mt'st  usetui  NLxT  phe¬ 
nomena  in  photorefractixe  media  We  then  discuss  their  appli¬ 
cations  in  optical  computing 

2.  PHOTORFFRACTIN  K  M ATF.RLALS 

The  photoretractixe  eftecl  is  a  phenomenon  in  xxhich  the  R'cal 
index  (>t  retraction  ts  changed  by  a  spatial  xtiriation  o!  the  light 
intensity  Such  an  eftecl  xxas  first  repi'rted  in  |h(ib  The  spatial 
inxtex  xanation  leasK  to  distortion  ot  the  xxaxetront,  and  such  an 
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effect  was  first  referred  to  as  'optica!  damage."  The  photore¬ 
fractive  effect  has  since  been  observed  in  many  electro-optic 
crystals  including  LiNbOi.  BaTiO^.  SBN,  BSO.  BGO.  GaAs. 
InP,  etc.  It  is  generally  believed  that  the  photorefractive  effect 
arises  from  optically  generated  charge  earners  that  migrate  when 
the  crystal  is  exposed  to  a  spatially  varying  pattern  of  illumi¬ 
nation  with  photons  having  sufficient  energv.  Migration  of  the 
charge  earners  due  to  dnft  or  diffusion  produces  a  space-charge 
separation,  which  then  gives  nse  to  a  strong  space-charge  field. 
Such  a  field  induces  a  refractive  index  change  via  the  Pockels 
effect.  This  simple  picture  of  the  photorefractive  effect  explains 
several  interesting  steady-state  optical  phenomena  in  these 
media. 

Althoueh  there  are  several  models  for  the  photorefractive 
effect.*”''  the  Kuk;  tarev  \  inetskii  model  is  the  most  widely 
accepted  one  *  In  this  model,  the  photorelractne  materials  are 
assumed  to  contain  donor  and  acceptor  traps  that  anse  from 
imperfections  in  the  cry  stal  These  traps  create  intermediate  elec¬ 
tronic  energy  states  in  the  bandgap  of  the  insulators.  When 
photons  with  sufficient  energy  are  present,  electronic  transitions 
due  to  photocxcitation  take  place.  .As  a  result  of  the  transitions, 
charge  earners  are  e.xcited  into  the  conduction  band  and  the 
ionized  donors  become  empty  trap  sites.  The  rate  of  carrier 
generation  is  isl  PkN'd  -  N'u  i.  whereas  the  rate  of  trap 
capture  is  ^rNN'u  •  Here,  s  is  the  cross  section  ot  photoion- 
ization.  I  IS  the  lighi  intensity .  p  is  the  rale  ot  thermal  generation. 
"VR  is  the  camer-ionized  trap  recombination  rate,  and  N  and 
N’u  stand  for  the  concentration  of  the  carriers  and  ionized  traps. 
No  IS  the  number  density  of  the  donor, 

Th.^  .,pacc-^l.ajgs  field  p.^>duv.sd  by  tne  migration  of  the  charge 
carriers  is  determined  by  the  follow  mg  scl  ot  equations* 

f’  1  ^ 

-N  =  -Nf.  -  J  . 

m  rtl  c 

-Sr  =  (sl  ~  Si  N-,  -  N;'.  ■  --  vkNNl'  . 

J  =  eAN'j  F.  -  ^VlogN  )  -  pi:  .  1.0 

V  (fE'  I  =  eiNA  -  N  -  Nf,  I  .  i-li 

where  c  is  the  unit  vector  along  the  c-avis  of  the  crystal.  N.^  is 
the  accepior  concentrainm.  p.  is  the  mobility.  T  is  the  absolute 
temperature,  k  is  Boltzmann's  constant,  n  is  the  index  of  re- 
tractii'n.  t  is  the  di  '  tru  tensor,  pi  is  the  photovoltaic  current, 
e  IS  the  electronic  coarge.  and  p  is  the  photovoltaic  consiant 
E''  stands  tor  the  space-charge  field  E  is  the  total  field,  which 
includes  E"  and  any  external  or  internal  fields  (such  as  chemical 
or  interna!  terri'clectnc  tielJsi 

As  a  result  of  the  presence  ol  the  space-charge  field,  a  change 
in  the  index  (fi  rctractuin  is  induced  by  means  ot  the  linear 
elcvtro-optis  eftCsf  ( fAickeh  ettcM  i 

a(  -'  1  -  r..K.  . 

rr 

where  r,.i  is  the  elesiro-opiu  ^oel  tic  lent  iwitti  i.t.k  ■  x.y.zi. 

.V  Nf»M.INK\R  OPTK  AI  PHENOMENA 

Ail  o;  c.e  M.t)  pherioni-cn,-.  rienPone.l  e.ulier  invlxe  the  tor 
ma’ion  I ’t  xoluiiis'  iriuCx  craiir.;-  iiisisiv'  itis-  ph..  torelr.istive  me¬ 


dium.  The  simplest  way  to  produce  a  volume  index  grating  is 
to  employ  two-beam  interference  inside  a  nonlinear  medium. 
This  process  takes  advantage  of  the  nonlinear  response  of  the 
material  to  illumination  by  electromagnetic  radiation.  Thus,  when 
two  beams  of  coherent  light  intersect  inside  a  photorefractive 
medium,  they  create  a  volume  index  grating.  When  the  two 
beams  propagate  through  the  self-induced  grating,  they  undergo 
Bragg  scattering.  One  beam  scatters  into  the  other  and  vice  versa 
in  a  process  termed  two- wave  mixing.  The  hologram  formed  by 
the  two-beam  interference  inside  the  photorefractive  medium 
can  also  be  erased  by  uniform  light  illumination.  Thus,  dynamic 
holography  is  possible  by  using  photorefractive  materials 

The  Bragg  scattering  involved  in  two-wave  mixing  is  very 
similar  to  the  readout  process  in  holography .  In  four-wave  mix¬ 
ing.  a  third  beam  is  used  to  read  the  hologram  formed  by  the 
two-beam  interference.  The  fourth  beam  is  generated  as  a  result 
of  the  Bragg  scattering.  To  satisfy  the  Bragg  condition,  the  third 
beam  must  be  counterpropagating  relative  to  one  of  the  two 
beams  that  are  involved  in  the  formation  of  the  index  grating 
In  two-wave  mixing,  the  Bragg  condition  is  automatically  sat¬ 
isfied.  In  w  hat  follow  s,  we  discuss  some  of  the  NLO  phenomena 
in  photorefractive  media  that  are  useful  for  optical  computing 
applications 

3.1.  Two-wave  mixing 

Two-wave  mixing  is  important  in  practically  all  of  the  useful 
NLO  phenomena  observed  in  photorefractive  media  The  process 
of  torming  an  index  v  ariation  pattern  inside  a  nonlinear  medium 
using  two-beam  interference  is  very  similar  to  that  of  hologram 
fomiation.  Such  an  index  vanation  is  often  penodic  and  is  called 
a  volume  grating  in  addition  to  the  holographic  prixress  known 
traditionally,  beam  coupling  occurs  simultaneously.  This  is  a 
unique  property  of  photorefractive  matenaN.  In  what  follows, 
we  briefly  review  the  coupling  and  energy  transfer  of  two  beams 
inside  a  photorefractive  medium. 

Consider  the  interaction  of  two  laser  beams  inside  a  pnoto- 
refractive  medium  (Eig,  li.  If  the  two  beams  are  of  the  same 
frequency,  a  .siaiionary  interference  pattern  is  formed.  Let  the 
electric  field  of  the  two  waves  be  wntten 


t:  --  Aevpliiiat  -  k -rij  .  j  =  ).2  ,  ifi' 

where  ,Ai..A:  are  the  wave  amplitudes,  w  is  the  angular  fre¬ 
quency.  and  ki.k;  are  the  wave  vectors  For  simplicity,  we 
assume  that  both  beams  are  polarized  perpendicular  to  the  plane 
of  incidence  ti  e  .  s-polarizedi 

V^lthln  a  scaling  factor,  the  intensitv  can  be  wntten 

t  L  •  =  L:  -  L;  '  l~l 

Lsing  Eq  tbl  tor  the  elcvtrk  field,  this  can  be  written 

I  -  .A,'  -  .A-'  *  .ALA.-exp:  -  iK  r-  -  .A;.A’evp' iK  r^  .  tSi 

w  here 

K  k  -  k.  (Sti 

and  K  2-  \.  where  \  is  the  p-criod  ot  the  fringe  pattern 
The  intensitv  |I.q  (ki|  h.i-  a  sp.itial  vari.ilion  inside  the  photo 
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refractive  medium  According  to  Kukhtarev's  model,  such  an 
intensity  pattern  will  generate  and  redistnbute  photocarriers.  cre¬ 
ating  a  space-charge  held  in  the  medium.  This  held  induces  a 
volume  index  grating  due  to  the  Pockels  effect.  In  general,  the 
index  grating  will  have  a  spatial  phase  shift  relative  to  the  in¬ 
terference  pattern.*'  The  index  of  refraction,  including  the  fun¬ 
damental  component  of  the  intensity -induced  gratings,  can  be 
written 


*2 


PHOTOREFflACTIVE  MEDtOV 


He. 


ni  ,.,.ArA: 
-^exp(iil)  1  — j — exp(  ■ 


iK  ri  c  c. 


Fig  1.  Schematic  diagram  illustrating  the  basic  idea  of  photorefrac- 
tive  two-wave  mixing. 


where 

i(i  =  Ii  i:  =  j.Ai'  'A-'  .  (Ill 

rui  IS  the  index  of  refraction  when  no  light  is  present,  (f)  is  real, 
and  ni  is  a  real  and  positive  number.  For  the  sake  of  simplicity . 
we  assume  a  scalar  grating.  The  phase  6  indicates  the  degree 
to  which  the  index  grating  is  shifted  spatially  relative  to  the 
intensity  pattern.  In  photoretractive  media  that  operate  by  dif¬ 
fusion  only  li  e.,  no  external  static  heldi.  for  example.  BaTiO;. 
the  magnitude  of  6  is  t:  2  w  ith  its  sign  depending  on  the  direction 
of  the  c-axis.  K  is  the  grating  wave  vector  gisen  by  Eq.  (9i 
The  parameter  m  depends  on  the  grating  spacing  and  direction 
and  on  the  material  properties  of  the  crystal,  e  g.,  the  electro¬ 
optic  coefficient.  Expressions  for  niexp(i6l  can  be  found  in 
Refs.  9  through  1 1 . 

The  spatial  phase  shift  between  the  interference  pattern  and 
the  induced  volume  index  grating  has  been  known  for  some 
time.*'  ''  The  phase  shitf  allows  for  nonreciprocal  steady-state 
transfer  of  energy  between  the  beams  To  investigate  the 

coupling,  we  substitute  Eq  (  Uli  lor  the  index  of  refraction  and 
E  =  El  *  E:  for  the  electric  held  into  the  wave  equation.  This 
leads  to  a  set  of  coupled  equations. 

The  si'lutions  lor  the  intensities  Ipzi  and  btzi  arc' 


I  -  m  c\p’y/‘ 

I  -  IT. 

I-(/i  =  1-(I); - e\pi-u/i  .  il.'i 

1  *  mexp'  -  y/ > 

where  i  is  normal  to  the  crystal  surfaces  and  m  is  the  input 
intensits  ratio 


In  the  absence  of  absorption  la  =  Oi,  lyzi  is  an  increasing 
function  of  /  and  Iq/I  is  a  decreasing  function  of  /.  proxided 
7  IS  positise  The  sign  of  y  depends  on  the  direction  of  the 
c-axis  As  the  result  of  the  coupling  for  y  >  fl  in  Fig.  1,  beam 
2  gains  energy  from  beam  I.  It  this  two-wa\e  mixing  gam  is 
large  enough  t(>  osercome  the  absorption  loss,  then  beam  2  is 
amplified  Such  an  amplification  is  responsible  tor  the  fanning, 
stimulated  scattering  and  oscillation  ot  lasc  beams  m  photore- 
fractisc  crystals"' 

Vk  e  mentioned  the  holographic  interpretation  of  two-wa\e 
mixing  in  photoretractise  media  earlier  Elaborating,  we  con¬ 
sider  the  formation  ot  an  index  pr.iiing  due  to  the  presence  ot 
two  coheren!  laser  beams  inside  a  photoretras ti\ e  crystal  This 


is  formally  analogous  to  the  recording  process  in  conventional 
holography .  Consider  Fig.  1 .  in  which  two  laser  beams  intersect 
and  fonti  an  induced  index  grating.  The  index  grating  [Eq.  (10)] 
contains  the  product  of  the  amplitudes  Ai  and  Ae.  This  grating 
is  a  hologram  formed  by  an  "object"  beam  A i  and  a  "reference" 
beam  Ae-  The  diffraction  component  of  the  transmission  function 
of  such  a  hologram  |i  =  expt  -  lAnL  A)]  is  given  approximately 
by 

I  -  An  -  ACA.'Cxp'  -iK  r)  -  Ai.A'expiiK-r)  .  (15i 

where  .Aj  and  A;  denote  the  complex  amplitudes  of  the  object 
and  reference  fields,  respeciixely ,  Equation  ll5l  assumes  that 
the  modulation  is  weak,  so  the  higher  order  terms  can  be  ne¬ 
glected. 

For  reconstruction  (see  Fig.  1 1.  the  hologram  is  illuminated 
b;  the  reference  beam  A:expi  -  ik;T).  The  diffracted  beam  can 
be  written 

T|.A!.A;.A:expi  -  ik;T  I  .  ilbi 

where  is  the  diffraction  efficiency  .  Notice  that  the  phase  of 
kz  cancels  out  and  the  diffracted  beam  is  a  reconstruction  of 
the  obicet  beam  .Aicxpi  -  ikon.  Similarly,  the  reference  beam 
,A:  can  be  reconstructed  by  illuminating  the  hologram  with  object 
beam  .Ai.  provided  beam  Ai  is  a  phase  object  (i.e..  A|  has  only 
phase  variation,  with  I.Ag  =  constant). 

In  addition  to  holographic  analogy .  two-wa\e  mixing  exhibits 
ampliticaiion.  which  is  a  unique  feature  not  available  in  con¬ 
ventional  holography.  Using  these  two  properties,  two-waxe 
mixing  can  be  used  for  beam  processing.  As  a  result  of  its  real¬ 
time  holographic  nature,  photorefractixe  two-waxe  mixing  ex- 
hibiis  nonrcciprocal  enerev  transfer  without  anx  phase  cross 
talk'" 

The  lack  of  phase  cross  talk  can  be  understood  in  terms  of 
the  diflraction  from  the  self-induced  index  grating  in  the  pht>- 
torelractixc  crystal.  Normally,  il  a  beam  that  contains  phase 
mlormaiion  ilxr.ti  is  diffracied  from  a  fixed  grating,  the  same 
phase  intormaiion  appears  in  the  diffracted  beam  For  a  self- 
induced  index  grating,  the  phase  information  il/ir.ti  is  impressed 
onK'  the  grating  in  such  a  way  that  diffraction  will  be  accom¬ 
panied  by  a  phase  shift  -ij/ir.ti.  Such  a  dynamic  hologram 
makes  scil-canccllation  of  phase  mlonriation  possible  when  the 
incident  beam  is  diftracted  from  the  grating  produced  bx  the 
incident  and  reference  beams  This  self-cancellation  ot  phase 
into'manon  is  actually  equivalent  to  the  reconstruction  of  the 
relerence  beam  when  the  hologram  is  reau  xiut  by  the  object 
beam 

F.nergx  iransler  xxiihixiit  phase  cn'ss  talk  can  be  employed  to 
compress  both  the  sp.iii.il  and  temporal  s|x'cua  c't  a  light  beam/'' 
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Fig  2.  Schematic  diagram  illustrating  the  basic  idea  of  photorefrac- 
tive  four-wave  mixing. 

In  other  words,  the  energy  transfer  without  phase  cross  talk  can 
be  utilized  to  clean  up  the  wavefront  of  a  laser  beam.  This  has 
been  demonstrated  expenmenially  using  a  photorefractivc  SB.N 
crystal''^'"''  to  correct  for  the  spatial-wavefront  and  temporal- 
wavefront  aberrations 

3.2.  Phase  conjugation  and  real-time  holography 

Referring  to  Fig  2.  consider  optical  four-wasc  mixing  in  a 
photorefractivc  medium  Beams  .^i  and  .A;  incident  from  the 
left  form  a  volume  index  grating  inside  the  medium  If  the  third 
beam  .A?  is  exactly  counterpropagating  relative  tc>  .A;,  then  a 
fourth  beam  is  generated  thai  pri'pagatcs  backward  relative  to 
.Ai  and  IS  a  time-reversed  replica  ot  that  beam.  This  is  what 
happens  m  the  readout  process  of  conventional  holography.'* 
except  that  in  Fvur-wave  mixing  the  writing  and  readout  pro¬ 
cesses  occur  simultaneously  Fv>r  this  reason,  phase  conjugation 
by  means  (if  optical  Uiur-wave  mixing  is  often  referred  to  as 
real-time  holographs 

If  we  consider  the  coml'inatu'n  (if  the  phois'iefractive  medium 
and  the  twci  counterpri'pagating  beams  i  .A;  and  .A-. )  as  an  optica! 
device,  then  such  a  device  behaves  like  a  retroreflective  mirror, 
,An  incident  beam  will  be  pha-e  conjugated  and  propagate  back¬ 
ward,  retracing  its  path,  Su  .in  optica!  device  is  called  a  phase- 
coniugate  mirror  (PC'Mi  ketiection  friim  a  PCM  exhibits  many 
interesting  properties,  including  wavefront  aberration  correc¬ 
tion.  piilari/alion  restoratuin."  ''  and  amplilicaiion,'  The  am- 
plitication  Is  possible  because  energy  can  he  coupled  trvim  the 
twii  counterpropagating  beams, 

.A  class  ot  phase  coniugators  that  have  received  considerab[e 
attention  recently  are  the  selt-pumped  phase  coniugators. 

In  these  coniugators.  there  are  no  externally  supplied  counter- 
propaeating  pump  beams  Thus,  no  critical  alignment  is  re¬ 
quired  The  phase-coniugaie  retlectiv  iiy  is  relatively  high  at  even 
l(iw  laser  power  These  devices  are  by  far  the  most  convenient 
PCMs  available  .Although  several  models  have  been  developed 
tor  self-pumped  phase  ci'niugators.'  the  phenomena  can  he 
easilv  understood  with  the  resonator  mode!  In  the  iwo- 

interaction  reeion  model.'  two  separate  tour-wave  mixings  are 
responsible  tor  the  phase  conjugation  In  the  2-K  grating  mi'del. 
selt-pumped  phase  coniug.ition  is  considered  as  a  priKCss  similar 
t(i  stimulate  '  Hnliouui  scattering  'SBSi  in  the  reson.itor  model, 
the  crvst.il  mibe  is  viewed  a-  an  opti.al  cavitv  that  supports  a 
multi'ude  ot  mode''  These  rimdes  are  trapped  inside  the  crystal 
owinc  to  tot.i!  interntil  rellection  at  the  surface 

When  .1  la-er  beam  en'ers  the  crvst.il.  sorne  ot  the  modes  m.iv 
be  exsited  .r-  a  resuli  ot  the  strone  par.iinetrR  gam  due  ti>  two 
w.K-e  mixine  In  p.iilkiii.ir.  ring  o-.mllaliiui  c.ii:  be  gener.iled 
ac.or  iine  t'  ,i  Itieorv  deveioped  in  Ret  '  '  W  iien  the  contie 
iir.it  ii  e,  .  if  Ine  re -on.in..  e  c  jv  I'v  rel.it  i v  e  to  (tic  iiKident  l.i-er  be.im 


supports  bidirectional  oscillation,  a  phase-conjugate  beam  is  gen¬ 
erated  by  means  of  four- wave  mixing. 

According  to  this  theory  the  frequency  of  oscillation  inside 
the  crystal  can  be  slightly  detuned  from  that  of  the  pump  beam. 
If  (u  is  the  frequency  of  the  incident  laser  beam,  the  frequency 
of  the  internal  oscillation  can  be  wntten 

w  =  w  -I-  b  .  (17) 

w  here  8  is  the  frequency  detuning  and  is  of  the  order  of  ±  I  Hz 
for  BaTiOx.  This  frequency  detuning  depends  upon  the  path 
length  of  the  nng  oscillation  inside  the  crystal.  The  bidirectional 
oscillation  provides  the  counterpropagating  pumping  beams.  By 
conservation  of  energy,  the  phase-conjugated  beam  has  a  fre¬ 
quency  of  CO  -*-  28. 

The  resonator  model  presents  a  simple  explanation  of  the 
frequency  shift  observed  in  BaTiO-.  self-pumped  phase  conju- 
gators  Experimental  evidence  indicates  that  internal  os¬ 

cillations  play  a  key  role  in  the  generation  of  phase-conjugate 
waves  ’ 

3.3.  Phase-conjugate  interferometry 

A  phase-conjugate  interl'erometer  is  simply  an  interferometer 
that  u>es  one  or  more  PCMs,  Several  contigurationv  of  phase- 
conjugate  interferometers  have  been  investigated  These  include 
a  Mach-Zehnder  interferometer  incorporating  a  phave  conjugator 
for  sell-referencing  phase  contounng  based  upon  the  interference 
between  an  incident  wave  and  its  phase-conjugate  counter¬ 
part,'^"'  a  .Michelson  interferometer  with  one  of  the  mirrors 
replaced  by  a  phase  conjugator."^"'  and  a  Michelson  interfer- 
(nneter  with  both  mirrors  replaced  by  one  or  two  phase  conju- 
gat(>rs. Several  optical  computing  applications,  including 
coherent  image  subtraction.'"^”"  "excluvive  OR"  logic  opera¬ 
tion. and  spatial  differentiation."  have  been  demonstrated 
Using  the  .Michelson  conliguration  with  both  mirrors  replaced 
by  a  single  self-pumped  BaTiO.i  phase  conjugator  (see  Sec,  3.2 1. 
.A  similar  conliguration  for  temporal  differentiation  lor  novelty 
tillering  I  has  also  been  conceived.'*''  In  this  section,  we  bnefly 
review  the  unique  features  of  a  phase-conjugate  Michelson  in- 
tenerometer  (F’CMli.  .Applications  in  optical  computing  are  dis- 
cuvsed  in  Sec  4.2 

One  of  the  most  important  and  interesting  phenomena  asso¬ 
ciated  with  PCMIv  IS  ri/ni’  reversuL  Referring  t('  Fig  3(ai.  con¬ 
sider  a  Michelson  interferometer  with  PC.Ms  .An  incident  laser 
beam  (with  intensity  1  ,i  is  split  by  beamsplit  er  BS  (assuming 
Kisslevsi  into  two  compiments.  which  are  retrorefiected  back  bv 
the  PC.Ms.  .According  to  Stokes'  principle  of  time  reversilulity  ." 
when  the  two  phase-con|ugate  beams  recombine  at  the  beam¬ 
splitter.  a  time-reversed  replica  of  the  incident  beam  is  generated 
This  beam  propagatcv  backward,  retracing  the  original  incident 
beam  path.  Thus,  there  is  no  light  at  the  output  piirt  (pi'rf  .A); 
the  total  darkness  is  a  result  of  time  reversal  Such  an  optical 
time  reversal  has  been  obser.ed  experimentally  using  twii  eiui- 
pled  BaTiO-.  crystals  tor  the  PCMs 

Referring  to  fug.  3ibl.  (igni'nng  transparencie-  Tqx.yi  and 
T'lx.vi,  the  etlects  id  which  will  be  discussed  in  Se^  4.2].  let 
t.r  and  t  .  r  be  the  amplitude  tratisinissuni  and  rellection  luk'I- 
fkients  lit  beamsplitter  HS  tor  beam  entering  tri'iii  the  right  and 
left,  respectively  T..^  input  beam  with  eomplex  liekt  .implitude 
I  is  split  bv  the  beamsplittei  inio  a  reilected  i.omponeni  rp.  and 
a  tr.iromitied  component  tl.  T  hese  f,' o  ^aimponents  arc  phase 
Loiiuieated  to  give  r'l  '  and  t'F.'  and  are  :  'nher  split  and  re- 

fiPTl,  k,  e  A;--  -Uiv  :s  \  4  3X1 
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Fig.  3.  (a)  Time  reversal  for  a  Michelson  interferometer  equipped  wrth 
PCMs,  (b)  Quantitative  analysis  of  the  time  reversal  of  a  PCMI. 

combined  with  r’E'r  *  fLi'f  at  port  B  and  r'E't  -  t*E’‘r' 
at  pon  A  The  principle  ol  re%ersibilit_\  requires  that 


( r  ■  '  fi'iE’  =  h"  ,  I  IKi 

nr*  ^  fr  lE'  =  d  .  (I9i 

Thus,  the  amplitude  cirefticients  satislv  the  relations 

j-  '1*1  =  I  .  Oil. 

tr*  -  I’r  =  (t  i21 1 


Equation  (2  1  i  establiNhcv  the  condition  lor  complete  cancellation 
of  the  held  at  pr'il  .A.  These  equation^  are  known  as  StoKe^' 
equations. 

In  practice,  the  cancellation  at  pLirt  .A  max  n<'t  be  complete 
bccuuxc  one  cannot  obtain  an  ideal  pha^e  coniueaiu'n.  The  mea- 
surement  of  time-rexerxal  lidelitx  umiil'  a  PCMI  wax  reported 
rcuentlx  bx  Eubank  and  \  a/que/  '  The  phase  instabilitx  arisine 
tr.'ni  the  uxe  o*  txxo  indepenJeni  ini'l  coupledi  self-pumped  phase 
coniUL’.iior'  IS  discussed  in  Ret  4,' 

4.  OPTK  AT  ( OMPl  TINf,  APPEIC.ATIONS 

The  Nl.ft  phenomena,  includine  two-  and  lour-waxe  mixine. 
optisal  phase  soniLicMtion.  anti  ph.isc  coniueale  interlerornetrx  . 
dessribed  in  the  prexious  setii.m  can  be  emploxeJ  m  a  wide 
ranc'  ot  aprTicalions  ti.r  optk.ii  eoniputinc  and  im.iee  proeess 
ire  A  Icxx  examples  are  eixen  in  this  scstion  lo  illustrate  some 
('I  the  men;-  an.l  limitations  of  these  Nl.<)  approaehcs 


Fig  4  Experimental  configuration  and  results  of  image  amplifica¬ 
tion  using  photorefractive  txa/o-wave  mixing 


4.1.  Image  amplification 

The  energy  transfer  in  photorefractive  two-wave  mixing  dis¬ 
cussed  in  Sec.  3.1  can  be  used  for  image  amplitication.' '  Spe- 
citically.  the  interaction  of  a  .strong  pump  beam  with  a  weak 
image-bearing  beam  inside  a  photorefractixe  medium  results  in 
an  amplification  of  the  image  information.  To  preserve  the  con¬ 
trast  of  a  gray  scale  transparency  ,  it  is  crucial  that  the  amplifi¬ 
cation  factor  be  unilonn  spatially  and  throughout  the  intensity 
dx  namic  ranee  of  the  imaee-carrMne  beam  This  issue  has  been 

*  <.  *.  <  T  “* 

analyzed  by  Fainman  and  Lee."  Their  results  indicate  that  a 
uniform  intensity  gain  of  more  than  40  dB  oxer  an  input  dynamic 
range  of  four  decades  li  e..  13  bus i  is  possible  in  photorefractixe 
crystals  such  as  BaTiO-  and  SBN  with  strong  coupling  I7L  ~ 
10,  where  y  is  the  coupling  coefficient  and  L  is  the  interaction 
length  I. 

The  basic  idea  described  above  is  illustrated  in  Fig.  4  along 
with  the  expenmental  results  Both  the  pump  and  the  image- 
carrying  beams  were  split  off  from  the  same  laser  output 
(514  5  nm).  The  signal  beam  was  first  attenuated  to  be  about 
SIX  orders  of  magnitude  weaker  than  the  pump  and  then  sent 
through  the  transparency  ti'  carry  the  image  information  pnor 
to  mixing  with  the  pump  beam  inside  die  crystal.  An  xyz-cut 
BaTiO:  crystal  was  used  as  the  ni'nlinear  medium.  The  close 
resemblance  between  the  am,'litied  image  (with  a  gam  of  ap¬ 
proximately  ZO.fKXIi  and  the  ongmal  image  indicates  that  the 
gain  IS  sufficiently  umliirm  oxer  the  dy  namic  range  of  a  typical 
gray  scale  picture  to  preserve  the  contrast. 

In  photorelraL live  crxstaN  such  as  BaTiO-,  and  SBN.  the 
phase  shift  between  the  inierlerence  tringes  and  the  index  grat¬ 
ing.  in  the  absence  I't  an  externally  applied  electric  held,  is 
approximatelx  -  2  (see  Se^  3  I  1  .As  a  result,  the  coupling  con¬ 
stant  and  hence  the  signal  gain  are  strimgest  when  the  optical 
Irequcncies  ot  the  pump  aiui  probe  beams  are  identical  |see,  foi 
example.  F-U  il4)  in  Kcl  54;  The  signal  gain  decreases 
(irasikaily '  '  as  the  optical  trcqiiencx  diliercncc  between  the 
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pump  and  probe  beams  exceeds  the  reciprocal  of  the  photore¬ 
fractive  response  time.  Since  the  typical  response  time  of  a  high 
gam  photorefractive  crystal  (e  g.,  BaTiO?  or  SBN)  is  of  the 
order  of  a  fraction  of  a  second  at  an  intensity  level  of  the  order 
of  1  W/cm‘,  the  frequency  difference  of  the  two  beams  has  to 
be  less  than  10  Hz.  tactically .  this  requires  that  the  two  beams 
originate  from  the  same  laser  To  overcome  this  limitation,  there 
have  been  numerous  investigations  worldwide  to  search  for  al¬ 
ternative  materials  or  processes  w  ith  much  faster  time  response. 

Several  other  image  processing  applications  employ  photo- 
refractive  two-wave  mixine  These  include  incoherent  to  co¬ 
herent  image  conversion.'  positive  and  negative  image  trans¬ 
fer.'"'^  dynamic  range  compression, NLO  range  imaging.'’' 
visualization  of  vibrational  mode  patterns.  '  and  novelty  fil¬ 
ters.'"*  In  Sec.  4  2  we  describe  a  unique  application  in  image 
subtraction  using  phase-conjugate  interferometry  . 

4.2.  Image  subtraction  and  “exclusive  OR"  logic 
operation  using  phase-conjugate  interferometry 

The  PCMI  examined  in  Sec.  3..^  can  be  readily  adapted  to  per¬ 
form  several  mathematical  operations  over  a  (or  a  pair  of)  two- 
dimensional  spatial  pattemi  s  i  of  amplitude  or  phase  modulation 
These  include  parallel  image  subtraction  and  addition,  negation 
and  XOR  logic  operations,  and  spatial  and  temporal  differen¬ 
tiations,  \\  e  bnefiy  review  some  of  these  applications  in  this 
section.  The  major  advantages  of  the  techniques  using  PC.Mls 
over  conventional  approaches  are  dynamic  stability  and  align¬ 
ment  insensitivity.  Details  are  discussed  in  Ref,  6(1. 

Optical  imam  suhtractum  Optical  image  subtraction*^'  has 
been  a  subject  of  considerable  interest  because  of  its  potential 
in  many  impt'rtanl  applications.  In  this  technique,  one  image  is 
subtracted  optically  from  another  to  detect  the  differences  be¬ 
tween  them 

The  PCMI  descnbed  earlier  can  be  used  as  a  real-time  image 
.subtractor.  Referring  to  Fig,  eibi.  consider  the  insertion  of  a 
pair  of  transparencies  (with  intensity  transmittance  functions 
T|(,x.y  I  and  T;(x.y  l)  in  the  beam  path  in  each  ann  of  the  inter¬ 
ferometer  equidistant  trom  beamsplitier  BS,  By  an  analysis  sim¬ 
ilar  to  that  in  the  previous  section,  it  can  be  shown  that  when 
the  two  phase-conjugated  beams  recombine  at  the  beamsplitter, 
the  image  intensity  at  output  plane  A  is  given  by 

Ui  \  y  1  ==  F:  ‘  p  ■  t'r  T;i  v.y  I  -  r'lTax.yi'  .  (22i 

where  p  is  the  reflection  coefticieni  of  the  phase  coniugators 

(assuming  equal  phase-coniugate  reflectivity  i  Using  Stokes 
equation  [Usj  (21)1.  this  tsccomc., 

Uix.yi  =  T, '  p 'RT  T|i  x.y  I  -  Tax.yi'  .  i2.''i 

where  R  r‘  and  T  =  t '  are  the  intensity  reflectance  and 

transmittance,  respectively,  tor  bcari'ispliiter  B.S.  .Note  that  the 
output  intensity  is  proportional  to  the  square  ol  the  dillercnce 
between  the  inten.ity  transmittance  functions 

Real-time  parallel  image  subtraction  using  a  PUSH  that  in¬ 
corporates  a  sell-pumped  BaTiU).  crvsia!  as  the  Pd'.M  lat  .“'14.5  nmi 
has  been  demonsiraicd  '  ""  .An  example  ot  the  exfserimenta! 
result'  is  given  in  Fig  5  The  horizontal  and  vertical  bars  are 
the  iitidccs  at  output  port  (see  f  ig  .Ubij  ol  traiisp.ircns le-  1 
and  2  respectively,  when  the  illuminating  be.im  for  trie  other 
arm  is  hlosked  The  ngh'  chcckerho.ir.;  p.ittern  represent'  the 
coherent  subtraction  ot  the  two  images  due  to  destructive  inter- 


Fig  5.  Experimental  results  for  image  subtraction  using  a  PCMI.  The 
horizontal  and  vertical  bars  are  the  images  of  transparencies  1  and 
2,  respectively,  when  the  illuminating  beam  for  the  other  arm  is 
blocked  The  checkerboard  pattern  is  the  intensity  distribution  of 
the  coherent  subtraction  of  the  two  images. 


ference  when  both  illuminating  beams  are  present.  .Note  that  the 
intensity  distribution  where  subtraction  takes  place  is  fairly  uni- 
lorm  and  is  very  close  to  that  c'f  the  true  dark  background  li.e.. 
dark  squares  where  the  dark  regions  of  the  bars  overlap) 

"Exclusive  OR"  loKic  operaium:  The  image  subtractor  can 
als(v  perform  logic  operation.  Consider  the  case  in  which  both 
transparencies  are  either  1  or  0.  According  to  Eq.  (23).  a  com¬ 
plete  cancellation  would  require  that  the  two  transparencies  be 
identical.  An  output  intensity  of  1  will  appear  at  port  .A  when 
only  one  of  the  transparencies  transmits  Thus,  such  an  image 
subtractor  can  act  as  an  XOR  gate.  Nk’hen  the  transparencies  are 
encoded  with  matrices  of  binary  data,  such  an  image  subtracior 
acts  as  a  two-dimensional  array  of  XOR  gates, 

Interisin  iruerswu  A  special  case  of  image  subtraction  is 
intensity  inversion,  which  is  obtained  by  renuming  one  of  the 
transparencies  such  that  the  transmittance  becomes  unity  in  one 
arm  E.xpenmental  results  are  shown  in  Fig.  6 

Amplitude  imai^c  subtraction:  The  image  subtraction  mcthixis 
descnbed  thus  far  result  in  subtraction  of  two  intensity  patterns 
due  to  double  passes  of  light  through  each  transparence .  A  new 
method  that  prov  ides  parallel  and  real-time  amplitude  subtraction 
of  two  images  by  using  holographic  interterence  in  photore¬ 
fractive  media  has  recently  been  demonstrated  *"  This  methcHl 
U'c's  a  double  Mach-Zehndcr  interferometer  w  ith  a  BaTiC);  PCM 
and  requires  only  a  single  pass  of  light  through  the  transpar¬ 
encies 

Spatial  diflrrentialion  A  PCMI  w  ith  a  pair  ot  identical  spatial 
patieni'  lone  in  each  arnu  can  be  used  to  perform  spatial  dif- 
tercntiation  Referring  to  Fig  3ibi.  if  one  ot  the  transparencies 
is  displaced  rel.ito- l  to  the  other  along  any  direction  x  perpen¬ 
dicular  to  the  optual  axis  by  an  amount  Ax.  Fxj  (2v)  becomes 

I.v'x.vi  -  [-.  ■  (I  R  r 'l,.i  \  ^  Ax, VI  -  'P.iv.yi  (24i 
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Fig.  6.  Experimental  results  for  the  intensity  inversion,  (a)  Intensity 
distribution  of  the  phase-conjugate  beam  in  the  first  arm  with  the 
transparency  removed  [h  ^Tilx.yl^  =  1],  (b)  image  of  the  trans¬ 
parency  in  the  second  arm  [I2  ^  T2lx,y)  ^],  and  (c)  the  intensity  in¬ 
version  of  (b)  [I  X  Ti(x,y)  -  T2(x,y)*  =  1  -  T2(x,y)^l 


where  To  is  the  iransparencs  intensity  transmittance.  This  equa¬ 
tion  can  be  interpreted  as  the  discrete  version  of  spatial  differ¬ 
entiation  with  respect  to  the  x  coordinate.  Note  that  such  an 
operation  leads  to  edge  enhancement  Using  a  slightly  different 
expenmental  configuration.  Kwong  et  al/"  demonstrated  both 
first-  and  second-order  spatial  differentiation  of  a  tv.o-dimen- 
sional  transmission  pattern 

Temporal  difierentiution:  .An  architecture  that  uses  a  two- 
dimensional  spatial  phase  modulator  in  one  arm  of  a  PCMl  lor 
the  temporal  differentiation  of  a  scene  has  been  proposed"  '  and 
demonstrated,'^ 

4.3.  Matrix  operations 

In  this  section,  we  review  the  use  of  four-wase  mixing  in  non¬ 
linear  media  to  perform  matrix-xector  and  matrix-matrix  mul¬ 
tiplications.  It  was  pointed  out  in  Sec.  3,3  that  four-wa\e  mixing 
IS  an  NLO  process  in  which  three  input  waves  mix  to  yield  a 
fourth  wave.  In  phase-matched  four-wave  mixing,  the  three  input 
waves  consist  of  two  counterpropagating  pump  waves  E;  and 
E:  and  an  arbitrary  probe  wave  E-.,  These  three  waves  are  cou¬ 
pled  through  the  third-order  susceptibility  \  "  of  the  medium 
A  fourth  wave  Ea  is  generated  that  can  be  written  as' 

Ea  ^  V'UiE.-E!  .  (251 

This  indicates  that  four-wave  mixing  can  be  used  for  the  mul¬ 
tiplication  of  Signals  In  aildition,  if  we  use  the  parallel  nature 
of  optisal  wave',  each  wave  can  carry  spatial  information  for 
the  purpose  of  image  processing.''  logic  operations.'”  ’'''  ''  nii- 
meru  pnKessmg.'^'  and  matrix  operations  ” '  ‘  In  what  follows, 

we  describe  some  unique  concept'  that  use  the  two  transverse 
dimensions  to  carry  matrix  infomiation  for  the  purpose  of  matrix 
multiplicatK'n 

Matrix  multiplication  K-tween  two  N  '  N  matrices  can  be 
written  as 

('  =  AB  iTfo 


w  here 

(  llA.b,  i:'i 


Eg  -  gHwt-k  ri 


Ibi 

Fig  7,  Schematic  diagram  illustrating  the  basic  idea  of  optical  ma¬ 
trix-matrix  multiplication  using  nonlinear  four-wave  mixing. 


spectively.  Beam  1  is  chosen  lo  propagate  along  the  y-axis  and 
beam  2  along  the  x-axis.  The  matrices  can  he  either  continuous 
or  discrete  In  the  discrete  case,  each  beam  consists  of  a  matrix 
of  beamlets.  as  shown  in  Fig  7(bi 

In  the  nonlinear  medium,  the  two  matrix-carry  ing  beams  form 
an  interference  pattern,  and  a  volume  grating  is  formed  This 
grating  contains  inlomiation  about  the  product  of  the  matnx 
elements  of  the  two  matrices  and  can  be  wntten  as 

An  =  n.'.Aix.zlB'iz.ylexptiKT  1  c.c  ,  (28i 

where  K  is  the  difference  ot  the  wave  vectors  of  the  matrix- 
carrying  beams  The  parameter  n;  is  the  Kerr  coefficient  and  is 
proportional  to  the  third-order  susceptibility  \' ''  of  the  medium 
Note  that  the  nonlinear  response  of  the  medium  performs  the 
function  of  parallel  multiplication 

The  volume  grating  is  read  out  by  a  third  beam,  which  can 
simpiv  be  a  plane  wave.  The  diffracted  beam  consist'  ol  the 
integrated  contribution  from  each  part  of  the  grating  along  the 
beam  path  and  can  be  written 


Ni'ic  th.it  mairix  multiplKatmii  consist'  ol  ivso  main  operations 
a  parallel  nuiltipli.  alion  and  .1  suinriiaiio: 

kelerrinj  to  Fig  ".ji,  eori'ider  .1  iou'  wave  mixing  config 
ur.iIioM  s.iitablc  tor  matnx  i;iiii:i;di.,aIioi,  H-.'ani'  I  and  2  contain 
the  intorni.ilion  about  the  'wo  n’.ilnecs  A  and  Bi/.vi.  re 


C'lx.yi’'  j  .A' x.ziB'iz.y  id.'  ,  (2Ui 

where  the  integraium  is  earned  out  along  the  beam  path  Note 
that  the  integration  completes  the  operation  cif  m.itnx  multipli¬ 
cation  The  inlormalion  about  the  product  ol  these  two  matrices 
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Fig  8  Schematic  diagram  illustrating  the  basic  idea  of  optical  matrix-vector  multiplication  using  nonlinear  four-wave  mixing 


is  nov.  impressed  on  the  tran^^e^se  spatial  distribution  of  the 
diffracted  beam.  In  ,  'actice.  the  dimensions  of  the  matnces  ma\ 
be  limited  b>  cross  talk  between  sarious  channels  due  to  dif¬ 
fraction  and  or  imperfection  in  the  optics, 

Owinc  to  the  phase-matchinc  requirement,  the  readout  beam 
must  be  incident  along  a  direction  that  satishes  the  Bragg  dif¬ 
fraction  condition  to  achieve  high  efticiencv ,  In  anisotropic  non¬ 
linear  media,  the  polarization  states,  as  well  as  the  direction  of 
propagation,  can  be  chosen  such  that  the  largest  of  the  nonlinear 
Susceptibilities  is  fulls  utilized 

Four-wave  mi.xing  can  be  either  degenerate  or  nondegenerate. 
In  the  degenerate  case,  all  of  the  beams  have  the  same  trequencv . 
In  the  nondegeneraie  case,  the  frequencies  of  the  beams  can  be 
slightlv  different.  This  mav  be  useful  for  the  purpose  of  sepa¬ 
rating  the  diflracted  beam  trvim  the  undiftracted  portion, 

.Another  scheme,  shown  m  Fig,  8.  is  suitable  for  matnx- 
vector  multiplication  Here,  as  an  example,  consider  a  discrete 
case  in  which  we  need  to  carrv  out  the  multiplication  of  an  N- 
element  vector  and  an  N  s  N  matri.x.  The  vector  is  fanned  out 
into  N  rows  of  identical  vectors.  These  N^N  beamlets  are 
directed  to  a  nonlinear  medium.  The  matrix,  which  also  contains 
N  r  N  beamlets,  is  also  directed  to  the  medium  in  such  a  wav 
that  each  beam  of  the  matrix  is  counterjiropagating  relative  to 
the  corresponding  beam  of  the  vectiT  Thus,  in  the  medium  there 
are  N  '  N  spaiiallv  separated  regions,  pumped  bv  a  pairot  coun- 
terpropjgating  beams  Now,  .S' v  N  probe  beams  are  directed 
mil'  the  medium  in  such  a  wav  that  each  probe  beam  can  prop- 
agaie  through  an  inlerscviion  region  The  probe  beams  will  be 
"['lane  wave  '  bcamlct'  propagating  in  parallel  As  a  result  ol 
the  tour-wave  miv  •■g.  each  probe  beam  generates  a  phase-con 
lu gated  beam  tha: .  ithin  a  proportional  it  v  factor,  can  be  w  ritten 
VlM.iiaoi  Bv  using  a  cvlindri^al  lens,  ue  obtain  a  summation 
('V  er  I  I  hi; ' .  w  e  h.i  v  e 

h-o  i.o.i  iTo 


where  atji  is  the  jr/i  element  of  the  vector  a  and  Mii.j)  is  the 
matrix  elemeni.  Such  a  scheme  for  inatnx-vector  multiplication 
can  also  be  used  for  matnx-matnx  multiplication  bv  decompos¬ 
ing  a  matrix  into  column  vectors  and  then  multiplying  tlie  matrix 
with  each  of  the  column  vecti>rs, 

•As  a  result  of  the  nature  of  the  ftvur-w  av  e  mixing  process  in 
the  medium,  this  matrix-vector  multiplier  operates  on  the  held 
amplitudes  and  thus  can  be  used  to  ’  andie  matrices  and  vectors 
with  complex  elements.  It  is  a  coherent  device  rather  than  an 
incoherent  one.  This  aspect  is  distincth  different  from  most  of 
the  earlier  approaches,  which  are  all  incoherent.  U'hcu  the  device 
IS  operated  in  the  coherent  mode,  the  phase  of  each  beamlct 
must  be  maintained  uniformlv  over  the  transverse  dimension  of 
the  beamiet.  In  addition,  such  phases  mu.'t  also  be  maintained 
fixed  in  the  summation  process.  If  these  phases  are  not  uniform 
over  the  beamlets,  the  final  step  becomes  an  incoherent  sum¬ 
mation  as  a  result  of  the  spatial  averaging.  Under  such  circum¬ 
stances.  this  matnx-vector  multiplier  operates  on  the  intensities 
and  thus  handles  onlv  positive  numbers 

Vve  now  describe  a  few  specihe  examples  of  optical  matrix- 
vector  multiplication  and  optica!  matnx-matnx  multiplication 
using  vanous  architectures  based  on  four-wave  mixing  and  phase 
conjugation  Other  approaches  to  optical  matrix  processors  are 
discussed  in  several  review  articles 


43.1 .  Optical  matrix-vector  multiplication  usinp 
photorefractive  four-wave  mixing 

In  the  hrst  experimental  demonstration  I't  optica!  iiiainx  vector 
multiplication  using  ti'ui-wave  mivng.'"'  ph.isc  coniugation  in 
a  phoiorefractive  BaFiO.  crvstal  w,.  iseJ  to  perform  pixel-bv- 
pixel  multiplication  The  sumnidtion  required  to  ohiain  matnx- 
vectoi  produc's  was  perfonned  subsequentiv  bv  a  cvlindrisal 
lens  external  to  the  crystal 

.A  schematic  diagram  of  the  experimental  contieuratioi,  and 
the  result  is  shown  in  I  ig  d  Mask  .M ,  imprints  the  matrix 
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into  the  two  counterpropagating  pumping  beams  A  collimated 
uniform  beam  is  used  as  the  probe.  The  phase-conjugate  output 
(31,  at  the  focal  plane  of  the  summing  cylindrical  lens  is  shown  in 
the  lower  right.  The  intensity  of  the  output  spatial  pattern  ts 
proportional  to  the  product  vector  v.'  given  by 


and  mask  M;  impnnts  the  vector  (repeated  on  the  mask  instead 
of  using  a  cylindncal  lens  to  fan  out  the  vector) 


(32) 


Fig  9  Experimental  configuration  and  the  result  of  optical  matrix- 
vector  multiplication  using  photorefractive  four-wave  mixing. 
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With  a  slight  mcxlification  of  the  experimental  geometry .  it 
has  been  successfully  demonstrated  '*'  that  both  the  pixel-b> -pixel 
product  and  the  summation  can  be  done  inside  the  nonlinear 
medium  to  achieve  the  desired  matnx-vector  product  without 
Use  of  an  external  cylindrical  lens. 

The  expenmental  geometry  is  illustrated  schematically  in 
Fig.  10  The  readout  vector,  which  consists  of  a  vertical  column 
of  beamlets  of  equal  intensity,  enters  the  crvsial  at  an  oblique 
angle  such  that  each  beamlel  traverses  all  of  the  counierpropa- 
gating  pumping  beamleis  in  the  same  hon/imtal  plane  inside  the 
crystal.  Since  the  proper  elements  of  the  matrix  and  vector  are 
encoded  in  the  two  counterpropagating  beams  v  la  the  appropriate 
masks,  the  phase-conjugate  output  of  each  heamlei  cimsist'.  of 
the  sum  of  the  prouuct  resulting  from  each  individual  encounter 

.An  expenmental  result  for  multiplication  of  a  4  >  4  binary 
matrix  and  a  4  x  1  binary  vector  is  shown  in  Fig  1 1 .  The  mag- 
nthed  images  of  the  masks  used  to  encode  the  matrix,  vector, 
and  probe  are  shown  in  Figv.  1  liai  through  1  lici.  respectivelv , 
The  image  representing  the  final  matrix-vector  product  is  shown 
in  Fie.  1  ltd  I 


4.J.2.  .\fatrix-matm  multiplication  by  means  of  color 
multiplexing 

1:  ha'  been  expenmentally  demonstrated  that  matrix-vector  mul¬ 
tiplication  Using  ['hotorefraciiv e  four-wave  mixing  can  be  ex¬ 
tended  to  matrix-matnx  multiplication  b'  color  multiplexing  ' 
The  basic  idea  is  to  decompose  the  problem  into  matrix-vector 
multiplications  and  Cuiry  out  all  oi  these  multiplications  simul¬ 
taneously  m  parallel  by  using  color  multiplexing  This  idea  of 
encoding  the  component  vectors  with  diflereni  colors  is  illus¬ 
trated  in  Fig-  12  for  the  case  ot  4  >  4  matrices  Each  row  vector 
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Fig  10  Modified  configuration  for  optical  matrix-vector  multiplicaiion  using  photorefractive  four-wave  mixing  In  this  approach,  both 
the  multiplication  and  summation  are  carried  out  inside  the  nonlinear  crystal 
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(a) 


(b) 


(c) 


(d) 


Fig.  11.  Experimental  results  of  optical  matrix-vector  multiplication  using  the  approach  illustrated  in  Fig  10. 


Fig  12  Schematic  di„gram  illustrating  the  basic  idea  of  optical  ma¬ 
trix-matrix  multiplication  using  color  multiplexing 


ot'  the  rrij'.n.v  Mi  illuminated  by  ime  eolor.  and  all  ol  the  color 
component^  are  then  comhined  h>  the  prism  taneular  multiplex- 
me  I  into  a  ^inele  rov.  prior  to  tunher  evpansU'.T  h>  anamorphic 
optics  (noi  shouni  to  match  ihe  mash  representing  the  second 
matrix  .After  proper  eleineni-b> -element  multiplications  and 
.summation  isumMing  optics  omitiedi  into  the  column,  the  re¬ 
sulting  multicolor  output  is  demultiplexed  into  diUcreni  color 
components  that  together  represent  the  tinal  product. 

4.3.3.  .Matrix-matrix  multiplication  usin^  spatial 
convolution 

In  addition  to  the  ahose  scheme,  another  method  tor  matrix 
matrix  multiplication  that  uses  spatial  ctinwilution  hy  means  ot 
fi'ur-wase  mixire  has  been  successtulK  demonstrated  '  The 
ke\  ditterence  between  this  approach  and  those  described  aK'se 
IS  that  the  nonlinear  crystal  is  kK'ated  at  the  common  bouner 
plane,  rather  than  at  the  common  image  plane,  of  the  input  matrix 
masks  The  encoding  scheme  is  also  dillerent  from  those  used 
m  the  earlier  approaches  Matnx-matrix  multiplication  in  lull 
parallelism  is  achie'ed  b\  space-multiplexing  by  mc^ns  I'l  spa¬ 
tial  convsdution  usine  degenerate  tour  wase  mixing 

Vi'hite  and  "laris’  ha\e  demonstrated  that  (spatial)  cor.vo- 
lution  and  sonelation  ol  two  (two-dimensional  i  pctiems  can  be 
ULhiescJ  in  real  ome  bs  touc-wase  mixing  in  the  common  bou 
ncr  plane  ('f  the  input  oattenis  and  rcsxirding  the  ph.ise-coniugated 
ou’i'Ul  at  the  correspt'ndmc  I'biect  plans  .A  typical  experimental 
I  ( "'ll is  s.)|'>wn  s^hernatisallx  in  big  I  s  I)  I  .'  is  a  small 
■  iK'.tur.  sinulating  a  point  soursc  sielta  tunctioni,  a  ratten: 
represent  me  'se  con  x  olun  .p  ol  ! '  1  and  I'd  ( I '  1  *  I  2 1  i  ob'-Txed 
IP  the  output  plane  The  design  ot  11  and  12  to  represent  ’wo 
m.itn.e-  si.  li..,:  tiieir  prodiim  o  repre'crted  by  (  l'i'2  is  noxx 
expl.tined 


Fig  13.  Schematic  diagram  illustrating  the  basic  idea  of  four-wave 
mixing  in  a  nonlinear  medium  located  at  the  common  -ourier  plane 
of  the  input  spatial  pattern. 


For  simp/icity .  consid-er  the  simplest  case  of  two  matnees  .4 
and  B  with  dimensions  2  r.  2  and  the  product  C"  =  .AH  gixen 
below : 


I  .^4 1 


[b:;  b, 


I  b;  b; 


(.TF| 


a.; hi-  -s  a:;b;;  aoh  ,  *  ai.-h.-; 

a;,h|i  *  a”b:i  a;.b  ;  ^  a.-T,-.- 


(,^hi 


Transparencx  L' 1  e  uisist-  of  lour  small  upenures.  each  with  an 
iniensiix  transmittance  pri'ponionti!  to  each  matrix  element  a, 
|see  big  I  bad  similar  transparency  I'd  conesponding  to 
the  tronsposf  ol  B  (i  e..  ro’xs  and  column  ,  interchanged '  is 
shiiwn  ill  big  14ih)  Noie  that  the  xcnical  scparalioii  in  I  '  is 
considerablx  larger  than  that  m  I'd.  x'.hhc  the  hon/onlal  elen.cnt 
separations  are  ideniical  in  11  and  Id  The  two-dimensional 
spatial  conxxilulion  ol  L' 1  and  Ud  is  shi'xxn  in  big  |4ui.  It 
consists  (it  a  total  x'f  U'  spi'ts  with  the  intetisit;.  ol  each  pr.'- 
poriional  to  the  cross  product  ol  all  ot  the  eleii'enis  in  .A  and  B 
The  tour  doubi-,-|s  Ixinc  along  me  y  axm  in  big  Nix  i  are  drawn 
siightlx  otisci  troin  ihc  axis  lo  -.-xpsise  ihc  mdixidual  xomponeni 
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C  =  AB 

^11  =  “Ill’ll  '  "12‘>21 

^12  =  «11‘>12  *  ®12t>22 

C21  =  »21*’n  "  ‘Z2^2t 

=  11  ®';2\  ^  *>12  j  C22  =  ®21*>12  *  »22*>22 

a2i  822  \  *>21  *>22  ' 


V 


a  (b,  |c; 


Fig  14  Schematic  diagram  illustrating  optical  matrix-matrix  multiplication  by  convolution  for  2  x  2  matrices  (a)  Matrix  A  and  mask  U1  to 
encode  A,  (bl  matrix  B  and  mask  U2  to  encode  B^,  and  (c)  spatial  pattern  resulting  from  convolution  of  U1  and  U2 


la'  (bi  Ic!  Idi 


2 


Fig  15  Experimental  results  for  optical  matrix-matrix  multiplication  by  convolution  by  means  of  four-wave  mixing  in  the  spatial  frequency 
domain  lai  anC  (bi  are  the  matrix  masks  used  for  encoding  the  matrices  A  and  B’,  (c)  shows  the  result  of  their  spatial  convolution,  and  Id)  is 
an  overexposed  version  of  (c'  to  reveal  the  cross  terms  that  do  not  contribute  to  the  matrix  product 


In  prji-ti^'c.  itvj'-i.'  tui'  conip'incrit'  spatialU  uvuiijp  on  the  \- 
j.\i^  Tha  microiiic'  n:  the  lour  do-uhictx.  iroiii  top  ii'  ho!n>ni 
on  thu  \  a\i'.  a!-,'  propomona;  to  ^  and  rC'-pc-.- 
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■Another  scheme  for  optical  niairix-xecior  (and  matrix-matrix  ) 
multiplications  that  uses  a  phase  conjueator  iw  ith  a  finite  storage- 
time  1  in  conjunction  \xith  a  spatial  light  modulator  iSLMi  to 
eliminate  the  pixel-bv-pixel  alignment  requirement  at  the  cost 
of  some  icdiiLtion  in  paralh  .ism  iN"  2  instead  of  N'l  has  been 
reponed  ''  Phase  abcrratioii  due  to  imperfection  in  the  optics  is 
self-conected  b\  the  phase-conjugation  process  The  optical  s\s- 
tem  involved  is  relati.elx  simple  civ',  ed  \xith  the  other  ap¬ 
pro-aches  Tigure  lb  IS  a  schematic  diagram  illustrating  the  basi^ 
concept  .An  SLM  is  used  li'  impress  the  matnx  and  vector 
information  m  sequence  tm  an  input  laser  beam  This  beam  is 
directed  toward  a  self-pumix-d  F>C'M,  which  has  a  hmie  time 
res[)onse  le  g,.  BaTiCt-i  .A  cvlmdncal  lens  is  used  in  the  phasc- 
eoniugate  output  beam  path  to  perlonii  the  summation 

The  principle  of  operation  is  as  follows  The  SLNf  first  im¬ 
presses  die  riutm  inIormaii.>n  onto  the  inpui  laser  beam  This 
beam  is  then  mvident  into  a  self  pum[X-d  f’CM.  which  stores  the 
matiix  mfonn.ilion  after  a  finite  grating  formation  time  When 
tile  matrix  mtormaiior,  i-  rem-wed  ftom  the  STM  and  all  of  the 
p-.xel-  are  lunied  ml-'  m.iximum  transmissi-tn  sondition.  the  phase 
."■ouc.i’e  fx-.i;i:  tf.a!  -.''nl.-.m  the  res  - -n 'tr.'ctioti  of  the  matrix 
ml- 'nii.o ;  -r  exisi-  l-n  ,,  r.inie  d.iitaii- T, ,  wt.iji  depends  on  the 
stiei.e:'  ■-!  t’.e  irijHi’  Me, of  t'-e.i::;  [i.i:;-.e  tho  tuiK  .  1!  tfie  next 
fr  ;■  ■  ■  '  me  Yl  ‘sf  ■  ir:  le-  ifu  v  e-  '>  -  t  1 1: !  -I  mat p.im!  lei  1  uji  - 
,.1-1  I-  [KM’.  Ii.e  iMi.-.se  ..  'oee.;'.  iX'ai:i  piopjg. 

■  .  -'.r  .  ,e'  M  M  Me-.  ,  tie  S  ,  t  :  -e;  te-etMeC  .<s  a  twe 
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Fig  16  Schematic  diagram  illustrating  the  concept  of  optical  ma¬ 
trix-vector  multiplication  using  a  photoref.  active  phase  conjiigator 
in  conj'jncticr-.  -^ith  an  SLM 

dimension  of  the  vector.  A  cvlindrical  lens  in  the  output  port 
perlomis  the  summation.  The  dark  siorae-e  time  during  which 
the  matnx  information  can  be  retnev-^d  is  determined  by  the 
photorelractive  material  and  the  pumping  conhguration.  It  ranges 
trom  secE'nd^  ti.'  mieroseconJs,  In  this  architecture,  the  same 
input  beam  is  used  tor  alternately  writing  and  reading  the  holo¬ 
gram. 

This  system  can  al-'O  perfiirm  matrix-matrix  multiplication  by 
time  multiplexing  In  this  case,  each  ci'lumn  vector  \\  (i  =  1 
to  N  '  constituting  the  second  matrix  .M;  is  sequentially  impressed 
onto  the  beam  ts.  r.iultiply  with  the  first  matrix  M]  according  to 
ihe  matnx-vostor  nuiltiplication  scheme  described  above  lii 
avoid  degradatk'H  o!  the  intormation  of  ,\1;  stored  in  the  phot(>re- 
frastive  hologram  during  the  readout,  it  o  necessary  to  refresh 
the  holographic  memory  with  .ST  to  restore  its  diftraction  ef¬ 
ficiency,  This  can  be  d''ne  by  rcimpressing  Ml  onto  the  beam 
after  each  readout  cycle  Consequentiy ,  a  tola'  of  2N  clo^k 
cycles,  consisting  o*  N  cycle'  ol  wntc  and  .\  cycles  tv|  read, 
will  be  required  to  carry  out  trie  nuiliinlk.iiion  of  two  .N  ■  ,\ 
matrices 

4.4.  Optical  interconnection 

(dpiktil  interconnection  linking  laser  arrays  and  detcctiT  arrays 
plays  a  key  role  in  optical  coniputirig  '  Conceptually  .  such 
an  ir;’,'',.>innecnon  can  he  adiieved  hy  using  iiptical  matrix- 

ve^ti'i  muliipli.  atioi,. 

V  M  V  I "  I 

where  v  is  the  input  veelor  representing  the  signals  earned  by 
an  arrav  o!  .N  lasers  and  v  ’  is  the  ouipul  vector  representing  the 
sienals  received  bv  the  array  ol  \  detectors  M  is  an  N  >  N 
binarv  matrix  represeniing  the  interconnection  patterr 

Vk  hen  an  SL.M  o  used  a'  the  matrix  mask,  the  .nterconneetion 
IS  re . I Tit.-ciirahle  Suppose  the  mtcrconnecTion  mask  consists  ol 
an  N  r  S  arrav  >  1  tran-iiiissioii  v  .ndows  The  input  vector  la 
row  id  N  elements  I  is  tanned  out  using  astigimiiie  optics  s(i  that 
the  lichl  trom  each  elemeni  is  hroadcasi  (ivcr  a  correspondme 
coluiii';  wind'.'W,  ,-Mier  traiismissiMn  through  the  iniercot! 
nc'.'.ioii  in.i-i  .  simil.i:  opiics  ,.ire  Use,l  to  eollcvl  the  light  Irom 
cM.i  r-w  ol  trie  windo'w.  .md  to  suin  ttu-  ouipul  into  a  column 
'o.-.'ia  o!  \  cleiiiei:;,  Sudi  all  ar.hi'.cwlure  prinules  Ihe  ncc 
e's.mc  \  p.ir.iilcli'm  tiowece''  a  laree  traction  ol  en-ergv  o 
jb-orn-ed  t' .  I’r.  ;  ■  .c.cp,  ..nri-  o;  Sl.M  \\  huip  used  as  a  cn'"b,i: 
"■cii,,',  -ij  t  ,11  ,1,' 0  i'.i:,  ,  ’  jr,^  [■>,,,  ,,p,  eT.ercv  e‘I)^ie.n,.v  ol  or:'. 

I  S  wtv:  N  !  ■  It,'.' d::!,.'o  'Mh,,rra‘.  '  I  in- o^^urs  fv-„.iiis,.- 
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vector  does  not  pass  through  the  crossbar  mask  The  energy  loss 
increases  as  the  dimension  of  the  array  increases  and  is  reterred 
to  as  the  "fan-out  energy  loss."  For  a  1000  x  1000  crossbar 
switch,  the  loss  due  to  fan-out  can  be  as  high  as  99.991.  This 
IS  not  acceptable  m  high  speed  computing  because  signals  are 
passing  through  Ihe  SLM  at  gigabit  rates  and  the  energy  loss 
can  be  enormous  In  addition  to  the  inherent  fan-out  energy  loss, 
all  SLMs  have  finite  insertion  loss  due  to  imperfect  transmissions 
such  as  those  that  occur  from  absorption  and  scattering.  Includ¬ 
ing  such  insertion  losses,  the  energy  efficiency  to:  a  crossbar 
switch  would  become  lN.  where  t  is  the  transmittance  (t  <  1) 
of  each  of  the  w  indow  s 

Recently,  a  scheme  has  been  developed^'' to  form  a  dy¬ 
namic  hologram  inside  a  photorelractive  crystal  such  that  the 
hologram  wiJ)  be  capable  of  performing  reconfigurable  optical 
interconnection.  Such  a  scheme  utilizes  nonreciprocal  energy 
transfer  in  two-wave  mixing  to  achieve  an  extremely  high  energy 
efliciency 

Relemng  to  Fig  17.  we  describe  this  new  meihcxj  of  recon- 
figurahle  optical  interconnection.  Figure  17iai  d^.MTibes  a  one- 
dimensional  case  to  explain  the  concept  .A  small  traction  of  the 
incideni  beam  is  split  off  using  beamsphtier  BS  Tne  reflected 
beam  (probe  beam  I  is  then  expanded  by  using  a  cylindrical  lens 
before  entering  the  SLM,  In  this  example,  the  input  beam  is  to 
be  connected  to  detectors  b  and  d  as  prescribed  by  the  SL.M. 
The  transmitted  probe  beam  is  then  recombined  with  the  main 
beam  inside  a  phoiiiretraciive  crystal.  As  a  result  of  ntvnrecip- 
rocal  energy  coupling,  almost  all  of  the  energy  in  the  main  beam 
IS  translerred  to  the  probe  beam,  which  carr.cs  the  interconnec¬ 
tion  pattern  The  result  is  an  optical  interconnection  with  a  very 
high  energy  efliciency 

Figure  17ib)  shows  the  reci'nligur.ible  inlerconnection  for 
laser  arrays  and  detector  arrays.  In  this  example  i4  >  4  mter- 
connectioni.  laser  1  is  to  he  connected  with  deiectv'rs  b  and  c. 
laser  2  with  detectors  a  and  d.  laser  .'  with  detectors  c  and  d. 
and  laser  4  with  detectors  a  and  c  In  terms  ol  matnx-vecii'r 
multiplication,  such  an  inter-.onneciK'n  can  be  written  us 
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where  v  i,  v;.  \  ..  and  v^  are  signals  carried  hv  ihe  laser  elements 
ol  the  array  \ 

,A  cylindrical  lens  is  used  to  Ukus  the  twii-dimensional  array 
of  beams  into  a  vector  Thus,  as  a  result,  detector  a  receives 
signals  trom  lasers  2  and  4.  detecior  b  receives  signals  trom 
laser  i,  detecior  c  receives  signal'  IrEim  lasers  I.  .T  and  4.  and 
deiectEU  d  receive'  signal'  trom  lasers  2  and  .v  Ihis  concept 
can  be  extended  to  uucrsonnect  N  lasers  with  M  detectors,  where 
.\  and  M  are  two  large  numbers 

The  two-wave  mixint:  descnbesl  ahove  ma\  be  viewed  as  a 
re.il-time  holography  in  wh.  h  the  reLortiing  and  readoai  iKcur 
simulianeouslv  inside  ihe  phoioretr.i.nve  crvsi.tl  The  beam¬ 
splitter  and  the  .SL.M  are  used  to  generale  the  retereiise  and  obiecl 
beams  to  record  a  volume  holo'cr.uii  lha;  represents  ihe  mtei 
coiinesiion  paiiern  i'  prescribed  bv  Ihe  SLM  Ihe  energv  eou 
pling  involved  in  iwo w.oe  niixing  ensu'C'  ih.i:  luc  liUra^lion 
etticiency  diirine  the  re.idoui  is  .ilni-"'  |f«'  this  requires  .i 
proper  orient, I'l.iii  oI  ihe  srvsl.;!  'o  ilia:  Ihe  eneicv  ol  the  le.nlou' 
he.iin  1'- gre.illv  depleted  llic  hicli  ere  c  v  ein.i-'iuv  lesull' Iioiii 
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Fig  17  Gchematic  drawings  of  optical  interconnections  using  dynamic  photorefractive  holograms  with  N  =  4.  (a)  1  x  N  optical  interconnection 
and  (b)  N  X  N  optical  interconnection. 


the  fa't  that  most  of  the  eneres  is  carried  h\  the  readout  beam. 
\Ahich  J(H.‘s  not  pass  through  the  SLM  but  is  diffracted  into  the 
interconnection  pattern  b\  the  hologram 

The  energ\  efticicncs  of  such  an  interconnection  pattern  can 
be  estimated  for  a  crossbar  suitch  as  follous.  Let  R  be  the 
reflectance  of  the  beamsplitter.  It  is  legitimate  t(i  assume  that 
the  beamsplitter  is  practicalK  lossless,  and  ue  may  assume  that 
the  surlace  of  the  photorelractise  crystal  is  antiretlection-coated 
so  that  the  hresnel  reflectnni  loss  can  be  neglected.  L  nder  these 
conditu'ns,  the  tuo  beams  that  arn\e  at  the  photorefractise  crys¬ 
tal  have  energies  (I  -  Ri  and  RlN.  respectively.  Inside  the 
crystal,  these  two  beams  undergt'  photorefractise  coupling.  As 
a  result,  most  of  the  energy  oi  the  pump  beam  i|  -  Ri  is 
transferred  to  the  probe  beam  RlN,  which  contains  the  inter- 
CvT.nection  pattern.  The  energy  et’ticiency  can  be  easily  derived 
and  IS  given  by'"' 
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Fig  18  Schematic  diagram  illustrating  the  relative  orientation  of  the 
pump  beam,  probe  beam,  and  c-axis  of  a  30' -cut  BaTiOj  crystal  for 
the  measurement  of  energy  efficiency  The  depleted  (transmitted) 
pump  beam  is  omitted  for  clarity. 
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where  m  is  the  beam  intcn'itv  ratio. 
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L  IS  the  interaction  length,  y  is  the  coupling  cotisiaiii.  and  a  is 
the  bulk  absi  irption  ci>cl  iicient  f-or  photoretra.live  cry  si.ds  such 
as  Bdlif),  and  SBN.  the  coupling  constant  |s  very  large  ti  e,. 
-yl  ■  •  li  The  cfliciency  can  be  writien  approximately  as 


.Notice  that  tor  larce  N.  the  energv  efticiency  is  limited  bv  the 
crystal  bulk  absorption  exp'  ol.i  and  is  maximi/ed  bv  using  a 
beamsplitter  with  a  very  small  reflectance  R  ti  c  .  R  (h 
In  evamining  die  energv  ellici'  icy.  wc  collimate  an  aruon 
ii'n  laser  beam  at  ^14  >  np,  mti.  2  nim  diameter  bv  ijse  ;  a  lens 
nt  local  lenctl)  t  2  III  -k  beam'plitler  is  used  to  redirect  y' - 
ol  the  cTiercc  throuch  the  Sl,\1.  est.iblisti,  the  probe  bc-arr. 
The  remaitider  ‘d  lt,e  energv  ir.insiniis  itirough  the  beamspliiier 
.ind  c '  s, Cu'e,  |p_,  pijt.iT.  Ir,  ou:  pielmiin.irv  experimeti'. 

Ih  ,'  Si  M  w  .1-  repia.  eil  w  nti  a  neutral  ilerisit .  I;  Her  w  itli  a  variah'e 
opto.,;  'ten-::,  i-  -imuLi'c  ih.  Ian,'i,t  enerev  loss  |  tie  Ivo. 
t'e  ■  i ■  t-,-:  -  ■  e - .  ;  |  a  ,  ,  r  - -la;  w  ill.  an  iuleracliiin  lencli. 


ol  approximately  5  mm  (see  Fig.  IS).  The  intensity  ol  the  pump 
beam  is  approximately  .k,2  W  cm",  whereas  the  intensity  of  the 
probe  beam  varies  from  U.16  V,  cm"  to  (klb  m\k  cm".  The 
amplitied  output  intensity  is  monitored  and  the  result  is  shown 
IP.  Fig  19.  where  the  energy  efticiency  (defined  as  the  ratio  of 
the  amplified  output  intensity  to  ihe  total  input  intensity  i  is  plot¬ 
ted  as  a  tunciion  of  1  .\  (where  N  is  the  number  ot  lan-out 
ehannels)  The  solid  line  represents  the  theoretical  ni  using  Exj 
(,'9 1  with  aL  and  yL  as  the  titling  parameters  The  results 
(ol.  =  9.2  and  nL  1 . 1 1  agree  with  those  obtained  from  ('ther 
independent  measurements  U'  within  IriN  Note  that  the  energv 
efficiency  is  better  than  24N  for  N  -  I  through  UKI.  In  other 
words,  if  the  tan-out  loss  js  9d''k  (e  g,,  a  UHi  >  !()(■  crosshari. 
the  energy  efhciency  of  th[s  new  scheme  tan  be  24  Kmes  betier 
than  that  of  the  direct  approach  Bulk  absorption  in  this  particular 
crvstal  accounts  for  a|iprc'ximately  65'",  of  the  energy  loss 
Fven  though  the  reconliguration  time  is  limited  by  the  phiUo- 
rcfractive  response  time,  which  is  typically  I'n  the  order  I'f  milli¬ 
seconds  at  modest  iniensuies,  the  phoioretraciive  interci'nnec- 
tion  sysiem  c;ui  aceepi  very  high  data  rate  signals  do  demonsiraie 
this  tact,  wc  frequency  modulate  a  laser  beam  largun, 
A  ■  .^I4,y  nmi  using  an  ac oiisto-iiptic  device  to  simulate  a  sig 

nal  that  u  to  he  inierconnccTcd  with  some  uulpui  The  signal 
Used  had  a  pulse  train  ol  IrequericV  -  (i  S."'.'  MH/  witti  eacli 
pulse  beine  approximalelv  I)  2  gis  w  ide  1  his  rate  is  dearlv  much 
hicher  than  the  rcciprowii  '!  the  ptiotoretracliv  e  resjionse  time 
I  lie  iiii'diilaled  l.tser  bc-am  was  then  spin  mto  two  beams  and 
mixed  III  the  erv  stai  a-  dcs^nived  belore.  and  Ilic  .implilied  probe 
hea;ii  was  mi.niioreC  vcilli  .i  ph' uodc  ice  lo:  I  t;c  oc^  illos.ejx' 
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Fig  19.  The  energy  efficiency  as  a  function  of  1/N,  where  N  is  the 
equivalent  number  of  fan  out  channels.  Open  circles  represent  the 
experimental  results,  and  the  solid  line  is  a  theoretical  fit  [using 
Eq  (39)1  with  aL  and  yL  as  the  fitting  parameters. 


Fig  20.  (a)  Signal  carried  by  the  probe  beam  and  (b)  amplified  signal 
carried  by  the  probe  beam  after  undergoing  two-wave  mixing 


trace  in  Fig,  20(.i)  shows  the  input  probe  signal,  and  the  trace 
in  Fig  20(b)  show,  the  amplitied  probe  signal.  The  results  show 
steads -state  resprrnse  in  which  the  temporalis  msKJulated  pump 
and  probe  beams  interact  simpls  b_s  diffracting  off  the  siationars 
index  grating  that  is  created  in  the  crsstal  alter  the  beams  interact 
for  a  tune  on  the  order  of  the  photorefractisc  response  time 
This  experiment  was  perlormed  merels  to  demonstrate  the  high 
signal  bandssidth  ol  the  ssst-em  (1ptimi/at;on  of  the  parameters 
was  not  done,  so  the  result'  shossn  in  Fig  20  can  ciearls  be 
impri  IS  ed 

5.  (  ()N(  I.l  SION 

Ue  have  de-'Dlx'd  snme  ol  the  inieresiine  and  important  pile 
non,'.' •  1.1  in  p:,.  ■;  Tctra..  i r. c  rued ..i  I’rioiorctr.i^tise  materials  sui.ii 
a-  ti.il  |(  1  StiN  and  HSf  I  ar-c  bs  tar  the  mos:  elticient  medi.i 


for  the  generation  of  holograms  and  phase-conjugate  waves  using 
relatively  low  optical  intensities.  In  addition,  two-wave  mixing 
in  these  media  exhibits  nonreciprocal  energy  transfer.  As  a  result 
of  the  strong  energs  coupling,  several  unique  photorefractive 
phenomena  occur  in  these  media.  These  include  self-pumped 
phase  conjugation.  nonreciprcKal  energy  transfer  without  phase 
cross  talk,  phase-conjugate  interferometry,  etc  These  unique 
phenomena  play  an  important  role  in  many  of  the  applications 
in  optical  computing  including  matrix-vector  multiplication,  par¬ 
allel  subtraction,  reconfigurable  interconnections,  etc.  Some  of 
these  applications  were  presented  and  discussed. 
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Using  Photorefractive  Holograms 
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ABSTRACT 

We  propose  and  demonstrate  the  use  of  Fourier  transform  to 
achicNC  maximum  energy  efficiencx'  in  a  pholorefi  acti\ e 
optical  interconnection.  The  results  of  experimental 
inxestigations  on  reconfigurable  optical  interconnections 
using  photorefractixe  holograms  in  a  barium  titanate  crystal 
are  presented  and  discussed.  High  energy  efficiency  is 
achiexed  bx'  matched  amplification  at  the  Fourier  plane. 
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1 .  Introduction 


Optical  interconnection  between  VLSI  circuits, 
computing  chip>  or  bo;^rd^  pla>s  an  important  role  in  optical 
computing  [1.2J.  Such  a  scheme  of  computing  proxides  the 
potential  of  achieving  extremely  high  speed  because  it 
utilizes  the  fast  sxxitching  of  electronics  and  the  wide 
bandwidth  of  optics  for  communication  [3].  A  generalized 
optical  crossbar  svxitch  can  provide  a  reconfigurable 
arbitrarx  interconnection  (4,5 1  betxxecn  an  arra>  of  N  lasers 
and  an  arrax  of  N  detectors.  Such  a  crossbar  can  be 
implemented  using  parallel  matrix-vector  multiplication 


which  provides  parallelism  15-8J. 


It  is  known  that  the  N-  parallelism  is  accompanied  with 
an  intrinsic  fan-out  energy  loss  [4,9],  The  fan-out  leads  to  an 
energy  efficiency  of  1/N.  We  recently  proposed  and 
demonstrated  a  novel  concept  of  reconfigurable  optical 
interconnection  19-11]  which  can  provide  high  energ\- 
efficiency.  Using  the  beam  coupling  in  photorefracti\ e 
crystals  [12,15],  such  a  scheme  of  interconnection  is  capable 
of  minimizing  the  fan-out  energy  loss  and  thus  achie\es 
extremely  high  energy  efficiency.  Although  the  basic  concept 
has  been  salidated  and  some  of  the  preliminarc'  results  ha\c 
been  reported  [9.11].  tlie  issue  of  energ\'  efficiencc’  remain^  a 
subject  of  experimental  irnesligation.  This  paper  reportv  the 
results  of  an  extensive  in\ estigation  on  the  energy  efficienc> 
of  reconfigurable  optical  interconnection  using 
phot  ore  fractise  holograms. 

In  what  follows,  we  will  briefl\  resiev.'  the  basic 
principle  of  operation  of  the  photorefractive  interconnection. 
We  will  point  out  the  use  of  Fourier  transform  to  achie\e 
maximum  beam  overlap  in  the  photorefractive  medium.  We 
then  discuss  the  results  of  our  experimental  in\ estigations. 


2.  Principle  of  Operation 


Figure  1  shows  the  basic  idea  of  reconfigurable 

interconnection  using  photorefractive  holograms.  A  beam 
splitter  is  used  to  split  off  a  small  portion  of  the  beam  which 
consists  of  beamlets  each  carrying  information  from  one  of 
the  sources.  This  small  portion  of  the  beam  will  be  called  the 
signal  beam.  The  beam  splitter  allous  the  majorits'  of  the 
energy  to  p^ss  through.  This  major  portion  of  the  beam  i^ 
called  the  pump  beam.  d'he  signal  beam  is  fanned  by  the 
cylindrical  lens  and  then  passes  through  the  spatial  light 
modulator( SLM )  which  contains  the  interconnect ioii  pattern, 
.After  passing  through  the  SL,\1,  the  signal  beam  is  imprinted- 
spatiall)'  uith  the  interconnection  pattern.  The  pump  beam 
does  not  pass  through  the  .spatial  light  modulator  and  thu^ 

suffers  no  energy  loss  due  to  the  fan-out  and  all  the  other 

loss  mechaniNms  associated  vsiih  the  SL.M. 

The  signal  bettm  which  contains  the  interconnection 
information  is  then  recombined  with  the  pump  beam  in  a 
photorefractive  crsstal.  Under  the  appropriate  conditions,  the 
signal  beam  can  be  amplified  at  the  expense  of  the  pump 
beam.  .Most  of  the  energy  of  the  pump  beam  will  be 
transferred  to  the  signal  beam  presided  the  length  cd 


interaction  is  large  enough  [14,15].  It  is  known  that 
photorefracti  ve  tv,o-wave  mixing  exhibits  energy  transfer 
without  phase  crosstalk  [16].  Thus  the  amplified  signal  beam 
will  also  bear  the  interconnection  pattern.  The  net  result  is 
an  optical  interconnection  with  high  energy  efficiency. 

Although  photorefractive  two-wave  mixing  is  in\ol\'ed, 
such  a  scheme  can  also  be  considered  as  a  real-time 
holographic  interconnection  with  a  very  high  diffraction 
elficiency.  The  combination  of  the  beant  splitter  and  the  SLM 
in  conjunction  with  the  pump  beam  is  used  to  record  a 
hologram  inside  the  photorefractive  crystal.  Such  a  hologram 
contains  the  interconnection  as  prescribed  b\’  the  SL.M.  Once 
the  hologram  is  formed,  the  pump  beam  can  be  diffracted  oft' 
the  hologr'mi  and  he  redirected  into  the  arrax  of  detectors. 
1  he  adxantage  here  is  that  the  recording  and  the  readout 
occur  simultaneously.  Th.  is  offers  the  possibilitx  of 
reconfigurab'e  interconnection  by  'ontrolling  the  SLM. 

do  achieve  maximum  energy  efficiencx,  it  is  im[iort.;:it 
that  the  signal  beam  and  the  pump  beam  oxerlap  completelx. 
Speci ficallx ,  beamlet  1  of  the  pump  beam  in  Fig.  1  must 
overlap  completelx  v\ith  column  1  of  the  signal  beam.  And 
siiiiihirlx .  beamlets  2,  .d,  and  4  must  also  overlap  ccmipletelx 
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uith  their  corresponding  columns  in  the  signal  beam.  It  is 
important  to  note  that  although  the  beamlets  in  the  signal 
beam  and  the  pump  beam  are  intrinsically  different  because 
of  the  interconnection  pattern,  the  individual  pixels  can  have 
identical  shape  (e.g.,  square).  Complete  overlap  is  possible  in 
the  Fourier  domain,  provided  all  the  pixels  are  identical.  Let 
s(x,y)  be  the  aperture  of  an  individual  pixel,  and  o(  u.\  )be 
the  Fourier  transform.  Fourier  transformation  of  a  beamlet  of 
the  pump  beam  and  the  corresponding  column  of  beamlets  of 
the  signal  beam  can  be  performed  by  a  lens  (see  Fig. 2).  The 
resulting  amplitude  distributions  at  the  rear  focal  plane  are 
g  i  \  e  n  b ; 


Pump:  P( u ,\  )  =  at  u ,\  )  — . —  ( 1 ) 

Signal:  Si  u.\  )  =  Z  Oi  u.v  )  exp  |  i  (u.v)  )  —  (2j 

n 


respect! vel>,  vshere  u,\  are  coordinates  in  the  Fourier  plane, 
'fhe  summation  in  equation  (2)  is  over  all  the  uindous  opeti 

in  the  column,  and  o^_(u,\  )  is  a  phase  \x'hich  depends  <n  the 


position  ot  the  windou. 


We  note  that  the  Fourier  transform  of  a  column  of  the 


signal  beam  consists  of  a  linear  superposition  of  identical 
patterns  each  with  a  different  phase  factor.  Such  phase 

factors  exp(io,i)  are  due  to  the  relative  positions  of  the 

windows  in  each  column.  Although  these  phase  factors  ma\ 
lead  to  interference  structures,  all  the  energies  are  confiticd 

under  the  same  envelop  a(u.\  ).  As  a  result  of  the  shitt 
in\ariance  it',  bouricr  transformation,  each  of  tiie  signttl 
beamlets  oserlaps  completel>  with  the  pump  beam  ;»t  the 
r-OLirier  plane.  Thus  mti.ximum  energ>’  coupling  is  aehiewd. 

iv _ FgjvriiTcntal  ln\  eviigati(nis 

The  energy  cffieiencs  of  p  h  o  t  ore  f  rac  1 1  \' e 
intereonnection  depends  on  several  ptirameters.  In  addition 
t(i  the  beam  o\erlap  mentioned  earlier,  it  depends  on  cryst.il 
(.orientation,  coupling  (.uiisiani.  lengtfi  of  interaclu'n.  beam 
iniensit)  ratio,  insertion  los:'  at  SlAl  11~.1S;,  and  (Mlier  loss 


mechanisms  such  as  Fresnel  reflection,  beam  fanning  [19] 
and  scattering,  etc.  We  have  carried  out  a  series  of 

experiments  to  investigate  the  energy  efficiency  of 
photorefractive  interconnection.  These  experiments  are 
designed  to  study  the  dependence  on  those  physical 
parameters  mentioned  above. 

Referring  to  Fig.  we  consider  the  cfficiencs  of 

energy  transfer  in  photorefractive  two-beam  coupling.  .An 
input  beam  of  power  Pj  is  split  by  a  beam  sp'itter  BS  into  a 

pump  beam  with  power  and  a  signal  beam  with  power  f\ 

(  P^  =  -  Pj,  ).  These  \\\o  beams  are  recombined  in  a 

photorefractive  BaTiO;  crystal.  I'he  crystal  is  oriented  such 

that  the  signal  beam  is  amplified.  Let  P^^  be  the  power  of  the 
amplified  signal  beam.  The  energy  efficients  t]  o  f 

photorefractive  two-beam  coupling  is  defined  as  ri  =  P^,  /  P^  . 
,*\  theoretical  expression  of  such  an  energy  cfficiencx  wa^ 


derived  earlier  in  reference  9.  In  our  experiment,  we 
examine  the  dependence  of  r|  on  the  pump-to-signal  intensity 
ratio  m,  defined  as  m  =  /  Ps  -  Such  a  parameter  can  be 

changed  by  using  a  variable  beam  splitter. 

Using  argon  ion  laser  as  the  source,  and  a  5  x  7  \  7 
mm-'  BaTiO;  crystal  as  the  photorefracti\'e  medium,  we  haNc 

measured  the  energ\’  efficiency  (n  )  as  a  function  of  the 
intensits'  ratio  m.  The  maximum  efficiency  is  approximate!} 
lU  when  the  beam  intensitx'  ratio  is  m  =  20.  In  addition  to 
these  measurements,  we  also  repeat  the  experiment  b\ 
inserting  a  neutral  density  filter  (NDF)  in  the  optical  path  of 
the  signal  beam.  Thi^  neutral  densit}  filter  is  used  to 
simulate  the  fanout  energy  loss.  The  idea  here  is  to  shoes  that 
the  amplified  signal  beam  power  does  not  decrease 
significantly  as  the  input  signal  beam  is  attenuated.  The 
experimental  results  are  shoevn  in  Fig. 4,  We  notice  that  an 
energ}  efficiencx  of  1.5^4  is  obtained  when  the  simulated 


fanout  loss  is  99.99c.  This  is  an  improvement  in  energy 
efficiency  by  more  than  a  factor  of  15.  As  discussed  earlier 
in  reference  9,  the  optimum  beam  splitting  ratio  m  (i.e.  the 
pump-to-signal  intensity  ratio  measured  immediately  after 
the  beam  splitter,  and  before  the  signal  beam  is  attenuated 

b>  the  NDF)  depends,  in  general.  on  the  photorefraciis  e 
coupling  constant  yL  and  the  fanout  loss  of  the  signal  beam. 
In  this  ptriicular  experiment.  the  crystal  is  oriented  such 
that  both  beams  enter  the  crystal  symmetrically  at  an  angle 

of  incidence  of  20'  as  shown  in  the  inset  in  Figure  4.  This 
configuration  is.  ho\se\er,  not  optimized  for  maximum 

energ)  efficiency  because  it  only  provides  a  coupling 

constant  yL  of  approximately  4.5. 

d  he  coupling  constant  yL  depends  on  the  cry  stal 
orientation  relative  to  the  beams.  In  another  experiment, 
the  same  crystal  is  oriented  such  that  the  signal  beam  enter 

th:  crystal  a:  an  angle  of  approximately  67'  uhich  is  near 
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Brewster's  angle.  This  incidence  configuration  minimizes  the 
Fresnel  reflection  loss.  The  pump  beam  enters  the  crystal  at 
an  angle  of  incidence  of  32°  such  that  the  grating  induced  in 
the  photorefractive  crystal  provides  a  coupling  constant  yL  of 
appro.ximately  9.2.  Again,  an  argon  ion  laser  is  used  as  the 
light  source.  The  beam  splitter  has  a  reflectance  of  95‘'.<' 
which  provides  a  beam  splitting  ratio  of  m  =  Fh  Lndcr  these 
conditions,  the  energy  efficiency  is  improved  signil'icaml\ , 
Figure  5  shows  the  e.xperimenta!  results  along  with  a 
theoretical  fit  (the  solid  line)  using  equation  (3)  in  reference 
9.  The  energy  efficiency  is  plotted  as  a  function  of  (1/N) 
which  is  the  transmittance  of  the  neutral  densits  filter.  We 
recall  that  (l/.N)  represents  the  maximum  energy  efficiency 
achievable  by  a  conventional  \  N  crossbar  as  a  result  of 

fanout  loss.  The  orientation  of  the  crystal  relative  to  the 
incoming  beatn  is  shown  in  the  inset  of  Figure  5.  We  notice 
that  an  improvement  by  a  factor  of  100  is  achieved  in  energs 
efficiency  w  hen  (1/N)  is  of  the  order  of  1  x  1(1' N  Further 
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improvement  in  energy  efficiency  can  be  achieved  by  coating 
the  crystal  surfaces  with  anti-reflection  dielectric  films  and 
by  using  crystals  with  higher  coupling  constants. 

In  the  experiments  discussed  above,  both  signal  and 
pump  beams  are  Gaussian  beams  with  no  spatial  intensit\ 
pattern.  In  the  interconnection  applications,  these  beams  are 
spatially  modulated.  As  we  mentioned  earlier,  the  energ\ 
efficiency  of  the  photorefractive  interconnection  as  depicted 
in  Fig.  1  depends  on  the  spatial  overlap  of  the  beams.  W'e 
novs  discuss  our  experimental  investigations  on  the 
utilization  of  Fourier  transform  to  achie\o  maximum  beam 
oserlap  and  thus  to  achieve  maximum  energy  efficienc\. 

The  experimental  configuration  is  illustrated  in  Fig. 6. 
The  output  from  an  argon  laser  (514. 5nm)  is  spatial-filtered 
and  expanded  to  form  a  collimated  beam  vsith  2cm  diameter. 
.'X  variable  beam  splitter  consisting  of  a  half-wave  plate  and 
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a  p'-larizing  beam  splitter  cube  is  used  to  split  the  incoming 
beam  into  two,  a  pump  beam  and  a  signal  beam.  Another 
half-wave  plate  is  used  to  rotate  the  polarization  of  the 
reflected  (signal)  beam  by  90°  from  the  direction 
perpendicular  to  the  plane  of  the  paper  to  the  in  plane 
direction.  .A  polarizer  placed  downstream  further  filters  out 
the  residual  perpendicular  component.  A  10x10  binarx 
matrix  mask  is  used  to  encode  a  spatial  pattern  op'o  the 
signal  beam.  The  objective  here  is  to  efficientlx  transfer 
energ)-.  via  photorefractive  txvo-xvave  mixing,  from  the  pump 
beam  to  the  signal  beam  xvhich  is  encoded  xxith  the 
interconnection  pattern.  .Anx-  one  or  more  (up  to  100)  of  the 
10x10  channels  can  be  selected  bx  using  an  appropriate 
mask. 

To  maximize  the  spatial  overlap  of  the  pump  and  the 
signal  hp;ims  inside  the  crystal  via  the  properties  of  Fourier 
transform  [see  equations  (1)  and  (2)|  described  in  Section  2. 
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a  pump  mask  with  a  single  aperture  identical  to  the  unit  cell 
of  the  signal  mask  is  used.  The  two  masks  are  placed  at  the 

front  focal  planes  of  tw'o  Fourier  transform  lenses  FLj  and 
FLt.  (of  identical  focal  length  f  =  25cm)  respectively.  Here  v\e 

use  two  lenses  of  identical  focal  length  because  optimum 
beam  coupling  occurs  at  an  angle  of  40“  between  the  beams 
Vvhich  are  beyound  the  numerical  aperture  of  a  single  lens. 
The  crystal  is  located  at  the  common  Fourier  plane  of  the  tuo 
masks  (see  Fig. 6).  The  shift-invariant  property  of  Fourier 
transformation  ensures  that,  apart  from  a  phase  factor,  the 
diffraction  pattern  from  each  signal  channel  oserlap'^  with 
that  from  the  pump  inside  the  crystal.  ’*  he  intensit)’  patterns 
of  the  pump  and  the  signal  at  the  image  plane  and  the 
Fourier  plane  are  shown  in  Fig. 7.  Note  that  the  tv.o  intensity 
distributions  at  the  Fourier  plane  differ  significantly  in  fine 
structures  due  to  multiple  beam  interference  among  the 
various  signal  channels.  A  magnified  version  that  re\eals 
these  fine  structures  can  be  found  in  Reference  [20i. 
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The  key  questions  we  address  in  this  section  are  the 
following;  (1)  How  efficiently  can  we  transfer  energ}'  from  a 
single  channel  to  multiple  channels  ?  (2)  How  does  the 
energy  efficiency  depend  on  the  number  of  signal  channels  ■? 
and  (3)  Hou’  is  the  energ\  distributed  among  the  signal 
channels  ? 

To  iru'estigate  how  the  energ>'  efficient'},  depends  on 
the  number  of  signal  channels,  we  used  a  xariahle 
rectangular  aperture  in  front  of  the  signal  mask  to  \'ar\  the 
number  cf  channels  from  lOO  to  1.  In  each  case,  ue  also 
adjusted  the  baser  power  and  the  \ariable  beam  splitter  so 
that  both  the  pump  power  and  the  total  signal  input  pouer 

are  fitscd  at  600  pW  and  6p\\’,  respectix  el} .  regardless  of  the 
number  of  signal  channels.  Figure  8  is  a  plot  of  energ} 
efficienc}  vs.  number  of  signal  channels.  The  energ} 
efficiency  turns  out  to  be  fairl}  insensitise  to  the  number  cd' 
signal  channels.  These  data  were  taken  with  the  geometr} 


shown  in  the  inset  of  Figure  4.  As  mentioned  before,  this 
particular  geometry  is  far  from  optimum,  and  one  can  expect 
an  energy  efficiency  of  30  to  409f  or  higher  once  optimized. 

With  all  tile  10x10  signal  channels  on.  the  intensity 
pattern  at  the  output  (image)  plane  is  sho\'>'n  in  Figure  ‘h  The 
ong.,.aI  output  pattern  with  the  photorefracti\ e  crxstal 
renKned  is  shown  in  Fig. 9  (a).  When  the  crNStal  is  in  place, 
the  output  pattern  is  slightly  degraded  as  shown  in  Fig, 9  (b). 
The  degradation  is  due  mainly  to  beam-fanning  [17],  This  can 
can  be  \erified  b>'  observing  the  output  patten,  using  a 
signal  Deam  which  is  ordiiianiy  poiarizeu  t o-polari/ed )  p; 
much  lower  in  power.  The  results  are  shown  in  Fig. 9  (c)  and 
id).  In  bold  casc^  the  criginad  intensity  nauern  is  restored  as 
the  beam  fanning  effect  diminished.  In  Fig. 9  (c).  significant 
reduction  of  beam  fanning  is  a  consequence  of  the  much 
smaller  electro-optic  coefficients  associate  with  the 
o-poL;rized  beam.  In  Fig. 9  id),  it  is  a  consequence  of  huiger 
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time  required  for  the  fanning  grating  to  build  up  uhen  the 
beam  intensity  is  weaker.  The  bottom  picture  in  F-ig.9(e) 
shows  the  amplified  signal  output  when  the  putnp  beam  i'^ 
turned  on.  The  presence  of  the  pump  beam  not  onl\ 
increases  the  output  energy  but  also  significant!)  reduce"  the 
signal  degradation.  To  monittm  the  cnerg)  d i ' m  1 1 o:. 
quantitalis  el) .  a  linear  detector  arra>  is  placed  along  one 
ro\s  of  the  signal  chatinels  at  the  image  plane,  I  lie 
oscil logratiis  representing  the  etiergx  distribution  among  ilie 
!('  channels  are  given  iti  Idg.K'  lor  the  lollowing  three  (.ascs: 
(a  I  vs  ith  the  crs  stal  retnc'sed;  (bi  vcith  the  cr_\stal  in  phuc 
and  the  pump  beam  c>ff;  ici  with  the  crvstal  in  place  and  the 
pump  beam  on.  In  all  cases,  the  f'laussian  '  top-hal  "  prolile 
originated  from  the  e.xpanded  laser  beam  is  more  or  less 
preser\  ed.  and  the  energy  distribution  is  su f  1 1 c  i e n tl  _\ 
uniform.  Note  that  the  vertical  scale  of  the  osciliogratii  in 
log.  10(b)  is  1 0m\7di\ ision  versus  the  1 0()m\7di\ isuMi  scale 
in  l  ig. 10(c).  d'he  signal  gain,  in  this  case  is  approximate!)  20. 


r^iiig  an  experimental  configuration  similar  to  the  one 


shovsn  in  Fig. 6.  uith  ail  the  Fourier  transform  lenses 
replaced  h\  astigmatic  optics  (which  image  along  the  x-a\is 
and  h'ouner  transform  along  the  y-axis  in  the  input  artn^.  and 
\ue  \ersa  in  the  output  arms),  ue  alsc)  demonstrated  a  hx'' 
geiiera!i/ed  ^  ro^-^har  sw  n  Ji ,  An  appropriate  pump  ma'-k  .md 
a  signal  ma'-k  for  a  specific  interconnection  pattern  are 
'tio'xn  in  ilie  leli  h.ind  side  r'f  I'lg.ll.  Tlie  cc'rre^l'iondin.: 
c;ensit\  [niiteMw  mside  tlte  cr\'-ta!  tit  the  Fourier  plane  are 
Aii'ssn  vUi  tlie  nglit.  IdealU,  xse  \xant  each  pum[^  eliatinel  ta 
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4.  SumtTiar\ 

In  sunimarx,  we  have  proposed  and  demonstrated  the 
use  ot  Fourier  transform  to  achieve  maximum  energ\ 
efficienc)  in  a  photorefracti ve  optical  interconnection.  We 
lui' e  shown  that  the  energ\  efficiencs  of  phoiorefraet  i  \  e 
dsnamie  tu'lograms  can  be  of  the  order  of  30^f  to  40^;,  In 
addition,  wc  ha\'e  also  \erified  that  the  energ\'  cfficienc\ 
in^ensiti\e  to  the  number  of  signal  channels,  and  that  the 
cnergx  distribution  among  the  signtil  channels  i^  suld’ieientl} 
n'otorm.  IdnalK  we  ha'.e  described  a  6x6  crossbar-NW  itch 
Using  this  approticl:, 

•Ml  the  experiments  described  in  this  paper  were 
carried  cmt  wiili  photographic  masks  replacing  the  SI. Ms.  do 
demonstrate  the  reconfigurability  and  to  studs  the 
reconI  iguration  time,  one  can  no  longer  avoid  the  use  of  SI. Ms 
which,  in  gcneriil,  ha\e  relativeU  poor  contrast  ratio  (of  the 


order  of  10  to  30)  relative  to  the  photographic  masks.  We 
are  currenth’  investigating  the  effect  of  limited  contrast  ratio 
of  SLMs  on  the  operation  and  the  reconfiguration  time.  The 
results  will  be  the  subject  of  future  publications. 
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HGURE  CAPTIONS 


Fig.l  A  schematic  diagram  illustrating  the  basic  idea  of  a 

reconfigurable  optical  interconnection  using  a 
photorefractive  crystal  in  conjunction  with  a  spatial  light 

modulator. 

Fig. 2  A  schematic  diagram  illustrating  the  basic  idea  of  using 

Fourier  transform  to  achieve  the  spatial  overlap  of  the  pump 
beam  and  the  signal  beam. 

Fig. 3  -An  experimental  configuration  for  measuring  the  energ) 

efficiencx’  in  photorefractive  two-beam  coupling. 

Fig. 4  The  energy  efficiency  (p  )  of  two-beam  coupling  in  a 

photorefractive  barium  titanate  crystal  as  a  function  of 
pump-to-signal  power  ratio  (m).  The  symbols  "  □ ’O". 
and  "A",  represent  the  cases  without  any  neutral  densits 
filter,  with  neutral  density  filters  NDFl.  NDF2.  and  NDF'.'. 
respective!) .  in  the  signal  input  arm.  The  percentage 

transmittance  of  the  three  neutral  densit)  filters  are  12.2''-;. 
0.749f,  and  0.129.  respectively.  The  orientation  of  the 
crxstal  relative  to  the  beam  is  given  in  the  inset. 

Fig. 5  Energy  efficiency  (ri)  of  photorefractive  two-beam  coupling 

in  a  barium  titanate  sample  as  a  function  of  the 

transmittance  of  a  neutral  density  filter  placed  in  the  signal 

input  arm.  The  transmittance  is  labeled  1/N  to  relate  it  to 
the  fanout  loss  of  a  NxN  permutation  crossbar  network.  The 


orientation  of  the  crystal  relative  to  the  beam  is  given  in  the 
inset. 

Fig. 6  An  e.xperimental  configuration  for  a  1-to-NxX  (for  N  =  10) 
broadcasting  network  using  photorefractive  holograms  at 
the  Fourier  domain. 

Fig. 7  The  intensity  patterns  of  the  masks  for  the  probe  and  the 

pump  beams  at  the  image  plane  and  the  Fourier  plane. 

Fig. 8  Energy  efficiencv  (ri  )  as  a  function  of  number  of  signal 

channels  (N)  in  a  1-to-NxN  broadcasting  configuration  using 
photorefractive  holograms  in  a  barium  titanate  sample. 

Fig.v  .A  iOxlO  intensity  pattern  at  the  output  image  plane  under 

sarious  conditions:  (a)  with  the  crystal  remoxed.  (b)  \siih 
the  crystal  in  place  and  pump  beam  off,  total  pouer  = 

40()p\V.  e-polarization,  (c)  same  as  in  (b)  but  with 
o-pol;irizaiion,  (d)  same  as  in  (b)  but  with  total  posser 
=  3p\\’,  (e)  amplified  output  signal,  total  signal  input  =  .'u\\'. 
pump  input  ==  bOOpW  . 

Fdg.lO  The  intensity  distribution  of  a  selected  row  of  i(>  out  td  the 
10x10  channels  shown  in  Fig. 9.  (at  with  the  crystal 
removed,  (b)  \sith  the  crystal  in  place  and  pump  beam  otf  . 
total  input  power  =10pW,  (c)  with  the  crystal  in  place  and 
pump  beam  on  .  input  pump  pow'er  =  400p\\.  The  vertical 
scale  (per  division)  in  the  pictures  are  5()m\'.  10m\’,  atid 
lOOmW  respecti\  el> . 


Fig.  11  The  intensity  patterns  of  the  masks  for  the  probe  and  the 
pump  beams  at  the  image  plane  and  the  crystal  plane  for  a 
6x6  generalized  crossbar  network. 
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Abstract 

Coherent  amplification  of  images  by  two  wave  mixing  in  photorefractive  crystals  is 
examined  with  attention  given  to  processing  in  the  Fourier  domain.  It  is  shown  that  the 
gain  that  is  experienced  as  the  probe  image  traverses  the  crystal  is  unifonn  across  the 
image.  TTie  gain  can  be  expressed  as  a  well  known  function  of  the  average  probe  to  pump 
intensity  ratio.  Experimental  verification  is  given  to  suppon  the  theory. 


The  phenomenon  of  energy  transfer  in  two  wave  mixing  in  nonlinear  optical  media 
has  been  used  effectively  for  the  coherent  amplification  of  images  [1.2].  In  diffusion 
driven  photorefractive  crystals,  for  example,  it  is  known  that  the  inherent  90°  phase  shift 
between  the  induced  refractive  index  grating  and  the  holographic  intensity  grating  is 
responsible  for  the  energy  coupling.  Shown  in  Fig.  1  is  a  typical  image  amplification 
experiment  in  which,  with  the  correct  geometry',  the  signal  image  wave  becomes  amplified 
at  'h?  “\-npn«;e  of  the  pump  wave  intensity  The  model  that  is  currently  being  usea  to 
describe  coherent  image  amplification  is  a  simple  one  developed  for  the  case  of  mi.xing  two 
plane  waves.  Since  the  signal  image  wave  can  be  expressed  as  a  superposition  of  many 
plane  waves  whose  propagation  directions  vary'  over  a  limited  range  (this  range  is  small 
when  compared  to  the  difference  between  the  average  signal  wave  direction  and  the  pump 
wave  direction),  the  simple  plane  wave  interaction  model  must  be  modified  to  be  more 
accurately  applicable.  In  this  letter,  we  present  a  simple  modification  of  Lhe  theory  of  two 
wave  mixing  involving  image  bearing  beams  and  discuss  in  particular  the  amplification 
process  when  the  crystal  is  placed  at  the  Fourier  transform,  plane  of  lens  LI .  Experimental 
evidence  is  also  given  to  support  the  theory. 

Shown  in  Fig.  2  is  the  image  amplification  geometry  in  which  the  amplitudes  of  the 

image  bearing  w'aves  are  represented  by  (Ej),  j=1.2 . N,  and  the  pump  plane  wave  is 

denoted  by  Ep.  In  photorefractive  crv'stals,  the  index  modulation  depends  on  the  intensity 
interference  pattern  so  that  the  present  analysis  begins  with  the  intensity  expression  which 
is  given  by 

1  r  N  N  N  T 

I(r)  —  I(t ,  1  +  T—  R e  Ep ^  E  j^exp[i(kp-km)’rj  +  ^  EqEj^exp[i(kq-km)  r]  ,,(1) 

^  m=l  q^m  m=l 

where  kpand  k^  are  the  wave  vectors  of  the  pump  and  the  mth  probe  beams  expressable, 
according  to  Fig.  2,  by 


kp  =  ^  ( -^sin9  +  £cos9  ) ,  kp,  =  ■^  (  ^sin9n,  +  fcos9m  )  ,  (2) 

^  /v  A 

N 

and  the  total  intensity  is  lo  =  !Epl-  +  X,  lEnJ-.  If  the  probe  waves  span  a  sufficiently  small 

m=l 

anguhtr  spectrum  (much  smaller  than  the  angle  between  the  pump  and  the  average  probe 
direction  such  that  N59  «  9),  then  the  interference  between  the  probe  waves  will  induce 
index  gratings  which  are  much  weaker  than  those  written  by  the  interference  between  the 
probe  waves  and  the  pump  wave  [3].  This  is  the  scenario  of  interest,  and  the  index 
distribution  can  then  be  approximated  as 

n  =  no-ReJi^  X  EpE^,  exp[-j(kp-kp,)-r]  1.  (3) 

^  m--; 


where  ni  is  the  index  modulation  coefficient  which  is  assumed  to  be  the  same  for  all  of  th.e 
gratings.  The  7t/2  phase  shift  is  due  to  the  diffusion  driven  process.  With  this  and  the 
assumption  of  slowly  vaio’ing  amplitudes  for  each  plane  wave  involved  in  the  interaction. 
Maxwell's  equations  yield  the  following  coupled  system: 


dEp 

dz 


IH,  V  IF  [2  F 


m  =  1 


ni>. 

4cos9 


dE,,, 

dz 


-i^lEJ-E„,,  m=l,: . N, 

1(1  '■ 


njX 

4cos9„, 


(4) 


The  standard  separation  of  intensity  and  phase  equations  is  useful,  and  in  particular,  the 
intensities  of  the  pump  and  probes,  Ip=IEpl^  and  Ini=lEmi2,  are  governed  by 


Ep  =  2Yp 


r  =  '>v 

*  m  iin- 


(5) 


In  the  assumed  diffusion  driven  scenario,  the  phases  of  die  beams  are  decoupled  so 
that  the  intensity  equations  describe  the  two  wave  mixing  process  completely.  To  obtain 
the  final  desired  result,  we  make  one  further  approximation  which  is  consistent  with  the 
assumptions  made  thus  far.  We  have  assumed  that  0m=  9,  so  that  it  is  justified  to 
approximate  ==  Fm-  With  this,  the  equations  become 


^=--7  I  I 
dz  1„  £ '»■ 

%=  f  m=l,2 . K. 


These  equations  are  easily  integrated  to  yield  the  following  solutions. 


R  hi  exnt-Fz) 

ln(Z,)  = - ^ ^ -  . 

1  +  P  exp(-Fz) 

h^/z) _ 1  +  P 

”  1  +  R  exp{-Fz) 
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The  above  is  the  main  result  of  this  letter.  The  solutions  indicate  that  the  probe  waves  are 
amplified  unifomily  with  a  gain  factor  which  is  a  function  of  the  ratio  of  the  pump  beam  to 
the  total  intensity  of  the  probes.  Small  deviations  from  this  are  expected  as  the  coupling 
coefficient  is  a  function  of  the  actual  probe  beam  angle,  and  also  as  harmonic  distortions  of 
the  gratings  result  from  the  modulation  depth  approaching  unity  [4-6]. 

The  uniformity  of  the  amplification  has  been  verified  by  the  following  experiment 
(see  Fig,  3).  An  output  beam  from  an  argon  laser  (514.5  nm)  is  spatially  filtered  and 
exptindcd  to  pass  through  a  rectangular  aperture  (RA:  2mm  x  6mm).  A  variable  beam 
.splitter  consisting  of  two  half-wave  plates  and  a  polarizing  beam  splitter  cube  is  used  to 
split  the  incoming  beam  into  two  (the  pump  and  'he  probe)  and  to  vary  the  intensity  ratio  of 
the  pump  and  the  probe  beams;  the  polarizer  positioned  after  the  half  wave  plate  in  the 


retlecied  path  ensures  that  both  transmitted  and  reflected  beams  have  the  same 
polari7ations.  The  probe  beam  is  transmitted  through  five  rectangular  windows  (().5mm  x 
2mm  each)  in  a  mask  (MAj)  to  form  5  probe  beamlets.  In  the  other  arm,  the  pump  beam 
illuminates  an  identical  rectangular  window  in  a  second  mask  (MA2)-  Two  spherical  lenses 
(FL]  and  FL2),  one  in  each  arm,  are  used  to  Fourier  transform  the  two  spatial  patterns  onto 
the  ciA'stal  plane.  A  spherical  lens  (FL3)  is  used  to  re-image  the  spatial  pattern  carried  by 
the  amplified  probe  beam  onto  the  detector  plane  where  the  optical  intensity  profile  is 
sampled  by  a  linear  detector  array  and  monitored  by  a  storage  oscilloscope.  We  have 
\  aried  the  pump-to-total  probe  power  ratio  in  the  range  1000  ~  1  and  obsers'ed  that  the 
energy  transferred  from  the  pump  is  equally  distnbuied  among  the  probe  beamlets  to  within 
2nT.  .-X  typical  example  of  the  intensity  profiles  is  shown  in  Fig.  4.  The  lower  trace  is  the 
intensity  profile  of  the  five  probe  beamlets  transmitted  through  the  crv  stal  when  the  pump  is 
off  The  upper  trace  represents  the  corresponding  profile  when  the  pump  beam  is  turned 
on.  For  this  specific  case,  the  optical  powers  of  the  pump  and  the  total  probe  are  0.6  mW 
and  0.2  mW.  respectively.  When  all  but  one  of  the  probe  beamlets  are  blocked,  almost  all 
of  the  depleted  pump  energy,  originally  distributed  among  the  five  channels,  redirects  itself 
into  the  single  unblocked  channel. 

TTe  dynamics  of  grating  adjustments  made  by  multiple  probe  beams  were  observed 
in  the  following  experiment.  A  BaTiOy  crystal  was  placed  at  the  intersection  of  two  weak 
probe  waves  and  a  strong  pump  beam  as  shown  in  Fig.  5.  As  shown,  the  two  probes  were 
arranged  so  that  they  propagate  with  a  small  angle  1  degree)  between  them  in 
comparison  with  the  angle  between  them  and  the  pump  beam  f~  30  degrees).  The  results 
are  shown  in  Figs.  6a-b.  The  input  intensities  of  the  two  probes  were  roughly  the  same 
fP-lOpW')  with  the  input  pump  power  set  at  P=5mW.  Shown  in  Fig.  6a  are  the  probe 
intensities  (after  passage  through  the  crystal  as  seen  by  detectors  1  and  2)  as  they  are 
strobed  on  with  shutters.  T'hc  amount  of  pump  power  that  is  scattered  into  the  probe 
detectors  was  small  compared  to  the  unamplified  probe  intensities  and  correspondingly 


becomes  negligible  when  compared  to  the  amplified  probe  beams.  Fig.  6b  shows  the 
amplified  probe  beams  as  the  probes  were  strobed  by  the  shutters  (top  trace=probel , 
bottom  trace=probe2;  scale  factor=l  V/mW).  The  power  shanng  is  clearly  evident  in  the 
mid  ponton  of  the  trace  where  both  probes  are  on.  When  the  order  of  events  w'ere  changed 
so  that  probe  2  was  turned  on  first,  the  results  w'ere  the  same  so  that  no  hysteresis  effects 
were  seen.  Note  that  probe  1  is  slightly  stronger  than  probe  2.  This  is  due  to  the  fact  that 
tile  two  beam  coupling  constants  (one  for  each  probe)  were  not  the  same. 

In  conclusion,  we  have  analyzed  the  nvo  wave  mixing  phenomenon  when  a 
multitude  of  probe  beams  are  used.  The  theoretical  results  shown  indicate  that  the  pump 
power  is  shared  by  the  probe  beams.  Experimental  evidence  was  given  to  validate  this 
power  shtinng  effect. 

We  acknowledge  useful  discussions  with  our  colleagues.  Fred  Vachss  and  Ragini 
Saxena.  Also,  at  the  time  of  the  submission  of  this  paper,  work  [7]  presented  at  the  OS  A 
annual  meeting  and  another  recently  published  |ki  which  dealt  with  similar  subjects  were 
brought  to  our  attention.  This  work  is  supported  primarily  by  the  Defense  .Advanced 
Research  .Aeenev. 
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Figure  Captions 

Typical  image  amplification  experiment  using  two  wave  mixing  in  photorefractive 
crystal'. 

Two  wave  mixing  with  the  ciw'stal  at  the  Fourier  plane  of  the  input  image.  The  pump 
Ep  is  a  plane  wave  and  the  image  effectively  consists  of  a  set  of  point  sources. 
Experimental  setup  used  to  verify  uniform  gain.  M.M  consists  of  five  identical 
apertures  and  comprises  the  probe  image  to  be  amplified.  Because  of  the  gaussian 
illumination,  the  intensity  pattern  at  M.Al  is  not  uniform.  M.A.2  is  the  pump  image 
apenure. 

Input  (lower  amplitude)  and  amplified  (higher  amplitude)  intensity  distributions.  Eac'^ 
pt'rtion  of  the  overall  input  distribution  receives  roughly  the  same  gain. 

Experimental  setup  used  to  record  the  dynamics  of  multiple  grating  buildup. 

[Dynamics  of  two  gratings  in  the  same  volume,  a)  input  probe  intensities  (upper 
corrc'ponds  to  “i  and  lo'wer-2  i.n  reference  to  Fig.  5).  b;  output  amplified  imeri'ine' 
il'istc'ry;  =1  turnei-l  on.  =1  turned  off,  #2  turned  on.  #2  turned  off.  -1  turned  on.  -2 
turned  on,  -'I  turned  ('ff. 
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ABSTRACT 


The  use  of  photorefractive  holograms  in  conjunction  with  a  spatial  light  modulator  (SLM)  to 
realize  a  reconfigurable  optical  interconnection  with  very  high  energy  efficiency  is  briefly 
discussed.  In  this  approach,  the  SLM  is  used  as  a  programmable  binary  matrix  mask  to  encode 
the  interconnection  pattern  to  a  coherent  laser  beam  whereas  the  photorefractive  crystal  is 
used  as  a  dynamic  holographic  medium  to  store  the  pattern  and  to  efficiently  diffract  the 
readout  beam  into  the  selected  channels  at  high  speed.  Vie  report  recent  experimental  results 
on  issues  such  as  energy  efficiency,  reconfiguration  time,  contrast,  and  uniformity. 

1.  INTRODUCTION 


Optical  interconnection  between  VLSI  circuits,  computing  chips  or  boards  plays  an 
important  role  in  parallel  optical  computing. ',2  Such  a  scheme  of  computing  provides  the 
potential  of  achieving  extremely  high  speed  because  it  uses  the  fast  switching  of  electronics  and 
the  wide  bandwidth  and  massive  parallelism  of  optics  for  communication. ^  A  generalized 
optical  crossbar  switch  can  provide  an  arbitrary  and  reconfigurable  interconnection'', ^  between 
an  array  of  N  lasers  and  an  array  of  N  detectors.  Conceptually,  such  an  interconnection  can  be 
achieved  by  using  optical  matrix-vector  multiplication,^-® 

V'zMV,  (J) 

where  V  is  the  input  vector  representing  the  signals  carried  by  an  array  of  N  lasers  and  V  is  the 
output  vector  representing  the  signals  received  by  the  array  of  detectors.  M  is  an  N  x  N  binary 
matrix  representing  the  interconnection  pattern.  Such  an  architecture  provides  the 
parallelism  required  for  arbitrary  interconnection.  Vl'hen  an  SLM  is  used  as  the  matrix  mask, 
the  interconnection  is  reconfigurable. 

The  parallelism  is  accompanied  with  an  intrinsic  fan-out  energy  loss;'',®  the  fan-out 
ieadb  to  an  energy  efficiency  of  1/N.  For  interconnection  with  large  N  (e.g.,  N  -  1000),  this 
large  amount  of  energy  loss  is  intolerable.  Vie  recently  proposed  and  demonstrated  a  novel 
concept  of  reconfigurable  optical  interconnection®, ’* 0  which  can  provide  very  high  energy 
efficiency.  Using  the  nonreciprocal  energy  coupling  in  photorefractive  crystals,''-*'*  such  a 
scheme  of  interconnection  is  capable  of  minimizing  the  fan-out  energy  loss;  thus,  it  achieves 
extremely  high  energy  efficiency.  Although  the  basic  concept  has  been  validated®,'®  using 
fixed  binary  matrix  masks  (or  transparencies)  and  the  energy  eHiciency  has  been  measured  and 
reported,'®,'®  the  demonstration  of  reconfigurability  with  high  energy  efficiency  and  high 
signal-to-noise-ratio  (S/N)  using  an  SLM  remains  a  subject  of  experimental  investigation.  This 
paper  reports  the  results  of  the  investigation  on  the  contrast  of  reconfigurable  optical 
interconnection  using  a  liquid  crystal  television  (LCTV)  in  conjunction  with  a  photorefractive 
barium  titanate  crystal.  The  effect  of  finite  contrast  of  the  SLM  on  the  performance  of  the 
holographic  interconnection  is  also  discussed. 

In  what  follows,  we  will  briefly  review  the  basic  principle  of  operation  of  the 
photorefractive  interconnection.  We  w'ill  summarize  some  of  our  earlier  results  on  energy 
efficiency  measurenient  and  on  the  use  of  Fourier  transforms  to  achieve  maximum  beam 
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overlap  in  the  phoiorefractive  medium. \^'e  then  discuss  the  results  of  our  experimental 
investigations  on  the  issue  of  contrast. 

2.  PRINCIPLE  OF  OPERATION 


Figure  1  shows  the  basic  idea  of  reconfigurable  optical  interconnection  using  photore- 
hologi uT.'. j.  A  bwain  ipliitci  is  aSCu  to  Split  off  a  small  portion  of  the  Deam  wnicn 
consists  of  beamlets  each  carrying  information  from  one  of  the  sources.  This  small  portion  of 
the  beam  will  be  called  the  signal  beam.  The  beam  splitter  allows  the  majority  of  the  energy  to 
pass  through.  This  major  portion  of  the  beam  is  called  the  pump  beam.  The  signal  beam  is 
fanned  by  the  cylindrical  lens  and  then  passes  through  the  spatial  light  modulator  (SLM)  which 
contains  the  interconnection  pattern.  After  passing  through  the  SLM,  the  signal  beam  is 
imprinted  spatially  with  the  interconnection  pattern.  The  pump  beam  does  not  pass  through  the 
spatial  ligt^r  modulator  and  thus  suffers  no  energy  loss  due  to  the  fan-out  and  all  <^ther  loss 
mechanisms  associated  with  the  SLM. 


V  Vv 


Fig.  1  A  schematic  d.ag'am  illustrating  the  basic  idea  of  a  reconfigurable  optical 
interconnection  using  a  photoref  ractive  crystal  in  conjunction  with  a  spatial  light  modulator. 

The  signal  beam  which  contains  the  interconnection  information  is  then  recombined  with 
the  pump  bean,  in  a  photorefractive  crystal.  Under  the  appropriate  conditions,  the  signal  beam 
can  be  amplified  at  the  expense  of  the  pump  beam.  Most  of  the  energy  of  the  pump  beam  will  be 
transferred  to  the  signal  beam  provided  the  length  of  interaction  is  large  enough.  ^  j,  >  “• 

Although  photorefractive  tw'o-wave  mixing  is  involved,  such  a  scheme  can  also  be 
considered  as  a  real-time  holographic  interconnection  with  a  very  high  diffraction  efficierxry. 
The  combination  of  the  beam  splitter  and  the  SLM  in  conjunction  with  the  pump  beam  is  used  to 
record  a  hologram  inside  the  photorefractive  crystal.  Such  a  hologram  contains  the 
interconnection  as  prescribed  by  the  SLM.  Once  the  hologram  is  formed,  the  pump  beam  can  be 
diffracted  off  the  hologram  and  be  redirected  into  the  array  of  detectors.  The  advantage  here 
is  that  the  recording  and  the  readout  occur  simultaneously.  This  offers  the  possibility  of 
reconfigurable  interconnection  by  controlling  the  SLM. 

To  achieve  maximum  energy  efficiency,  it  is  important  that  the  signal  beam  and  the  pump 
bearri  overlap  completely.  Specifically,  beamlet  1  of  the  pump  beam  in  Fig.  1  must  overlap 
completely  with  column  1  of  the  signal  beam.  And  similarly,  beamlets  2,  3,  and  4  must  also 
T  overlap  completely  with  their  corresponding  columns  in  the  signal  beam.  Although  the  beamlets 

in  the  signal  beam  and  the  pump  beam  are  intrinsically  different  because  of  the  interconnection 
pattern,  the  urdi'.iduai  pixels  can  have  loentical  shape  (e.g.,  square).  Complete  overlap  is 
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possible  in  the  Fourier  domain,  provided  all  the  pixels  are  identical.  The  results  of  our 
experimental  investigations  of  beam  overlap  and  the  energy  efficiency  were  presented 
recently'^  and  discussed  in  detail  in  a  separate  paper,  Jhe  key  results  pertinent  to  the 
energy  efficiency  are  summarized  in  the  next  section. 

3.0  ENERGY  EFFICIENCY 


f 


# 


'^e  first  investigated  the  energy  efficiency  due  to  beam  coupling  in  a  photorefractive 
crystal.  Referring  to  Fig.  2,  we  define  the  efficiency  <h)  of  energy  transfer  in  photorefractive 
holograms  as  n  ^  optical  power  of  the  input  beam  and  the 

amplified  signal  beam,  respectively.  A  theoretical  expression  of  such  an  energy  efficiency  was 
derived  earlier  in  reference  9  Using  an  uncoated  barium  titanate  crystal  (-  5  «  7  .  7  mm, 
30‘’-cut,  from  Sanders  Associates)  as  rhe  photorefractive  medium  and  an  argon  ion  laser 
(514.3  nm)  as  the  source,  we  have  achieved  an  energy  efficiency  as  high  as  30%.  Further 
improvement  in  energy  efficiency  can  be  achieved  by  coating  the  crystal  surfaces  with 
antireflection  dielectric  films  and  by  using  crystals  with  highei  crA'pl...^.  eoiiiidiiti.  Ultimately, 
the  energy  efficiency  is  limited  by  the  bulk  absorption  in  the  crystal. 


INPUT  BEAM 

f 


^  PUMP  BEAM 


AMPLIFIED 
SIGNAL  OUTPUT 

ENERGY  EFFICIENCY  (r;)  . 

PROBE-TO-PUMP  POWER  RATIO  (r)  ■ 


Fig.  2  An  experimentai  configuration  for  measuring  the  energy  efficiency  in  photorefractive 
two-beam  coupling. 

Using  an  experimental  configuration  as  schematically  illustrated  in  Fig.  3,  we  have  also 
successfully  demonstrated  ‘  ‘  ^  that  : 

•  the  energy  efficiency  is  insensitive  to  the  lateral  position  of  the  signal  w'indow  in  the 
SUM  (shift-invariant). 

•  the  energy  efficiency  is  insensitive  to  the  number  of  signal  channels 

•  the  energy  distribution  among  all  the  signal  channels  is  fairly  uniform. 

The  properties  listed  above  are  due  to  matched  amplification  at  the  Fourier  plane. As 
indicated  in  Fig.  3,  each  individual  pixel  of  the  signal  beam  is  identical  to  that  of  the  pump 
beam  to  ensure  max. mum  beam  overlap  at  the  Fourier  plane.  Another  desirable  characteristic 
of  two-wave  mixing  at  the  Fourier  plane  is  the  preservation  of  contrast  of  the  input  pattern.  In 
the  next  section,  we  discuss  the  issue  of  contrast  and  a  simple  method  to  enhance  the  contrast. 
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Fig.  3  An  expenriiental  configuration  for  a  1-to-N  »  N'  (for  N  =  10)  brcaF'-^sTing  nPT.u'orW 
photorefractive  holograms  at  the  Fourier  domain. 

4.0  CONTRAST 


All  of  our  experimental  results  reported  earlier  and  summarized  in  the  previous  section  are 
t  obtained  with  photographic  masks  replacing  the  SLMs.  To  demonstrate  the  reconfigurability 

and  to  study  the  reconfiguration  time,  an  SLM  in  which  the  interconneuliori  Oaiiern  caii  be 
varied  either  electrically  or  optically  can  be  used.  In  our  experiment,  we  used  liquid  crystai 
television  (LCT\)  which  has  a  relatively  poor  contrast  ratio  (of  the  order  of  -  to  3C)‘N‘8 
compared  with  photographic  masks.  In  the  paraHel  matrix-vector  multiplication  approach,  the 
maximum  number  of  fan-out  of  the  optical  cross-bar  can  be  limited  bv  the  finite  contrast  of  the 
SLM. 

As  a  sim.ple  examiple,  consider  an  M  »  M  permutation  matrix  mask  with  contrast  ratio 
C  =  Tj/Tq,  where  T;  and  Tf,  are  the  iraensity  transmittance  of  the  "ON"  state  and  the  "OFF" 
state,  respectively.  Referring  to  Fig.  4  (illustrated  for  the  case  of  M  =  4),  let  us  assume  that 
the  signal  from  each  source  is  intensity  modulated  between  zero  and  peak  value  P.  The  signal 
level  received  by  each  detector  through  the  "ON"  cell  is  PT),  whereas  the  noise  level  received 
through  the  (M-1)  "OFF"  cell  is  P(M-l)  Tq.  The  signal  to  noise  ratio  S/N  is  therefore  given  by 

S/iN  =  Tj/(M-I)  Tq  ^  C/(M-I)  =  C/M  for  M  »  1.  (2) 

The  condition  S/N  >  1  is  satisfied  provided  t.iat  the  number  of  fan-out  channels  (M)  is  less  than 
the  contrast  ratio  (C). 

t  A  simple  method  to  improve  the  contrast  is  by  double-passage  through  the  same  SLM  via 

retroreflection  or  phase  conjugation.  In  principle,  if  the  contrast  ratio  for  single  passage  is  C, 
the  contrast  ratio  for  double  passage  will  become  C2.  The  experimental  configuration  used  to 
^  verify  this  fact  is  shown  schematically  in  Fig.  5.  An  expanded  and  collimated  laser  beam  from 

an  argon  laser  (314.3  nni)  is  transmiitted  through  a  circular  aperture  (diameter  ■=  14  mm)  located 
at  the  input  plane.  Half  of  the  aperture  is  masked  by  a  neutral  density  filter  (NDF)  to  form  a 
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2  D  SPATIAL  LIGHT  MODULATOR 


Fig. 4  Reconf igurabk  interconnection  using  an  SLM  with  finite  contrast. 


Fig.  5  An  experimental  configuration  for  contrast-enhancement  by  double-passage. 

binary  amplitude  mask  at  the  input  plane.  The  transmitted  beam  is  collected  by  a  lens  (L)  and 
directed  into  a  barium  titanate  crystal.  The  crystal  is  oriented  so  that  the  incoming  beam  is 
retro-rellected  via  self-pumped  phase  conjugation.  ‘  ^  After  traversing  through  the  aperture, 
the  phase  conjugated  beani  is  sampled  by  a  beam  splitter  through  imaging  optics  which  reimages 
the  binary  mask  on  the  output  plane.  The  input  intensity  distribution  (i.e.  single-passage 
through  the  mask)  and  the  output  intensity  distribution  (i.e.,  double-passage  through  the  i:.ask) 
are  measured  by  scanning  a  detector  with  a  small  aperture  (diameter  e  Imm'  across  the  beam 
diameter  at  the  input  and  the  output  planes,  respectively.  The  experimental  .esults(for  a  mask 
with  optical  density  OD  =  0.5)  are  shown  in  Fig.  6(a)  and  (b).  Apart  from  some  imperfection  at 
the  edge  of  the  output  image,  the  output  contra'*  ratio  is  approximately  the  square  of  that  of 
the  input  as  expected.  Instead  of  a  phase-conjugate  mirror,  a  plane  mirror  can  also  be  placed  at 
the  input  plane  immediately  after  the  mask  to  retro-reflect  the  beam  back  through  the  binary 
i  mask.  The  resulting  intensity  distribution  is  shown  in  Fig.  6(c).  e  have  repeated  the  same 

measurement  with  a  different  NDF  (OD  =  1.0)  and  also  with  a  LCTV  at  the  input  plane;  the 
results  are  tabulated  in  Table  1.  Photographs  of  the  output  image  of  an  arbitrary  binary  pattern 
written  on  a  LCTV  are  shown  in  Fig.  7(a)  for  mirror  reflection  and  (b)  for  seif-pumped  phase 
conjugation.  The  phase  conjugation  has  the  key  advantage  of  being  intrinsically  self-aligned. 
The  grating  formation  time  for  self-pumped  phase  conjugation  using  a  typical  barium  titanate 
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I'lg.b  Intensii;.  pi-oiiici  arid  cor.trasi  o:  a  binary  iiiiage:  (a)  ai  the  input  plane,  (b)  at  the  output 
plane  after  nouD.e-Das>age  '.la  phase-conjugate  reflection,  (c)  at  the  output  plane  after  double¬ 
passage  v:a  nurror  reflection.  The  arrov^s  indicate  the  intensities  (of  the  bright  and  daik 
regions)  used  for  uaiculating  the  contrast  ratio. 
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C>.rtrast  Lnnaricenient  of  Binary  Images  by  Double-Passage 
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crystal  a;  a  fcvs  teris  of  milliwatt  of  optical  pov.er  is,  however,  relatively  slov.  (  a  fraction  of  a 
seco’  rj  to  several  seconds).  Note  that  the  output  obtained  by  mirror  reflection  is  fairly  uniform 
(see  Fig.fafc)  and  Fog. 7(a))  and  the  reflection  is  practically  instantaneous. 

The  pn^torefractive  iii terconnect ion  holograms  used  to  diffract  the  pump  Dearn  into  the 
sigr.ai  c’idrcie.5  specified  bv  the  mat:  .>.  mask  do  not  degrade  the  contrast  provided  that  the 
T  ive  crystal  is  (tlai  ed  at  the  f-ourier  plane. 'S'  ‘  In  contrast,  photorelractive  image 
amp! ;  f  Km.  L  ion  via  tv-o-uave  rmxing  in  the  iniage  plane  tends  to  decrease  the  t  ontrast  of  the 
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Fig. 7  Double-passage  image  of  a  binary  matrix  pattern  on  a  liquid  crystal  T.V.  (a)  via  mirror 
reflection,  (b)  via  phase-conjugate  reflection. 

input  image. -  An  expenmenta!  result  varifying  this  contrast  preserving  property  is  given  ;n 
Fig.  8.  Intensits  profiles  of  a  binary  image  and  its  amplified  output  are  shov.’n  in  (a)  and  (b), 
respectiveK .  Note  that  the  contrast  ratio  of  the  input  and  the  output  image  are  almost 
identical  and  that  the  verticai  scale  in  (b)  is  20  times  that  of  (a). 

SC5001 7 


(a) 


(b) 


Fig. 8  Experimental  results  showing  contrast-preserving  image  amplif ication  by  photorefractive 
two-wave  mixing  at  the  Fourier  domain;  (a)  input  intensity  profile,  (b)  amplified  output  intensity 
profile. 

To  extrapolate  froni  the  above  discussion,  if  an  SLM  with  a  contrast  ratio  of  the  order  of 
100  (for  example,  a  ferroelectric  liquid  crystal  SLM)  is  used  in  the  double  passage  mode,  the 
maximum  number  of  channels  (M)  of  a  M  »  M  crossbar  switch  limited  by  the  contrast  ratio  of 
the  SLM  can  be  as  high  as  lO**! 

5.0  EXPERIMENTAL  DEMONSTRATION  OF  RECONFIGUR ABILITY' 


Using  an  experimental  configuration  similar  to  the  one  illustrated  in  Fig.  3,  with  the 
photograpnic  mask  (transparency)  '■eplaced  by  an  SLM,  we  have  demonstrated  the  reconfigur¬ 
ability  of  the  interconnection  by  writing  new  interconnection  patterns  on  the  LCTV.  The 
patterns  are  generated  by  an  IBM  PC  and  sent  to  the  LCTV.  With  optical  power  of  the  order  of 
0  tens  of  milliwatts,  the  franie  rate  is  currently  limited  by  the  photorefractive  response  time  to  a 
few  frames  per  second.  The  data  rate,  however,  is  not  limited  by  the  photorefractive  response 
time  and  can  be  higher  than  several  megahertz.^  To  enhance  the  contrast  of  the  binary  pattern, 
^  the  signal  beam  is  passed  through  the  SLM  twice  via  a  retroref lecting  mirror  (as  described  in 
Section  4.0)  prior  to  interacting  with  the  pump  beam  in  a  photorefractive  barium  titanate 
crystal. 
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6.  SUMMARY 

We  have  described  a  novel  method  for  reconfigurable  interconnections  using 
photorefractive  holograms  in  conjunction  with  an  SLM.  Several  important  features,  including 
high  energy  efficiency  (>  30%),  uniform  energy  distribution,  and  high  data  rate  transmission 
(>  several  megahertz)  of  this  dynamic  holographic  interconnection  are  discussed.  We  have  also 
described  the  problem  associated  with  the  finite  contrast  of  the  SLM  and  have  shown  that  the 
contrast  can  be  greatly  enhanced  (i.e.,  scpjared)  by  double  passage  through  the  SLM  either  by  a 
retro-reflecting  mirror  or  by  a  phase-conjugator.  We  also  report  the  demonstration  of 
reconfigurability  by  using  a  LCTV  in  conjunction  with  a  photorefractive  barium  titanate  crystal. 
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